CERN-TH/2003-307 



The Magnetized Universe 



Massimo Giovannini^ 



CERN, Theoretical Physics Division, CH-1211 Geneva 23, Switzerland 

Abstract 

Cosmology, high-energy physics and astrophysics are converging 
on the study of large-scale magnetic fields. While the experimen- 
tal evidence for the existence of large-scale magnetization in galaxies, 
clusters and superclusters is rather compelling, the origin of the phe- 
nomenon remains puzzling especially in light of the most recent obser- 
vations. The purpose of the present review is to describe the physical 
motivations and some of the open theoretical problems related to the 
existence of large-scale magnetic fields. 
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1 A triple point 



Why do we observe magnetic fields in the sky? Why do we hve in a magne- 
tized Universe? A variety of observations imply that stars, planets, galaxies, 
clusters of galaxies are all magnetized. The typical magnetic field strengths^ 
range from few n G (in the case of galaxies and galaxy clusters) to few G 
(in the case of planets, like the earth) up to 10^^ G (in the case of neutron 
stars). Physical phenomena are characterized depending upon the typical 
time and length scales where they take place. Magnetic fields of stars and 
planets are related to length-scales which are much smaller than the diameter 
of the Milky Way (of the order of 30 Kpc) or of the local superchistcr. In 
this sense, a magnetic field of the order of G is minute on the terrestrial 
scale but it is sizable over the scale of the supercluster. 

In the present review only large-scale magnetic fields will be consid- 
ered, i.e. magnetic fields whose typical scale exceeds the AU ( lAU = 
1.49 X lO^^cm). Large-scale magnetic fields must be understood and treated 
as an essential part of the largest structures observed in the sky. Legitimate 
questions arising in the study of magnetized structures concern their forma- 
tion and evolution. The answers to these questions are still not completely 
settled. Specific observations will eventually allow to discriminate between 
the different competing explanations. 

Magnetic fields of distant spiral galaxies are in the fi G range. There 
is also compelling evidence of large-scale magnetic fields which are not as- 
sociated with individual galaxies. This empirical coincidence reminds a bit 
of one of the motivations of the standard hot big-bang model, namely the 
observation that the light elements are equally abundant in rather different 
parts of our Universe. The approximate equality of the abundances implies 
that, unlike the heavier elements, the light elements have primordial origin. 
The four light isotopes D, ^He, ^He and ^Li are mainly produced at a specific 
stage of the hot big bang model named nucleosynthesis occurring below the 
a typical temperature of 0.8 MeV when neutrinos decouple from the plasma 
and the neutron abundance evolves via free neutron decay. The abundances 
calculated in the simplest big-bang nucleosythesis model agree fairly well 
with the astronomical observations. In similar terms it is plausible to argue 

^In this review magnetic fields will be expressed in gauss. In the SI units IT = lO^G. 
For practical reasons, in cosmic ray physics and in cosmology it is also useful to express 
the magnetic field in GeV^ (in units ?i = c = 1). Recalling that the Bohr magneton is 
about 5.7 X IQ-^MeV/T the conversion factor will then be IG = 1.95 x IQ-^OCeV^. 
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that large-scale magnetic fields have similar strengths at large scales because 
the initial conditions for their evolutions were the same, for instance at the 
time of the gravitational collapse of the protogalaxy. The way the initial con- 
ditions for the evolution of large-scale magnetic fields are set is generically 
named magnetogenesis. 

There is another comparison which might be useful. Back in the seven- 
ties the so-called Harrison-Zeldovich spectrum was postulated. Later, with 
the developments of inflationary cosmology the origin of a flat spectrum of 
curvature and density profiles has been justified on the basis of a period of 
quasi-de Sitter expansion named inflation. It is plausible that in some infla- 
tionary models not only the fluctuations of the geometry are amphfled but 
also the fluctuations of the gauge flelds. This happens if, for instance, gauge 
couplings are effectively dynamical. As the Harrison-Zeldovich spectrum can 
be used as initial condition for the subsequent Newtonian evolution, the pri- 
mordial spectrum of the gauge flelds can be used as initial condition for 
the subsequent magnetohydrodynamical (MHD) evolution which may lead, 
eventually, to the observed large-scale magnetic flelds. 

Cosmologists and theoretical physicists would like to understand large- 
scale magnetization in terms of symmetries which are broken. There are 
other two different motivations leading, independently, to the study of large- 
scale magnetic flelds. For instance it was observed long ago by Fermi that 
if cosmic rays are in equilibrium with the galaxy, their pressure density is 
comparable with the one of the magnetic fleld of the galaxy. Large-scale 
magnetic flelds are also extremely relevant for astrophysics. Magnetic flelds 
in galaxies are sufficiently intense to affect the dynamics of interstellar gas 
both on the galactic scale and also on smaller scales characteristic of star 
formation processes. This is the third motivation leading to the study of large 
scale-magnetic flelds: the astrophysical motivation which tries to combine 
our knowledge of the Universe in powerful dynamical principles based on the 
microscopic laws of nature. 

The cross-disciplinary character of the physical phenomena addressed in 
this review is apparent from the table of content. In Section 2 a brief his- 
torical account of this flfty years old problem will be given. Then, to clarify 
the nature of the various theoretical constructions, the empirical evidence of 
large scale magnetic flelds will be discussed in Section 3. Section 4 contains 
some background material on the evolution of magnetic flelds in globally 
neutral plasmas in flat space. Section 5 and 6 address the problem of the 
evolution of large scale magnetic flelds in curved backgrounds and of their 
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origin. In Section 7 the possible effects of magnetic fields on the thermody- 
namical history of the Universe are scrutinized. Section 8 will be concerned 
with the effects of large scale magnetic fields on the Cosmic Microwave Back- 
ground radiation (CMB) and on the relic background of gravitational waves. 
In Section 9 the role of gravitating magnetic fields in cosmological solutions 
will be swiftly pointed out. 

The perspective of present review is theoretical. In this sense various (very 
important) experimental results (for instance concerning optical polarization) 
will receive only a swift attention. 

2 From Alfven to ROSAT 

In January 1949 large-scale magnetic fields had no empirical evidence. The 
theoretical situation can be summarized as follows: 

• the seminal contributions of Alfven P convinced the community that 
magnetic fields can have a very large life-time in a highly conducting 
plasma; 

• using the new discoveries of Alfven, Fermi j2] postulated the existence 
of a large-scale magnetic field permeating the galaxy with approximate 
intensity of /i G and, hence, in equilibrium with the cosmic rays; 

• rather surprisingly Alfven [Hj did not appreciate the implications of the 
Fermi idea and was led to conclude (incorrectly) that cosmic rays are 
in equilibrium with stars disregarding completely the possibility of a 
galactic magnetic field; 

• in 1949 Hiltner jl] and, independently. Hall 'P\ observed polarization 
of starlight which was later on interpreted by Davis and Greenstein jU] 
as an effect of galactic magnetic field aligning the dust grains ^. 

According to the presented chain of events it is legitimate to conclude that 

''it should be noticed that the observations of Hiltner ^ and Hall j^l took place from 
November 1948 to January 1949. The paper of Fermi "5' was submitted in January 1949 
but it contains no reference to the work of Hiltner and Hall. This indicates the Fermi was 
probably not aware of these optical measurements. 
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• the discoveries of Alfven were essential in the Fermi proposal (who was 
already thinking of the origin of cosmic rays in 1938 before leaving 
Italy); 

• the idea that cosmic rays are in equilibrium with the galactic magnetic 
fields (and hence that the galaxy possess a magnetic field) was essential 
in the correct interpretation of the first, fragile, optical evidence of 
galactic magnetization. 

The origin of the galactic magnetization, according to [2j, had to be somehow 
primordial. This idea was further stressed in two subsequent investigations 
of Fermi and Chandrasekar [71 |H] who tried, rather ambitiously, to connect 
the magnetic field of the galaxy to its angular momentum. 

In the fifties various observations on polarization of Crab nebula sug- 
gested that the Milky Way is not the only magnetized structure in the sky. 
The effective new twist in the observations of large-scale magnetic fields 
was the development (through the fifties and sixties) of radio-astronomical 
techniques. From these measurements, the first unambiguous evidence of 
radio-polarization from the Milky Way (MW) was obtained [TU] . 

It was also soon realized that the radio-Zeeman effect (counterpart of 
the optical Zeeman splitting employed to determine the magnetic field of 
the sun) could offer accurate determination of (locally very strong) magnetic 
fields in the galaxy. The observation of Lyne and Smith ^T] that pulsars 
could be used to determine the column density of electrons along the line 
of sight opened the possibility of using not only synchrotron emission as a 
diagnostic of the presence of a large-scale magnetic field, but also Faraday 
rotation. In the seventies all the basic experimental tools for the analysis of 
galactic and extra-galactic magnetic fields were present. Around this epoch 
also extensive reviews on the experimental endeavors started appearing and 
a very nice account could be found, for instance, in the review of Heiles ^2] • 

It became gradually clear in the early eighties, that measurements of 
large-scale magnetic fields in the MW and in the external galaxies are two 
complementary aspects of the same problem. While MW studies can provide 
valuable informations concerning the local structure of the galactic magnetic 
field, the observation of external galaxies provides the only viable tool for the 
reconstruction of the global features of the galactic magnetic fields. As it 
will be clarified in the following Sections, the complementary nature of global 
and local morphological features of large-scale magnetization may become, 
sometimes, a source of confusion. 
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Since the early seventies, some relevant attention has been paid not only 
to the magnetic fields of the galaxies but also to the magnetic fields of the 
clusters. A cluster is a gravitationally bound system of galaxies. The local 
group (i.e. our cluster containing the MW, Andromeda together with other 
fifty galaxies) is an irregular cluster in the sense that it contains fewer galaxies 
than typical clusters in the Universe. Other clusters (like Coma, Virgo) are 
more typical and are then called regular or Abcll clusters. As an order of 
magnitude estimate, Abell clusters can contain 10^ galaxies. 

In the nineties magnetic fields have been measured in single Abell clus- 
ters but around the turn of the century these estimates became more reliable 
thanks to improved experimental techniques. In order to estimate magnetic 
fields in clusters, an independent knowledge of the electron density along 
the line of sight is needed (see Sec. 3). Recently Faraday rotation mea- 
surements obtained by radio telescopes (like VLA have been combined 
with independent measurements of the electron density in the intra-cluster 
medium. This was made possible by the maps of the x-ray sky obtained with 
satellites measurements (in particular ROSAT ^). This improvement in the 
experimental capabilities seems to have partially settled the issue confirming 
the measurements of the early nineties and implying that also clusters are 
endowed with a magnetic field of /xG strength which is not associated with 
individual galaxies. 

The years to come are full of interesting experimental questions to be 
answered. These questions, as it will be discussed in the following, may have 
a rather important theoretical impact both on the theory of processes taking 
place in the local universe (like ultra-high energy cosmic ray propagation) and 
on the models trying to explain the origin of large scale magnetic fields. Last 
but not least, these measurements may have an impact on the physics of CMB 
anisotropics and, in particular, on the CMB polarization. In fact, the same 
mechanism leading to the Faraday rotation in the radio leads to a Faraday 
rotation of the CMB provided the CMB is linearly polarized. One of the 
important questions to be answered is, for instance, the nature and strength 
of the supercluster magnetic field and now more careful statistical studies are 
starting also along this important direction. Superclusters are gravitationally 
bound systems of clusters. An example is the local supercluster formed 

^The Very Large Array Telescope, consists of 27 parabolic antennas spread over a 
surface of 20 km^ in Socorro (New Mexico) 

^Thc ROcgtcn SATcUite (flying from June 1991 to February 1999) provided maps of the 
x-ray sky in the range 0.1-2.5 keV. A catalog of x-ray bright Abell clusters was compiled. 
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by the local group and by the VIRGO cluster. Together ROSAT, various 
observations with the EINSTEIN, EXOSAT, and GINGA satelhtes showed 
the presence of hot diffuse gas in the Local Supercluster. The estimated 
magnetic field in this system is, again, of the order of the n G but the 
observational evidence is still not conclusive. Another puzzling evidence is 
the fact that Lyman-a systems (with red-shifts z ~ 2.5) are endowed with a 
magnetic field [TB] . 

From the historical development of the subject a series of questions arises 
naturally: 

• what is the origin of large-scale magnetic fields? 

• are magnetic fields primordial as assumed by Fermi more than 50 years 



• even assuming that large-scale magnetic fields are primordial, is there 
a theory for their generation? 

3 Large-scale magnetic fields observations 

Before describing the theory of something (not of everything) it is highly 
desirable, in natural sciences, to understand in some detail the empirical evi- 
dence of the subject under discussion. While technical accounts of the various 
experimental techniques exist already [121 UHl UH Uil UHl UZl UHI it seems use- 
ful, in the present context, to give an account of the experimental situation 
by emphasizing the physical principles of the various measurements. There 
are valuable books dealing directly with large-scale magnetism |19| QUI 1^ . 
Furthermore excellent reviews of the morphological features of large scale 
magnetism can be found in El ■ Even if the present review will 

not be concerned with planetary magnetic fields, analogies with physical sit- 
uations arising closer to the earth are often useful and, along this perspective, 
magnetic field observations in the solar system have been recently reviewed 
[22\ . 

Magnetic fields observed in the Universe have a homogeneous (or uniform) 
component and a non-uniform component. It is then useful, for the purposes 
of this Section, to write, in a schematic notation, that 
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where B is the uniform component and SB the non- uniform component. 
Both components are phenomenologically very relevant. Different experi- 
mental techniques can probe different components of large-scale magnetic 
fields (i.e. either the total magnetic field or its homogeneous component). 
Another important distinction is between B± and B\\: there are measure- 
ments (like synchrotron emission) which are sensitive to B± i.e. the trans- 
verse component of the magnetic field; there are other measurements (like the 
Faraday Rotation measure or the Zeeman splitting of spectral lines) which 
are sensitive to B\\ i.e. the magnetic field along the line of sight. 



3.1 Zeeman splitting 

In order to measure large scale magnetic fields, one of the first effects coming 
to mind, is the Zeeman effect. The energy levels of an hydrogen atom in 
the background of a magnetic field are not degenerate. The presence of a 
magnetic field produces a well known splitting of the spectral lines: 

Auz = (3.2) 

/vrme 

where B\^ denotes the uniform component of the magnetic field along the line 
of sight. From the estimate of the splitting, the magnetic field intensity can 
be deduced. Indeed this technique is the one commonly employed in order 
to measure the magnetic field of the sun [221 • The most common element in 
the interstellar medium is neutral hydrogen, emitting the celebrated 21-cm 
line (corresponding to a frequency of 1420 MHz). If a magnetic field of /i G 
strength is present in the interstellar medium, according to Eq. ()3.2|) . an in- 
duced splitting, Auz ~ 3Hz, can be estimated. Zeeman splitting of the 21-cm 
line generates two oppositely circular polarized spectral lines whose apparent 
splitting is however sub-leading if compared to the Doppler broadening. In 
fact, the atoms and molecules in the interstellar medium are subjected to 
thermal motion and the amount of induced Doppler broadening is roughly 
given by 



where fth is the thermal velocity oc yT/m where m is the mass of the atom 
or molecule. The amount of Doppler broadening is Ai/^op ~ 30 kHz which is 
much larger than the Zeeman splitting which should be eventually detected. 
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Zeeman splitting of molecules and recombination lines should however be 
detectable if the magnetic field strength gets larger with the density. Indeed 
in the interstellar medium there are molecules with an unpaired electron spin. 
From Eq. ()3.2|) it is clear that a detectable Zeeman splitting (i.e. comparable 
or possibly larger than the Doppler broadening) can be generically obtained 
for magnetic fields where ^ 10^^ G, i.e. magnetic fields of the order of 
the m G. Molecules with an unpaired electron spin include OH, CN, CH and 
some other. In the past, for instance, magnetic fields have been estimated in 
OH clouds (see jT^j and references therein). Magnetic fields of the order of 
10 m G have been detected in interstellar H2O maser clumps (with typical 
densities O{10^^cm~^)) More recently attempts of measuring magnetic 
fields in CN have been reported [21]. The possible caveat with this type 
of estimates is that the measurements can only be very local: the above 
mentioned molecules are much less common than neutral hydrogen and are 
localized in specific regions of the interstellar medium. In spite of this caveat, 
Zeeman splitting measurements can provide reliable informations on the local 
direction of the magnetic field. This determination is important in order to 
understand the possible origin of the magnetic field. This aspect will be 
discussed, in more detail, when describing the magnetic field of the Milky 
Way. 

3.2 Synchrotron emission 

The first experimental evidence of the existence of large-scale magnetic fields 
in external galaxies came, historically, from the synchrotron emission of Crab 
nebula The emissivity formula (i.e. the energy emitted from a unit 

volume, per unit time, per unit frequency interval and per unit sohd angle) 
for the synchrotron is a function of B± and and of the relativistic electron 
density, namely 



where u is the frequency. Wo is a proportionality constant depending only 
upon the spectral index a which also determines the (isotropic) relativistic 
electron number density distribution for electrons with energies in the range 



(3.4) 



dE: 



nJE = riQE-'^dE. 



(3.5) 
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A useful estimate for the maximal emission from electrons of energy E is 
given by |2I] 



In generic terms the synchrotron emission is then sensitive to the total trans- 
verse magnetic field. In oversimplified terms the measurement proceeds in 
three steps: 

• estimate of the emission intensity; 

• from the estimate of the emission intensity, defining as L the typical 
size of the source, the quantity 5j"+_L^/^WoL can be deduced; 

• the magnetic field can then be obtained by specifying the electron den- 
sity distribution of the source. 

While the second step is rather well defined [a can indeed be determined 
from the observed emission spectrum ^) the third step has to be achieved 
in a model-dependent way. The relativistic electron density is sometimes 
estimated using equipartition, i.e. the idea that magnetic and kinetic energy 
densities may be, after all, comparable. Equipartition is not always an em- 
pirical evidence, but can certainly be used as a working hypothesis which 
may or may not be realized in the system under consideration. For instance 
equipartition probably holds for the Milky Way but it does not seem to be 
valid in the Magellanic Clouds |26j. The average equipartition field strengths 
in galaxies ranges from the 4/i G of M33 up to the 19 jj G of NGC2276 ^ EH] . 

In order to illustrate the synchrotron emission take for instance the MW. 
Away from the region of the galactic plane synchrotron emission is the domi- 
nant signal of the galaxy for frequencies below 5 GHz. Between 1.4 and 5 GHz 
the spectral dependence is C(z/~°-^). It should be noticed that the spectral 
index of synchrotron emission may change in different frequency ranges and 
also in different regions of the sky due to the variation of the magnetic field 
pUj . This aspect is crucial in the experimental analysis of CMB anisotropy 
experiments in the GHz region. In fact, in the context of CMB studies the 
synchrotron is a foreground which should be appropriately subtracted. The 
way this is done is by extrapolation of lower frequency measurements. The 

^Typical galactic values of a are between 2.4 and 2.8. An important caveat to this 
statement, relevant for foreground extraction (for instance in studies of CMB anisotropics) 
is that the synchrotron spectrum may have a break. 




(3.6) 
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problem of proper subtraction of synchrotron foreground is particularly acute 
in the case of the PLANCK mission 

The synchrotron has an intrinsic polarization which can give the orienta- 
tion of the magnetic field, but not the specific sign of the orientation vector. 
The strength of synchrotron polarization is proportional to |-B±/-Btot,±P, i-e. 
the ratio between the magnetic energy densities of the uniform and total 
magnetic field. Thus, according to Eq. (j3.1|) . the synchrotron polarization is 
sensitive to the random component of the large-scale magnetic field. 



3.3 Faraday Rotation Measurements 

To infer the magnitude of the magnetic field strength Faraday effect has been 
widely used. When a polarized radio wave passes through a region of space 
of size 6i containing a plasma with a magnetic field the polarization plane of 
the wave gets rotated by an amount 

50occuBMj 6i, (3.7) 

which is directly proportional to the square of the plasma frequency ^ Up (and 
hence to the electron density) and to the Larmor frequency ub (and hence 
to the magnetic field intensity). A linear regression connecting the shift in 
the polarization plane and the square of the wavelength A, can be obtained: 

= 00 + RMAl (3.8) 

By measuring the relation expressed by Eq. ()3.8|1 for two (or more) sepa- 
rate (but close) wavelengths, the angular coefficient of the regression can be 
obtained and it turns out to be 

in units of rad/m^ when all the quantities of the integrand are measured in 
the above units. The explicit dependence of the red-shift can be also easily 
included in Eq. fl3.9|l . Notice, in general terms that the RM is an integral 

■^The satellite mission PLANCK EXPLORER will provide, after 2008, full sky maps in 
nine frequency channels ranging from 30 to 900 GHz. 
®See Eq. H4.5|l in the following Section. 
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over distances. Thus the effect of large distances will reflect in high values of 
the RM. Furthermore, the Faraday effect occurs typically in the radio (i.e. 
cm < A < m), however, some possible applications of Faraday effect in the 
microwave can be also expected (see Section 8 of the present review). 

The shift in the polarization plane should be determined with an accuracy 
greater than 5(f) ~ ±7r. Otherwise ambiguities may arise in the determination 
of the angular coefficient appearing in the linear regression of Eq. (j3.8|) . This 
aspect is illustrated in Fig. [T] which is rather standard but it is reproduced 
here in order to stress the possible problems arising in the physical determi- 
nation of the RM if the determination of the shift in the polarization plane 
is not accurate. 

The RM defined in Eq. (|3.9p not only the magnetic field (which should be 
observationally estimated), but also the column density of electrons. From 
the radio-astronomical observations, different techniques can be used in order 
to determine Ue- One possibility is to notice that in the observed Universe 
there are pulsars. Pulsars are astrophysical objects emitting regular pulses 
of electromagnetic radiation with periods ranging from few milliseconds to 
few seconds. By comparing the arrival times of different radio pulses at 
different radio wavelengths, it is found that signals are slightly delayed as 
they pass through the interstellar medium exactly because electromagnetic 
waves travel faster in the vacuum than in an ionized medium. Hence, from 
pulsars the column density of electrons can be obtained in the form of the 
dispersion measure, i.e. DM oc / n^M. Dividing the RM by DM, an estimate 
of the magnetic field can be obtained. Due to their abundance, pulsars lead 
to the best determination of the magnetic field in the galactic disk ^]. In 
Fig. 121 (adapted from [3T]) a map of the antisymmetric RM sky is reported. 
In the picture the open circles denote negative RM while filled circles denote 
positive RM. The size of the circle is proportional to the magnitude of the 
RM. The convention is, in fact, to attribute negative RM to a magnetic field 
directed away from the observer and positive RM if the magnetic field is 
directed toward the observer. 

As in the case of synchrotron emission also Faraday rotation measure- 
ments can be used as a diagnostic for foreground contamination. The idea 
would be, in this context to look for cross-correlations in the Faraday rota- 
tion measure of extra-galactic sources and the measured microwave signal at 
the same angular position. A recent analysis has been recently reported 

If magnetic field or the column density change considerably over the in- 
tegration path of Eq. (j3.9|) one should probably define and use the two-point 
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Figure 1: The possible ambiguities arising in practical determinations of 
the RM are illustrated. The RM is the angular coefficient of the linear 
regression expressed by Eq. ()3.8p . Clearly it is not necessary to know the 
initial polarization of the source to determine the slope of a straight line in 
the (0, A^) plane, but it is enough to measure at two separate wavelength. 
However, if the accuracy in the determination of is of the order of vr the 
inferred determination of the angular coefficient of the linear regression (|3.8|) 
is ambiguous. 
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SP: b=-90° 



Figure 2: The filtered RM distribution of extragalactic radio sources. Tlie 
antisymmetric distribution is clear especially from the inner galactic quad- 
rant. This picture is adapted from 

function of the RM, i.e. 

C(r) = (RM(f + f)RM(f)). (3.10) 

The suggestion to study the mean-squared fluctuation of the RM was pro- 
posed [HSl IHlj • More recently, using this statistical approach particularly ap- 
propriate in the case of magnetic fields in clusters (where both the magnetic 
field intensity and the electron density change over the integration path), 
Newmann, Newmann and Rephaeli quantified the possible (statistical) 
uncertainty in the determination of cluster magnetic fields (this point will 
also be discussed later). The rather ambitious program of measuring the 
RM power spectrum is also pursued j2S|E3|- the analysis of correla- 

tions in the RM has been discussed. 

3.4 Magnetic fields in the Milky Way 

Since the early seventies |12^ IHH] the magnetic field of the Milky way (MW) 
was shown to be parallel to the galactic plane. RM derived from pulsars allow 
consistent determinations of the magnetic field direction and intensity [101 
1^ . In the Milky Way, the uniform component of the magnetic field is in the 
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plane of the galactic disk and it is thought to be directed approximately along 
the spiral arm. There is, though, a slight difference between the northern 
and southern hemisphere. While in the southern hemisphere the magnetic 
field is roughly 2/iG, the magnetic field in the northern hemisphere is three 
times smaller than the one of the southern hemisphere. The magnetic field 
of the MW is predominantly toroidal ^. However, looking at the center of 
the galaxy, a poloidal component (typically 1/3 of the toroidal one) has 
been detected. In fact, in the central 100 pc of the MW there are in fact 
pooloidal (dipole) fields whose origin is probably primordial see also 
[43, 44J)- The reason for this statement is that through the usual plasma 
physics mechanisms (like the dynamo theory to be discussed later in Section 
4) it is hard to account for a sizable poloidal component, even if localized in 
the central region of the MW. 

The common practice is to classify large-scale magnetic fields in spiral 
galaxies according to the parity of the toroidal field under reflections of the 
azimuthal angle [^El- According to some authors this distinction palys a 
crucial role in order to assess the primordial (or non-primordial ) nature of 
the observed field (see, for instance, |EIIIH|)- In short, the virtues and limi- 
tations of the previous distinction can be summarized as follows. Suppose to 
plot the toroidal component of the measured magnetic field (for instance of 
the MW) as a function of the azimuthal angle 9. This is, in fact, equivalent, 
to plot the RM{6) for instance near the equatorial plane of the galaxy. Then 
two qualitatively different situations can arise. In the first case RM{6) upon 
rotation of tt around the galactic center is even and this corresponds to ax- 
isymmetric spiral galaxies (ASS). In the second case RM{9) upon rotation 
of TT around the galactic center is odd and this corresponds to bisymmetric 
spiral galaxies (BSS). This sharp distinction is invalidated by two experimen- 
tal facts. First of all there may be extra phases so that it is not clear, in 
practice, if RM{6) is, predominantly, even or odd. In spite of possible am- 
biguities (which are, however, an experimental problem) the distinction may 
be very useful: if the magnetic fields would originate through some dynamo 
amplification the preferred configuration will be of the ASS-type and not of 
the BSS-type. There are strong indications that the MW is a BSS spiral 
PUI I31j . A related issue is the reversal of the magnetic field from one spiral 
arm to the other. In [M] 63 rotation measures from polarization observations 

^Recall that the toroidal field is defined as the vector sum of the radial and azimuthal 
components, i.e. i?toroidai = Beee + Brer- 
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of southern pulsars have been reported. For the galactic disk, convincing ev- 
idence of field reversal near the Perseus arm is presented. In the solar circle 
three reversals are observed: near the Carina- Sagittarius arm, near the Crux- 
Scutum arm, and possibly a third one in the Norma arm. These reversals 
are claimed to be consistent with BSS models. 

Since this is a relevant point, the possible controversies arising in the 
analysis of the magnetic field of the MW will now be swiftly mentioned. 
It is appropriate, for this purpose, to give a geometrical characterization of 
the arm structure of the MW [45j- The spiral structure of the galaxy can 
be described in terms of a two-dimensional coordinate system {x, y) whose 
origin (0, 0) is at the galactic center. In this coordinate system the sun is 
located at (0, 8) where the coordinates are expressed in kpc. For each spiral 
arm the equation reduces to four curves, each rotated by 7i/2 of the form: 

x = rcos6', y = rsm6, r = rQe'^^^^^°\ (3-11) 

In Eq. (Km ro ~ 2.3kpc, Oq ~ 0, -7r/2, -tt, -3/27r, k ~ 0.21. In Fig. Hthe 
map of the spiral arms is illustrated according to the model of Eq. (j3.1ip . 
Fig. El is adapted from the recent paper of Vallee [121 • The region with 
dashed lines is poorly known from the observational point of view ("zona 
galactica incognita" in the terminology of Vallee 03]). In terms of Fig. El the 
basic question related to the magnetic field structure concerns the relative 
orientation of the magnetic field direction between one arm and the nearest 
one. As mentioned abovej^ni IHI] some studies suggest that three reversals 
are present. Some other studies [U] based on a statistical re- analysis of 
the most recent data from 1995 to 2001 suggest that only one reversal is 
observed. The two large arrows in Fig. El represent the field directions in 
the Carina-Sagittarius and in the Perseus arms. In jlH] it has been argued 
that the magnetic field of the MW has the structure of axisymmetric spiral. 
However, the presence of one field reversal seems to be not fully consistent 
with this interpretation. 

Recently jlHl Han, Manchester, Lyne and Qiao, using the Parkes multi- 
beam pulsar survey, provided further observational support for the detection 
of a counterclockwise magnetic field in the Norma spiral arm. The morpho- 
logical properties of the magnetic field in the Carina-Sagittarius and in the 
Crux-Scutum arms have been confirmed jlHl- These results were obtained 
from the analysis of pulsar rotation measures. 

If the magnetic field of the MW flips its direction from one spiral arm 
to the other, then, as pointed out by Sofue and Fujimoto jlT] the galactic 
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Figure 3: The map of the MW is illustrated. Following [IHI the origin of the 
two-dimensional coordinate system are in the Galactic center. The two large 
arrows indicate one of the possible (3 or 5) field reversals observed so far. 
The field reversal indicated in this figure is the less controversial one. 

magnetic field should probably be associated with a BSS model. In the Sofue- 
Fujimoto model the radial and azimuthal components of the magnetic field 
in a bisymmetric logarithmic spiral configuration is given through 



where tq ~ 10.5 kpc is the galactocentric distance of the maximum of the 
field in our spiral arm, (3 = 1/ tanp and p is the pitch angle of the spiral. 
The smooth profile /(r) can be chosen in different ways. A motivated choice 
is llHllli (see also jSl) 



where ri = 8.5 kpc is the distance of the Sun to the galactic center and r2 = 2 
kpc. The original model of SF does not have dependence in the z direction, 




(3.12) 




(3.13) 



18 



however, the z dependence can be included and also more complicated models 
can be built 07j. Typically, along the z axis, magnetic fields are exponentially 
suppressed as exp [— z/zq] with ~ 4 kpc. 

In the BSS model there are then three issues to be determined in order 
to specify the parameters (SO]: 

• the number and location of field reversals; 

• the value of the pitch angle p introduced in Eq. ()3.12|) : 

• the field strength in the disk. 

The answers to these three questions are, in short, the following: 

• 3 to 5 field reversals have been claimed; 

• the pitch angle is determined to be p = —8.2'^ ± 0.5° jHH ED]; 

• the strength of the regular field is, as anticipated 1.8 ± 0.3 /iG. 

It is interesting to notice that the total field strength, in the notation of 
Eq. p.ip can reach even 6 /i G indicating, possibly, a strong stochastic 
component. Differently from other spirals, the Milky Way has also a large 
radio halo. The radio-halo indicates the large-scale height of the magnetic 
field and its origin is unclear As far as the stochastic component of the 
galactic magnetic field is concerned, the situation seems to be, according 
to the reported results, still unclear [HIlEllESl- It is, at present, not fully 
understood if the stochastic component of the galactic magnetic field is much 
smaller than (or of the same order of) the related homogeneous part as 
implied by the estimates of the total field strength. 

The structure of magnetic fields can be relevant when investigating the 
propagation of high-energy protons [SH ESI as noticed already long ago [HE] • 
This aspect leads to crucial (and structural) ambiguities in the analysis of 
the propagation of charged particles of ultra-high energy. 

The current observations of ultrahigh energy cosmic rays (UHECR) do 
not lead to a firm evidence of the existence of a cutoff of the cosmic ray 
spectrum between 10^^ and 10^*^ eV. This is the celebrated Greisen, Zat- 
sepin, Kuzmin (GZK) cutoff jSZlEH]- If the source distribution of UHECR is 
isotropic and homogeneous because photopro duct ion interactions on the mi- 
crowave background, then the GZK cutoff should be present. The isotropiza- 
tion of UHECR can be explained by scattering in large magnetic fields. This 
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suggestion can be achieved in the presence of a sizable galactic halo. Mag- 
netic fields may change not only the local intensity of UHECR but also their 
energy spectrum JE^. Various papers analyzed the acceleration and prop- 
agation of UHECR in magnetized structures (see, for instance, jHUl IHH 112] 
and references therein). 

3.5 Magnetic fields in external galaxies 

An excellent presentation of the evidence of cosmic magnetism in nearby 
galaxies has been given by Beck et al.p3] and, in a more concise form, by 
Widrow [63J. The past few years have witnessed a long discussion on the 
primordial nature of magnetic fields in nearby galaxies whose features are 
partially reported in P3]. In analogy with the problem discussed in the case 
of the MW two different opinions have been confronted. Theorists believing 
the primordial nature of the galactic magnetic field supported the BSS model. 
On the contrary, dynamo theorists supported the conclusion that most of the 
nearby galaxies are ASS. As previously described, the variation of the RM 
as a function of the azimuthal angles has been used in order to distinguish 
ASS from BSS [64^i65j. The simplistic summary of the recent results is that 
there are galaxies where BSS structure dominates (like M81,M51) jM]- At 
the same time in some galaxies the ASS dominates [U7j . 

3.6 Magnetic fields in Abell clusters 

Already the short summary of the main experimental techniques used for the 
detection of large scales magnetic fields shows that there may be problems in 
the determination of magnetic fields right outside galaxies. There magnetic 
fields are assumed to be often of nG strength. However, due to the lack of 
sources for the determination of the column density of electrons, it is hard 
to turn the assumption into an experimental evidence. Magnetic fields in 
clusters have been recently reviewed by Giovannini [UHl and, more extensively, 
by Carilli and Taylor |69j . 

Since various theoretical speculations suggest that also clusters are mag- 
netized, it would be interesting to know if regular Abell clusters posses large 
scale magnetic fields. Different results in this direction have been reported 
[701 [711 [72] (see also [ZH]). Some studies during the past decade [7ni[7T] dealt 
mainly with the case of a single cluster (more specifically the Coma cluster). 
The idea was to target (with Faraday rotation measurements) radio sources 
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inside the cluster. However, it was also soon realized that the study of many 
radio sources inside different clusters may lead to experimental problems due 
to the sensitivity limitations of radio-astronomical facilities. The strategy is 
currently to study a sample of clusters each with one or two bright radio- 
sources inside. 

In the past it was shown that regular clusters have cores with a detectable 
component of RM KT^J [72| • Recent results suggest that Gauss magnetic 
fields are indeed detected inside regular clusters ITS] . Inside the cluster means 
in the intra-cluster medium. Therefore, these magnetic fields cannot be as- 
sociated with individual galaxies. 

Regular Abell clusters with strong x-ray emission were studied using a 
twofold technique [731171]. From the ROSAT full sky survey, the electron den- 
sity has been determined [73]. Faraday RM (for the same set of 16 Abell clus- 
ters) has been estimated through observations at the VLA radio-telescope. 

The amusing result (confirming previous claims based only on one cluster 
[701 I7T] ) is that x-ray bright Abell clusters possess a magnetic field of fi 
Gauss strength. The clusters have been selected in order to show similar 
morphological features. All the 16 clusters monitored with this technique 
are at low red-shift {z < 0.1) and at high galactic latitude (|6| > 20°). 

These recent developments are rather promising and establish a clear 
connection between radio-astronomical techniques and the improvements in 
the knowledge of x-ray sky. There are various satellite missions mapping 
the x-ray sky at low energies (ASCA, CHANDRA, NEWTON 1°). There is 
the hope that a more precise knowledge of the surface brightness of regular 
clusters will help in the experimental determination of large scale magnetic 
fields between galaxies. It should be however mentioned that evidence for 
the presence of relativistic electrons and magnetic fields in clusters was di- 
rectly available even before, from measurements of extended regions of radio 
synchrotron emission (for frequencies 10~^ < u < 1 GHz) [7B] . 

It is interesting to notice that intra-cluster magnetic fields of G strength 
can induce Faraday rotation on CMB polarization. By combining infor- 
mations from Sunyaev-Zeldovich effect and X-ray emission from the same 
clusters, it has been recently suggested that a richer information concerning 
electron column density can be obtained [^. In Fig. |3]the results reported 

^°ASCA is operating between 0.4 AND 10 keV and it is flying since February 1993. 
CHANDRA (NASA mission) and NEWTON (ESA mission) have an energy range com- 
parable with the one of ASCA and were launched, almost simultaneously, in 1999. 
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in [7^1 are summarized. In Fig. |3]the RM of the sample of x-ray bright Abell 
clusters is reported after the subtraction of the RM of the galaxy. At high 
galactic latitude (where all the observed clusters are) the galactic contribu- 
tion is rather small and of the order of 9.5rad/m^. 

The results reported in In Fig. HI the open points represent sources 
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Figure 4: From Ref. the RM deduced from a sample of 16 X-ray bright 
Abell clusters is reported as a function of the source impact parameter. 

viewed through the thermal cluster gas, whereas the full points represent 
control sources at impact parameters larger than the cluster gas. The excess 
in RM attributed to clusters sources is clearly visible. 

Using the described techniques large scale magnetic fields can be observed 
and studied in external galaxies, in clusters and also in our own galaxy. While 
the study of external galaxies and clusters may provide a global picture of 
magnetic fields, the galactic observations performed within the Milky Way 
are more sensitive to the local spatial variations of the magnetic field. For 
this reasons local and global observations are complementary. The flipped 
side of the coin, as we will see in the second part of the present Section, is 
that the global structure of the magnetic field of our galaxy is not known 
directly and to high precision but it is deduced from (or corroborated by) 
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the global knowledge of other spiral galaxies. 
3.7 Magnetic fields in superclusters 

Magnetic fields in the local supercluster are notoriously difficult to measure. 
In the absence of direct measurements it is the opinion of the author that 
these fields cannot be assumed to be small, i.e. -C fiG. In this situation it 
seems useful to consider the possibility that these fields are large. In [23] 
this point of view has been taken in connection with the propagation of 
the most energetic cosmic rays. This approach is desirable since the (few) 
reliable claims of measurements of large-scale magnetic fields in the Local 
Supercluster seem to indicate that these fields are rather strong and more in 
the range of the n G than in the range of the n G. 

Recently [7H] the possible existence, strength and structure of magnetic 
fields in the intergalactic plane, within the Local Supercluster, has been 
scrutinized. The local supercluster is centered approximately at the VIRGO 
cluster, about 18 Mpc from the Local Group. A statistically significant 
Faraday screen acting on the radio-waves coming from the most distance 
sources has been found. This analysis supports the existence of a regular 
magnetic field of 0.3/x G in the Local Supercluster. 

In the past detection of radio emission in the Coma supercluster has been 
reported [THI. The plane of the Coma supercluster is defined by the Coma 
cluster and by the Abell cluster 1367. The observed magnetic field has been 
extimated to be of the order of 0.5/x G. 

In jHUj using numerical simulations the spectrum of protons above 10^^ 
eV has been determined under the assumption that the injection spectrum 
is determined by a discrete set of continuously emitting sources. The sources 
follow the profile of the local Supercluster. It has been shown that if |-B| < 
0.05/iG the source distribution assumed in the Supercluster is inconsistent 
with the observations. On the other hand, if \B\ < 0.3/iG, 10 sources in the 
local Supercluster would lead to spectra consistent with the observations. 
This analysis seems to fit with [7HI and with the general arguments put for- 
ward in [55j. In a complementary perspective, the authors of claim to be 
able to simulate the evolution of large-scale magnetic fields from reasonably 
high red-shifts. The obtained results, according to [HIl, correctly reproduce 
the present large-scale structure of cluster and supercluster magnetic fields 
both in their observed and unobserved features. 
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4 Globally Neutral Plasmas 



The simplified discussion of tlie plasma evolution is often related to magneto- 
hydrodynamics (MHD) which is a one-fluid plasma description holding under 
very specific assumptions. The goal of the first part of the present Section 
(together with the related Appendix A) will be to present the reduction of 
the common two-fluids plasma dynamics to the simplified case of one-fluid 
description provided by MHD. 



4.1 Qualitative aspects of plasma dynamics 

Consider, as a starting point, a globally neutral plasma of electrons and 
ions (in the simplest case protons) all at the same temperature Tq and with 
mean particle density uq. Suppose, for simplicity, local thermodynamical 
equilibrium The charge densities will be given, respectively, by 

rze(0)=noe'^*/^°^(l + ^), 

np(0) = noe-<^/^«^(l-^), (4.1) 

where is the electrostatic potential. In Eqs. ()4.1|) it has been assumed that 
the plasma is weakly coupled, i.e. |e0/To| <^ 1. A test charge qt located in 
the origin will then experience the electrostatic potential following from the 
Poisson equation, 

= 47reK(0) - np(0)] - 47rq,6^^\x), (4.2) 
or, using Eqs. pTjl into Eq. (jOI) 



V'0-40=-47rgt5(=^)(f), (4.3) 



where 



Ad = i/^^, (4.4) 

is the Debye length. For a test particle the Coulomb potential will then have 
the usual Coulomb form but it will be suppressed, at large distances by a 
Yukawa term, i.e. e~'^^'*'°. 



'^This assumption is violated in the realistic situations. 
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In the interstellar medium there are three kinds of regions which are 
conventionally defined, namely H2 regions (where the Hydrogen is predom- 
inantly in molecular form), H° regions (where Hydrogen is in atomic form) 
and H"*" regions (where Hydrogen is ionized). Sometimes H"*" regions are 
denoted with HI and H2 regions with HII. In the H+ regions the typical tem- 
perature To is of the order of 10-20 eV while for no let us take, for instance, 
no ~ 3 X 10-2cm-3. Then Ad ~ 30km. 

For r Xb the Coulomb potential is screened by the global effect of the 
other particles in the plasma. Suppose now that particles exchange momen- 
tum through two-body interactions. Their cross section will be of the order 
of ttgjn/To and the mean free path will be f^fp ~ ^o^/lc^em^o); i-e. recalling 
Eq. (j4.4p Ad ^ ^mfp- This means that the plasma is a weakly collisional 
system which is, in general, not in local thermodynamical equilibrium. This 
observation can be made more explicit by defining another important scale, 
namely the plasma frequency which, in the system under discussion, is given 
by 

/4vrnoe2 

uJp = \ , 4.5 

V "^e 

where rrie is the electron mass. The plasma frequency is the oscillation 
frequency of the electrons when they are displaced from their equilibrium 
configuration in a background of approximately fixed ions. Recalling that 



'To/nie is the thermal velocity of the charge carriers, the collision 
frequency Uc — fther/^mfp is always much smaller than Up ~ fther/^D- Thus, 
in the idealized system described so far, the following hierarchy of scales 
holds: 

Ad < Cfp, t^c < ujp, (4.6) 

which means that before doing one collision the system undergoes many os- 
cillations, or, in other words, that the mean free path is not the shortest scale 
in the problem. Usually one defines also the plasma parameter Af = no ^A^^, 
i.e. the number of particles in the Debye sphere. In the approximation of 
weakly coupled plasma A/" ^ 1 which also imply that the mean kinetic energy 
of the particles is larger than the mean inter-particle potential, i.e. |e0| <^ Tq 
in the language of Eq. 1)4.11) . 
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4.2 Kinetic and fluid descriptions 

From the point of view of the evolution of large-scale magnetic fields, the 
galaxy is globally neutral system of charged particles with typical rotation 
period of 3 x 10^ yrs. Two complementary descriptions of the plasma can 
then be adopted. The first possibility is to study the full system of kinetic 
equations (the Vlasov-Landau description [S21IHS1E11 for the one particle dis- 
tribution function). In the kinetic description the "observables" of the plasma 
dynamics are related to the various moments of the distribution function. In 
a complementary perspective, it is also possible to study directly the evolu- 
tion equations of the various moments of the distribution function (charge 
and matter densities, momentum transfer etc.). The resulting description is 
an effective one since not all the moments of the distribution function can be 
kept. The necessity of a consistent truncation in the hierarchy of the various 
moments is usually called closure problem. MHD is an effective description 
which holding for 

L > Ad, t > uj-\ (4.7) 

where L and t are the typical length and time scales of the problem. The 
MHD approximation holds for typical length scales larger than the Debye 
scale and typical frequencies much smaller than the plasma frequency. In 
particular, the spectrum of plasma excitations obtained from the full kinetic 
description matches with the spectrum of MHD excitations but only at suf- 
ficiently small frequencies (typically of the order of the ionoacoustic waves) 



If the plasma is weakly collisional and out of thermodynamical equilib- 
rium, an approximation where the evolution equations are still tractable is 
given in terms of the one-particle distribution function fi{x,v,t) where the 
subscript denotes the given charge carrier (for instance electrons and ions). 
The one-particle distribution obeys the Vlasov-Landau equation for each par- 
ticle species of charge qi and of mass mi: 

^+v-V,fi + ^{E + vxB)-V,n=(^) . (4.8) 



at rrii \ at J 



coll 



The term appearing at the right-hand side of Eq. ()4.8p is the binary col- 
lision term. In a weakly coupled plasma, the long-range force is given by 
the Coulomb interaction, while the short-range component arises thanks to 
binary collisions. In other words the electric and magnetic fields appearing 
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in Eq. (|4.8|) are mean fields obeying mean Maxwell's equations: 




i 



(4.9) 




Vx5 = — + 47r^ TiiQi / d^vvfi{x, v, t), 



(4.11) 



(4.10) 







(4.12) 



where ni is the mean particle density. In the context of the kinetic approach 
the evolution equation should be solved self-consistently with the Maxwell's 
equations. Clearly, in a system where the density of particles is small the 
mean field given by the distant particles is more important than the the force 
produced by the closest particles. In practice the Vlasov-landau system can 
be linearized around some equilibrium value of the one-particle distribution 
function and specific examples will be discussed later on. 

The one-particle distribution function is not directly observable. The 
directly observable quantities are the the various moments of the one-particle 
distribution, for instance, the particle density. 



and the related matter and charge densities for each particle species (i.e. 
Pm,i{x,t) = mini{x,t) and Pq^i{x,t) = gini(x, t)) or summed over the differ- 
ent species (i.e. = J2iPm,i and Pq = J2iPq,i)- However, the dynamical 
informations on the plasma evolution are not only encoded in the moment 
of order zero but also in the higher moments (the first, the second and even 
higher) . 

For some class of problems the kinetic approach is to be preferred. For 
instance the kinetic approach is mandatory for all the high frequency phe- 
nomena (like, for instance, equilibration of charge and current density fluc- 
tuations). As far as the low- frequency phenomena are concerned one can 
imagine to obtain a reduced (fluid-like) description. 

4.3 From one-fluid equations to MHD 

The evolution equations for the one-fluid variables are obtained from the two- 
species kinetic description by some algebra which is summarized in the first 




(4.13) 
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part of Appendix A. The bottom line of the derivation is that the moments of 
the one-particle distribution function obey a set of partial differential equa- 
tions whose solutions can be studied in different approximations. One of 
these approximations is the MHD description. The set of one-fluid equations 
(derived in Appendix A) can be written as 

^^"^ +V-{pr.v) = 0, (4.14) 



dt 

+ V-/=0. (4.15) 



dt 

'dv 



Pn 



v = JxB-VP, (4.16) 



E + VX B = -J+ —{fx B -VPe). (4.17) 



cr eriq 



supplemented by the Maxwell's equations 

V-^ = 47rpq, V-5 = (4.18) 

f)B 

VxE + — = 0, (4.19) 

^ dE 

V X 5 = — + 47r J (4.20) 

at 

Eqs. pmj) and (prT5|) are the mass and charge conservation, Eq. ()4.16|1 is 
the equation for the momentum transfer and Eq. ()4.17j) is the generalized 
Ohm law. In Eqs. ()4.14j) - ()4.17j) v is the one-fluid velocity field, P is the total 
pressure and Pg is the pressure of the electrons. 

If the plasma is globally neutral the charge density can be consistently 
neglected. However, global neutrality is plausible only for typical lengths 
much larger than Ad where the collective properties of the plasma are the 
leading dynamical effect. Under this assumption, Eqs. (|4.15|) and (j4.18|) 
imply 

V-J = 0, V-^ = 0, (4.21) 

which means that MHD currents are, in the first approximation, solenoidal. 

The second step to get to a generalized form of the MHD equations is to 
neglect the terms responsible for the high frequency plasma excitations, like 
the displacement current in Eqs. ()4.20|) whose form becomes then 

V X B = 4nJ. (4.22) 
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Given the two sets of consistent approximations discussed so far, the form 
taken by the generahzed Ohm law may not be unique and, therefore, differ- 
ent MHD descriptions arise depending upon which term at the right hand 
side of Eq. ()4.17|) is consistently neglected. The simplest (but often not real- 
istic) approximation is to neglect all the terms at the right hand side of Eq. 
(j4.17p . This is sometimes called ideal MHD approximation. In this case the 
Ohmic electric field is given by ^ —v x B. This ideal description is also 
called sometimes superconducting MHD approximation since, in this case, 
the resistive terms in Eq. (j4.17|) . i.e. J /a go to zero (and the conductivity 
a — i> oo). In the early Universe it is practical to adopt the superconducting 
approximation (owing to the large value of the conductivity) but it is also 
very dangerous since all the 0{l/a) effects are dropped. 

A more controllable approximation is the real (or resistive) MHD. The 
resistive approximation is more controllable since the small expansion pa- 
rameter (i.e. l/o") is not zero (like in the ideal case) but it is small and finite. 
This allows to compute the various quantities to a a given order in l/a even 
if, for practical purposes, only the first correction is kept. In this framework 
the Ohmic current is not neglected and the generalized Ohm law takes the 
form 

J^a{E + vxB). (4.23) 

Comparing Eq. ()4.23|) with Eq. ()4.17|) . it is possible to argue that the resis- 
tive MHD scheme is inadequate whenever the Hall and thermoelectric terms 
are cannot be neglected. Defining ujb = eB/m as the Larmor frequency, 
the Hall term J x B can be neglected provided Luj'^/ujbVq ^ 1 where L us 
the typical size of the system and vq the typical bulk velocity of the plasma. 
With analogous dimensional arguments it can be argued that the pressure 
gradient, i.e. the thermoelectric term, can be consistently neglected pro- 
vided LvQUJEfne/Te ^ 1. Specific examples where the Hall term may become 
relevant will be discussed in the context of the Biermann battery mechanism. 

The set of MHD equations, as it has been presented so far, is not closed 
if a further relation among the different variables is not specified pZj . Typi- 
cally this is specified through an equation relating pressure and matter den- 
sity (adiabatic or isothermal closures). Also the incompressible closure (i.e. 
V ■ = ) is often justified in the context of the evolution of magnetic fields 
prior to recombination. It should be stressed that sometimes the adiabatic 
approximation may lead to paradoxes. It would correspond, in the context of 
ideal MHD, to infinite electrical conductivity and vanishingly small thermal 
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conductivity (i.e. the electrons should be extremely fast not to feel the re- 
sistivity and, at the same time, extremely slow not to exchange heat among 
them). 

4.4 The Alfven theorems 

The ideal and resistive MHD descriptions are rather useful in order to il- 
lustrate some global properties of the evolution of the plasma which are 
relevant both in the evolution of large-scale magnetic fields prior to recom- 
bination and in the discussion of the gravitational collapse in the presence 
of magnetic fields. These global properties of the plasma evolution go under 
the name of Alfven theorems [SZl (see also [HH] ) ■ 

In the ideal limit both the magnetic flux and the magnetic helicity are 
conserved. This means, 

/ B-dt = ^— [ VxVx5-rfE, (4.24) 

where S is an arbitrary closed surface which moves with the plasma If 
cr — cx) the expression appearing at the right hand side is sub-leading and 
the magnetic flux lines evolve glued to the plasma element. 

The other quantity which is conserved in the superconducting limit is the 
magnetic helicity 

Hm= I d^xA ■ B, (4.25) 
Jv 

where A is the vector potential whose presence may lead to think that the 
whole expression ()4.25|1 is not gauge-invariant. In fact A ■ B is not gauge 
invariant but, none the less, TYa/ is gauge- invariant since the integration 
volume is defined in such a way that the magnetic field B is parallel to the 
surface which bounds V and which we will call dV. If n is the unit vector 
normal to dV then B ■ n = on dV and the gauge dependent contribution 
to the integral appearing in Eq. ()4.25|) vanishes. In physical terms the 
integration may always be performed imagining that the space is sliced in 
flux tubes of the magnetic field. This procedure is correct provided the 
magnetic flux is, at least approximately, conserved as implied by Eq. ()4.24|) . 
The magnetic gyrotropy 

B -V X B (4.26) 

-^^Notice that in Ref. [HE] Eq. I|4.24(l has been derived without the An term at the right 
hand side because of different conventions. 
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it is a gauge invariant measure of the diffusion rate of T^m- In fact, in the 
resistive approximation jHHl 



As in the case of Eq. ()4.25|1 . for cr — oo the magnetic hehcity is approxi- 
mately constant. The value of the magnetic gyrotropy allows to distinguish 
different mechanisms for the magnetic field generation. Some mechanisms are 
able to produce magnetic fields whose flux lines have a non-trivial gyrotropy. 
The properties of turbulent magnetized plasmas may change depending upon 
the value of the gyrotropy and of the helicity [57j . 

The physical interpretation of the flux and magnetic helicity conservation 
is the following. In MHD the magnetic field has to be always solenoidal 
(i.e. V ■ -B = 0). Thus, the magnetic flux conservation implies that, in the 
superconducting limit (i.e. a oo) the magnetic flux lines, closed because 
of the transverse nature of the field, evolve always glued together with the 
plasma element. In this approximation, as far as the magnetic field evolution 
is concerned, the plasma is a collection of (closed) flux tubes. The theorem 
of flux conservation states then that the energetical properties of large-scale 
magnetic fields are conserved throughout the plasma evolution. 

While the flux conservation concerns the energetical properties of the 
magnetic flux lines, the magnetic helicity concerns chiefly the topological 
properties of the magnetic flux lines. If the magnetic field is completely 
stochastic, the magnetic flux lines will be closed loops evolving independently 
in the plasma and the helicity will vanish. There could be, however, more 
complicated topological situations jHHl where a single magnetic loop is twisted 
(like some kind of Mobius stripe) or the case where the magnetic loops are 
connected like the rings of a chain. In both cases the magnetic helicity will 
not be zero since it measures, essentially, the number of links and twists 
in the magnetic flux lines. Furthermore, in the superconducting limit, the 
helicity will not change throughout the time evolution. The conservation 
of the magnetic flux and of the magnetic helicity is a consequence of the 
fact that, in ideal MHD, the Ohmic electric field is always orthogonal both 
to the bulk velocity field and to the magnetic field. In the resistive MHD 
approximation this is no longer true. 

If the conductivity is very large (but finite), the resistive MHD approxi- 
mation suggests, on one hand, that the magnetic flux is only approximately 
conserved. On the other hand the approximate conservation of the magnetic 




(4.27) 
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helicity implies that the closed magnetic loops may modify their topological 
structure. The breaking of the flux lines, occurring at finite conductivity, 
is related to the possibility that bits of magnetic fluxes are ejected from a 
galaxy into the inter-galactic medium. This phenomenon is called magnetic 
reconnection and it is the basic mechanism explaining, at least quahtatively, 
why, during the solar flares, not only charged particles are emitted, but also 
magnetic fields. In the context of large scale magnetic fields the approximate 
(or exact) magnetic flux conservation has relevant consequences for the role 
of magnetic fields during gravitational collapse. 

4.5 Magnetic diffusivity equation 

From Eqs. ()4.22|) and ()4.23|) the Ohmic electric field can be expressed as 

E=^VxB-vxB, (4.28) 
which inserted into Eq. ()4.19|1 leads to the magnetic diffusivity equation 

dB ^ -* 1 n - 

^ = ^x(vxB) + -—V^B. (4.29 
ot Ana 

The first term of Eq. ()4.29p is the dynamo term. The second term of Eq. 
(I4.29|l is the magnetic diffusivity term whose effect is to dissipate the mag- 
netic field. By comparing the left and the right hand side of Eq. ()4.29j) . the 
typical time scale of resistive phenomena is 

^ 47raL^ (4.30) 

where L is the typical length scale of the magnetic field. In a non-relativistic 
plasma the conductivity a goes typically as T^/^ In the case of plan- 

ets, like the earth, one can wonder why a sizable magnetic field can still be 
present. One of the theories is that the dynamo term regenerates contin- 
uously the magnetic field which is dissipated by the diffusivity term pUj . 
In the case of the galactic disk the value of the conductivity is given by 
(T ~ 7 X 10"''Hz. Thus, for L ~ kpc — 10^(L/kpc)^sec. 

^■^It is common use in the astrophysical applications to work directly with rj = (4:Tt<t)~^. 
In the case of the galactic disks = lO^^cm^ Hz. The variable rj denotes, in the present 
review, the conformal time coordinate. 
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In Eq. (l4.3Up the typical time of resistive phenomena has been introduced. 
Eq. (|4.3(J|) can also give the typical resistive length scale once the time-scale of 
the system is specified. Suppose that the time-scale of the system is given by 
tu ~ Hq^ ~ lO^^sec where Hq is the present value off the Hubble parameter. 
Then 

K = (4.31) 

leading to L„ ~ AU. The scale ()4.31|) gives then the upper limit on the 
diffusion scale for a magnetic field whose lifetime is comparable with the age 
of the Universe at the present epoch. Magnetic fields with typical correlation 
scale larger than L^j are not affected by resistivity. On the other hand, 
magnetic fields with typical correlation scale L < L^r are diffused. The value 
L„ ~ AU is consistent with the phenomenological evidence that there are no 
magnetic fields coherent over scales smaller than 10~^ pc. 

The dynamo term may be responsible for the origin of the magnetic field 
of the galaxy. The galaxy has a typical rotation period of 3 x 10® yrs and 
comparing this number with the typical age of the galaxy, (9(10^'^yrs), it can 
be appreciated that the galaxy performed about 30 rotations since the time 
of the protogalactic collapse. 

From Eq. (|4.29|) the usual structure of the dynamo term may be derived 
by carefully averaging over the velocity filed according to the procedure of [23 
inO] . By assuming that the motion of the fluid is random and with zero mean 
velocity the average is taken over the ensemble of the possible velocity fields. 
In more physical terms this averaging procedure of Eq. (j4.29j) is equivalent 
to average over scales and times exceeding the characteristic correlation scale 
and time tq of the velocity field. This procedure assumes that the correlation 
scale of the magnetic field is much bigger than the correlation scale of the 
velocity field which is required to be divergence- less (V ■ v = 0). In this 
approximation the magnetic diffusivity equation can be written as: 

^ = a(VxB) + _V=5. ,4.32) 

where 

a = -'^{v-V X v), (4.33) 

is the so-called dynamo term which vanishes in the absence of vorticity. In 
Eqs. fl4.32|) - ()4.33|) B is the magnetic field averaged over times longer that tq 
which is the typical correlation time of the velocity field. 
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It can be argued that the essential requirement for the consistence of the 
mentioned averaging procedure is that the turbulent velocity field has to be 
"globally" non-mirror symmetric JHl- If the system would be, globally, in- 
variant under parity transformations, then, the a term would simply vanish. 
This observation is related to the turbulent features of cosmic systems. In 
cosmic turbulence the systems are usually rotating and, moreover, they pos- 
sess a gradient in the matter density (think, for instance, to the case of the 
galaxy). It is then plausible that parity is broken at the level of the galaxy 
since terms like Vpm ■ V x t? are not vanishing [IP^. 

The dynamo term, as it appears in Eq. (|4.32j) . has a simple electrody- 
namical meaning, namely, it can be interpreted as a mean ohmic current 
directed along the magnetic field : 

J = -aB (4.34) 

This equation tells us that an ensemble of screw-like vortices with zero mean 
helicity is able to generate loops in the magnetic flux tubes in a plane or- 
thogonal to the one of the original field. Consider, as a simple application of 
Eq. ()4.32|) . the case where the magnetic field profile is given by 

B^{z, t) = f{t) sin kz, By = f{t) cos kz, B,{k, t) = 0. (4.35) 

For this profile the magnetic gyrotropy is non-vanishing, i.e. B ■ V x B = 
kp{t). From Eq. (|4.32|) . using Eq. (j4.35j) f{t) obeys the following equation 

f={ka-^)f (4.36) 



dt \ Aixa, 

admits exponentially growing solutions for sufficiently large scales, i.e. k < 
4:7r\a\a. Notice that in this naive example the a term is assumed to be con- 
stant. However, as the amplification proceeds, a may develop a dependence 
upon i.e. a ^ ao{l-^\B\^)ao[l-^f{t)]. In the case of Eq. ()4.36|) this 
modification will introduce non-linear terms whose effect will be to stop the 
growth of the magnetic field. This regime is often called saturation of the 
dynamo and the non-linear equations appearing in this context are some- 
times called Landau equations jT^] in analogy with the Landau equations 
appearing in hydrodynamical turbulence. 

In spite of the fact that in the previous example the velocity field has 
been averaged, its evolution obeys the Navier-Stokes equation 

+ {v- V)v- v\J^i 



P 



dt 



-Vp+JxB, (4.37) 
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where v is the thermal viscosity coefficient. Since in MHD the matter current 
is solenoidal (i.e. V-(p'u) = 0) the incompressible closure Vp = 0, corresponds 
to a solenoidal velocity field V-tT= 0. Recalling Eq. (j4.22|) . the Lorentz force 
term can be re-expressed through vector identities and Eq. (j4.37|) becomes 



P 



— + (t; ■ V)v - uV^v 









= -V 


p 4- 





{B-V)B. (4.38) 



In typical applications to the evolution of magnetic fields prior to recombi- 
nation the magnetic pressure term is always smaller than the fiuid pressure 
, i.e. p ^ l-^P- Furthermore, there are cases where the Lorentz force term 
can be ignored. This is the so-called force free approximation. Defining the 
kinetic helicity as a; = V x -u, the magnetic diffusivity and Navier-Stokes 
equations can be written in a rather simple and symmetric form 



dB 1 

= V X (t; X 5) + V^B 

at Ana 



^ = V X ({/ X cJ) + z/V^cJ. (4.39) 
at 

In MHD various dimensionless ratios can be defined. The most frequently 
used are the magnetic Reynolds number, the kinetic Reynolds number and 
the Prandtl number: 

R^ = vLsa, (4.40) 

R=^, (4.41) 

Pr = ^ = z,af^y (4.42) 



R \L 



where Lb and are the typical scales of variation of the magnetic and veloc- 
ity fields. In the absence of pressure and density perturbations the combined 
system of Eqs. fl4.22|) and ()4.38|) can be linearized easily. Using then the 
incompressible closure the propagating modes are the Alfven waves whose 
typical dispersion relation is uj"^ = clk^ where Ca = |-B|/\/47rp. Often the 

-"^^Recall that in fusion studies the quantity (3 — 8Tr\B\'^ /p is usually defined. If the 
plasma is confined, then (3 is of order 1. On the contrary, if /3 ^ 1, as implied by the 
critical density bound in the early Universe, then the plasma may be compressed at higher 
temperatures and densities. 
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Lundqvist number is called, in plasma literature [221 IHZ| magnetic Reynolds 
number. This confusion arises from the fact that the Lunqvist number, i.e. 
CaLa, is the magnetic Reynolds number when v coincides with the Alfven 
velocity. To have a very large Lundqvist number implies that the the conduc- 
tivity is very large. In this sense the Lunqvist number characterizes, in fusion 
theory, the rate of growth of resistive instabilities and it is not necessarily 
related to the possible occurrence of turbulent dynamics. On the contrary, 
as large Reynolds numbers are related to the occurrence of hydrodynamical 
turbulence, large magnetic Reynolds numbers are related to the occurence of 
MHD turbulence 

4.6 The dynamo mechanism 

According to the naive description of the dynamo instability presented above 
the origin of large-scale magnetic fields in spiral galaxies can be reduced to 
the following steps: 

• during the 30 rotations performed by the galaxy since the protogalactic 
collapse, the magnetic field should be amplified by about 30 e-folds; 

• if the large scale magnetic field of the galaxy is, today, C(yuG) the 
magnetic field at the onset of galactic rotation might have been even 
30 e-folds smaller, i.e. O{10~^^G); 

• assuming perfect flux freezing during the gravitational collapse of the 
protogalaxy (i.e. a oo) the magnetic field at the onset of gravita- 
tional collapse should be (9(10~^^) G over a typical scale of 1 Mpc. 

This picture is oversimplified and each of the three steps mentioned above 
will be contrasted with the most recent findings. 

The idea that a celestial body may acquire a magnetic field by differential 
rotation can be traced back to the paper of Larmor of 1919 [HI]. One of the 
early ancestors of the dynamo mechanism, can be traced back to the model 
of Fermi and Chandrasekar [7| [S]. In [3 |H| the attempt was to connect 
the existence off the galactic magnetic field with the existence of a galactic 
angular momentum. Later on dynamo theory has been developed in greater 
detail (see [201) its possible application to large-scale magnetic fields has 
been envisaged. 

The standard dynamo theory has been questioned in different ways. Pid- 
dington [^21 ESI pointed out that small-scale magnetic fields can grow large 
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enough (until equipartition is reached) to swamp the dynamo action. The 
quenching of the dynamo action has been numerically shown by Kulsrud and 
Anderson More recently, it has been argued that if the large-scale mag- 
netic field reaches the critical value R^^'^v the dynamo action could also 
be quenched [^3 EE] ■ 

Eq. ()4.29p is exact, in the sense that both v and B contain long and short 
wavelength modes. The aim of the various attempts of the dynamo theory 
is to get an equation describing only the "mean value" of the magnetic field. 
To this end the first step is to separate the exact magnetic and velocity fields 
as 

B={B)+ b, 

v={V) + V, (4.43) 

where (B) and (V^) are the averages over an ensemble of many realizations 
of the velocity field v. In order to derive the standard form of the dynamo 
equations few important assumptions should be made. These assumptions 
can be summarized as follows: 

• the scale of variation of the turbulent motion V should be smaller than 
the typical scale of variation of (B). In the galactic problem {V) is 
the differential rotation of the galaxy, while V is the turbulent motion 
generated by stars and supernovae. Typically the scale of variation of 
V is less than 100 pc while the interesting scales for (B) are larger than 
the kpc; 

• the field 6 is such that \b\ < 

• it should happen that {V ■ V x V) ^ 0. 

• magnetic flux is frozen into the plasma (i.e. magnetic flux is conserved). 

From the magnetic diffusivity equation (j4.29|) . and using the listed assump- 
tions, it is possible to derive the typical structure of the dynamo term by 
carefully averaging over the velocity field according to the procedure outlined 
in [ini 1201 EH- Inserting Eq. into (ji^ and breaking the equation 

^^v is the velocity field at the outer scale of turbulence. 
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into a mean part and a random part, two separate induction equations can 
be obtained for the mean and random parts of the magnetic field 

= V X (^{V) X (B)^ + V X X 6), (4.44) 

^ = V X (V X {B)) + V X {{V) X b) 
at 

+V X (f X 6) - V X X b), (4.45) 



where the (magnetic) diffusivity terms have been neglected. In Eq. ()4.44|1 . 
{V X b) is called mean field (or turbulent) electromotive force and it is the 
average of the cross product of the small-scale velocity field V and of the 
small scale magnetic field b over a scale much smaller than the scale of (B) 
but much larger than the scale of turbulence. Sometimes, the calculation 
of the effect of {V x b) is done in the case of incompressible and isotropic 
turbulence. In this case {V x b) = 0. This estimate is, however, not realistic 
since {B) is not isotropic. More correctly {V x b) should be evaluated 
by using Eq. (j4.45j) which is usually written in a simplified form 

1^ = V X (y X (5)), (4.46) 

where all but the first term of Eq. ()4.45|1 have been neglected. To neglect 
the term V x {(y) x b) does not pose any problem since it corresponds to 
choose a reference frame where iy) is constant. However, the other terms, 
neglected in Eq. ()4.4fi|l . are dropped because it is assumed that \b\ <^ 
This assumption may not be valid all the time and for all the scales. The 
validity of Eq. ()4.46p seems to require that a is very large so that magnetic 
diffusivity can keep always b small [98]. On the other hand [HZl one can argue 
that b is only present over very small scales (smaller than 100 pc) and in this 
case the approximate form of eq. (j4.46|) seems to be more justified. 

From Eqs. (j4.44|) - (j4.46|) it is possible to get to the final result for the 
evolution equation of {B) jHTj as it is usally quoted 



dt 

where 



V X («(5)) + (3V\B) + V X (V) X (5) , (4.47) 



a = -^(y-VxF), (4.48) 
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P = j{V'), (4.49) 

where a is the dynamo term, (3 is the diffusion term and tq is the typical corre- 
lation time of the velocity field. The term a is, in general, space-dependent. 
The standard lore is that the dynamo action stops when the value of the 
magnetic field reaches the equipartition value (i.e. when the magnetic and 
kinetic energy of the plasma are comparable). At this point the dynamo 
"saturates". The mean velocity field can be expressed as {V) ~ i7 x r* where 
|fi(r)| is the angular velocity of differential rotation at the galactocentric ra- 
dius r. In the case of flat rotation curve |^^(r)| = Q{r) ~ which also 
implies that d\n{r)\/dr < 1. 

Eq. ()4.47|) can then be written in terms of the radial and azimuthal com- 
ponents of the mean magnetic field, neglecting, for simplicity, the diffusivity 
term: 

The second equation shows that the a effect amplifies the radial com- 
ponent of the large-scale field. Then, through Eq. ()4.50p the amplification 
of the radial component is converted into the amplification of the azimuthal 
field, this is the Q effect. 

Usually the picture for the formation of galactic magnetic fields is related 
to the possibility of implementing the dynamo mechanism. By comparing the 
rotation period with the age of the galaxy (for a Universe with ~ 0.7, ~ 
0.65 and ~ 0.3) the number of rotations performed by the galaxy since 
its origin is approximately 30. During these 30 rotations the dynamo term 
of Eq. ()4.47|) dominates against the magnetic diffusivity. As a consequence 
an instability develops. This instability can be used in order to drive the 
magnetic field from some small initial condition up to its observed value. 
Eq. (j4.47|) is linear in the mean magnetic field. Hence, initial conditions 
for the mean magnetic field should be postulated at a given time and over a 
given scale. This initial mean field, postulated as initial condition of ()4.47|] 
is usually called seed. 

Most of the work in the context of the dynamo theory focuses on repro- 
ducing the correct features of the magnetic field of our galaxy. The achievable 
amplification produced by the dynamo instability can be at most of 10^^, i.e. 
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e^°. Thus, if the present value of the galactic magnetic field is 10"^ Gauss, 
its value right after the gravitational collapse of the protogalaxy might have 
been as small as 10~^^ Gauss over a typical scale of 30-100 kpc. 

There is a simple way to relate the value of the magnetic fields right after 
gravitational collapse to the value of the magnetic field right before gravita- 
tional collapse. Since the gravitational collapse occurs at high conductivity 
the magnetic flux and the magnetic helicity are both conserved. Right before 
the formation of the galaxy a patch of matter of roughly 1 Mpc collapses by 
gravitational instability. Right before the collapse the mean energy density 
of the patch, stored in matter, is of the order of the critical density of the 
Universe. Right after collapse the mean matter density of the protogalaxy 
is, approximately, six orders of magnitude larger than the critical density. 

Since the physical size of the patch decreases from 1 Mpc to 30 kpc the 
magnetic field increases, because of flux conservation, of a factor (pa/pb)^''^ ~ 
10"^ where pa and pb are, respectively the energy densities right after and right 
before gravitational collapse. The correct initial condition in order to turn 
on the dynamo instability would be \B\ ~ 10^^^ Gauss over a scale of 1 Mpc, 
right before gravitational collapse. 

This last estimate is rather generous and has been presented just in order 
to make contact with several papers (concerned with the origin of large scale 
magnetic fields) using such an estimate. The estimates presented in the 
last paragraph are based on the (rather questionable) assumption that the 
amplification occurs over thirty e-folds while the magnetic flux is completely 
frozen in. In the real situation, the achievable amplification is much smaller. 
Typically a good seed would not be 10~^^ G after collapse (as we assumed 
for the simplicity of the discussion) but rather 

\B\ > IQ-^^G. (4.52) 

The possible applications of dynamo mechanism to clusters is still under 
debate and it seems more problematic [731 [Tlj. The typical scale of the grav- 
itational collapse of a cluster is larger (roughly by one order of magnitude) 
than the scale of gravitational collapse of the protogalaxy. Furthermore, the 
mean mass density within the Abell radius ( ~ 1.5h~^ Mpc) is roughly 10^ 
larger than the critical density. Consequently, clusters rotate much less than 
galaxies. Recall that clusters are formed from peaks in the density field. The 
present overdensity of clusters is of the order of 10^. Thus, in order to get 
the intra-cluster magnetic field, one could think that magnetic flux is exactly 
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conserved and, then, from an intergalactic magnetic field > 10 ^ G an 
intra cluster magnetic field \B\ > 10~^ G can be generated. This simple esti- 
mate shows why it is rather important to improve the accuracy of magnetic 
field measurements in the intra-cluster medium discussed in Section 3: the 
change of a single order of magnitude in the estimated magnetic field may 
imply rather different conclusions for its origin. Recent numerical simulations 
seem to support the view that cluster magnetic fields are entirely primordial 

inni. 

5 Magnetic Fields in Curved Backgrounds 

As the temperature increases above 1 MeV, a relativistic plasma becomes 
almost a perfect conductor. There are two complementary effects associated 
with large-scale magnetic fields in curved backgrounds: 

• large-scale magnetic fields may evolve over a rigid background space- 
time determined, for instance, by the dynamics of barotropic fluids 
and, in this case the energy density of the magnetic field must always 
be smaller than the energy density of the fluid sources; 

• magnetic fields may be so intense to modify the structure of the space- 
time and, in this case, their energy density is comparable with the 
energy density of the other sources of the geometry. 

The first effect will be examined in the present Section, while the discussion 
of the second effect, which is more speculative, will be confined to Section 9. 

5.1 The standard cosmological model 

The standard model of cosmological evolution rests on three important as- 
sumptions |lUUj . The first assumption is that over very large length scales 
(greater than 50 Mpc) the Universe is described by a homogeneous and 
isotropic Friedmann-Robertson- Walker (FRW) line element: 

ds^ = G^^dx^'dx" = a\r])[dr]^ - dx^], (5.1) 

where a{r]) is the scale factor and rj the conformal-time coordinate (notice 
that Eq. ()5.1|) has been written, for simplicity, in the conformally flat case). 
The second hypothesis is that the sources of the evolution of the background 
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geometry are perfect fluid sources. As a consequence the entropy of the 
sources is constant. The third and final hypothesis is that the dynamics of 
the sources and of the geometry is dictated by the general relativistic FRW 



equations 



a 



2 



M^n' = jp, (5.2) 

MUn'-n') = ^ip + p), (5.3) 

p' + 3n{p + p) = 0, (5.4) 

where Ti = (In a)' and the prime denotes the derivation with respect to t]. 
Recall also, for notational convenience, that aH = Ti where H = a/a is the 
conventional Hubble parameter. The various tests of the standard cosmo- 
logical model are well known |1UH llU2j . Probably one of the most stringent 
one comes from the possible distortions, in the Rayleigh- Jeans region, of the 
CMB spectrum. The absence of these distortions clearly rules out steady- 
state cosmological models. In the standard cosmological model one usually 
defines the proper distance of the event horizon at the time ti 

de = a{ti) / — — , (5.5) 

this distance represents the maximal extension of the region over which causal 
connection is possible. Furthermore, the proper distance of the particle hori- 
zon can also be defined: 

d,{t2) = ait,) / \ (5.6) 
Jt^in a[t) 

If the scale factor is parametrized as a{t) ~ t", for < a < 1 the Universe 
experiences a decelerated expansion, i.e. a > and a < while the curvature 
scale decreases, i.e. H < 0. This is the peculiar behaviour when the fluid 
sources are dominated either by dust (p = 0) or by radiation (p = p/3). 
In the case of the standard model, the particle horizon increases linearly in 
cosmic time (therefore faster than the scale factor). This implies that the 
CMB radiation, today observed with a temperature of 2.7K over the whole 

i^Units Mp = (87rG')"i/2 = 1.72 x IQi^GeV will be used. 

^^The overdot will usually denote derivation with respect to the cosmic time coordinate 
t related to conformal time as dt — a{r])dr]. 
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present horizon has been emitted from space-time regions which were not 
in causal contact. This problem is known as the horizon problem of the 
standard cosmological model. The other problem of the standard model is 
related to the fact that today the intrinsic (spatial) curvature /c/a^ is smaller 
than the extrinsic curvature, i.e. H^. Recalling that k/{a?H'^) = k/a? it 
is clear that if d < 0, 1/a? increases so that the intrinsic curvature could 
be, today, arbitrarily large. The third problem of the standard cosmological 
model is related to the generation of the large entropy of the present Universe. 
The solution of the kinematical problems of the standard model is usually 
discussed in the framework of a phase of accelerated expansion |lU3j . i.e. 
d > and d > 0. In the case of inflationary dynamics the extension of the 
causally connected regions grows as the scale factor and hence faster than 
in the decelerated phase. This solves the horizon problem. Furthermore, 
during inflation the contribution of the spatial curvature becomes very small. 
The way inflation solves the curvature problem is by producing a very tiny 
spatial curvature at the onset of the radiation epoch taking place right after 
inflation. The spatial curvature can well grow during the decelerated phase 
of expansion but is will be always subleading provided inflation lasted for 
sufficiently long time. In fact, the minimal requirement in order to solve 
these problems is that inflation lasts, at least, 60-efolds. The final quantity 
which has to be introduced is the Hubble radius H^^{t). This quantity is 
local in time, however, a sloppy nomenclature often exchanges the Hubble 
radius with the horizon. Since this terminology is rather common, it will also 
be used here. In the following applications it will be relevant to recall some of 
the useful thermo dynamical relations. In particular, in radiation dominated 
Universe, the relation between the Hubble parameter and the temperature 
is given by 



where is the effective number of relativistic degrees of freedom at the 
corresponding temperature. Eq. ()5.2|) . implying H'^Mp = p/3, has been 
used in Eq. fl5.7|l together with the known relations valid in a radiation 
dominated background 




(5.7) 




n{T) 



C(3) 




(5.8) 
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where ({3) ~ 1.2. In Eq. (|5.8|) the thermodynamical relation for the number 
density n{T) has been also introduced for future convenience. 

5.1.1 Inflationary dynamics and its extensions 

The inflationary dynamics can be realized in different ways. Conventional 
inflationary models are based either on one single inflaton field |lU4HlU5lllU6| 
llU7j or on various fields |lU8j (see |lU9j for a review). Furthermore, in the 
context of single-field inflationary models one oughts to distinguish between 
small-field models |lU6j (like in the so-called new-inflationary models) and 
large-field models |lU7j (like in the case of chaotic models). 

In spite of their various quantitative differences, conventional inflationary 
models are based on the idea that during the phase of accelerated expansion 
the curvature scale is approximately constant (or slightly decreasing). After 
inflation, the radiation dominated phase starts. It is sometimes useful for 
numerical estimates to assume that radiation suddenly dominates at the end 
of inflation. In this case the scale factor can be written as 

7] ^ 



^iiv) = , V < -Vu 

ariji) = ^""^ V > -Vi, (5-9) 

where a is some effective exponent parameterizing the dynamics of the pri- 
mordial phase of the Universe. Notice that if a = 1 we have that the pri- 
mordial phase coincides with a de Sitter inflationary epoch. The case a = 1 
is not completely realistic since it corresponds to the case where the energy- 
momentum tensor is simply given by a (constant) cosmological term. In 
this case the scalar fluctuations of the geometry are not amplified and the 
large-scale angular anisotropy in the CMB would not be reproduced. The 
idea is then to discuss more realistic energy-momentum tensors leading to 
a dynamical behaviour close to the one of pure de Sitter space, hence the 
name quasi-de Sitter space-times. Quasi-de Sitter dynamics can be realized 
in different ways. One possibility is to demand that the inflaton slowly rolls 
down from its potential obeying the approximate equations 



3 



dV 

3i70~- — , (5.10) 

0(p 
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valid provided ei = —H/H"^ <^ 1 and 62 = (p/{Hip) <^ 1. There also exist 
exact inflationary solutions like the power-law solutions obtainable in the 
case of exponential potentials: 



2 f 



V = Voe V-^^p, 



2pMp 
~t ' 



a(t)~tP. (5.11) 

Since, from the exact equations, 2MpH = —ifP' the two slow-roll parameters 
can also be written as 

In the case of the exponential potential ()5.11|) the slow-roll parameters are 
all equal, ei = £2 = 1/p. Typical potentials leading to the usual inflation- 
ary dynamics are power-law potentials of the type V{ip) ^ y^", exponential 
potentials, trigonometric potentials and nearly any potential satisfying, in 
some region of the parameter space, the slow-roll conditions. 

Inflation can also be realized in the case when the curvature scale is 
increasing, i.e. H > and a > 0. This is the case of superinfiationary 
dynamics. For instance the propagation of fundamental strings in curved 
backgrounds |lim llllj may lead to superinfiationary solutions |113| I114j . A 
particularly simple case of superinfiationary solutions arises in the case when 
internal dimensions are present. 

Consider a homogeneous and anisotropic manifold whose line element can 
be written as 

ds"^ = G ^,ydx^ dx" = a?{ri)[dif - '^ijdx'dx^] - b^{i])'yabdy''dy'', 
fi,^ = 0, D — 1 = d + n, i,j = l,...,d, a,b = d + 1, ...,d + n. (5.13) 

[t] is the conformal time coordinate related, as usual to the cosmic time t = 
J a{vi)drj ; 7jj(x), Jabiy) are the metric tensors of two maximally symmetric 
Euclidean manifolds parameterized, respectively, by the "internal" and the 
"external" coordinates {x*} and {|/°}]. The metric of Eq. (j5.13|) describes 
the situation in which the d external dimensions (evolving with scale factor 
a{r])) and the n internal ones (evolving with scale factor b{ri)) are dynamically 
decoupled from each other |115j . 
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A model of background evolution can be generically written as 



a{r]) = aJ-—) , b{ri)=bi(~—) 
a{r]) = ai J^!LL^ j ^ b{7]) = bi, 



X 



V < -Vi, 




-7]i<7]< T]2, 



r]>r]2. (5.14) 



In the parameterization of Eq. ()5.14|) the internal dimensions grow (in con- 
formal time) for A < and they shrink for A > 

Superinflationary solutions are also common in the context of the low- 
energy string effective action |119t 11201 1121] . In critical superstring theory 
the electromagnetic field F^^ is coupled not only to the metric (g^u), but also 
to the dilaton background (if). In the low energy limit such interaction is 
represented by the string effective action |119[ I12(J1 1121j . which reads, after 
reduction from ten to four expanding dimensions, 



were = — In Vg = ln((y'^) controls the tree-level four-dimensional gauge 
coupling (0 being the ten-dimensional dilaton field, and Vq the volume of the 
six-dimensional compact internal space). The field a is the Kalb-Ramond 
axion whose pseudoscalar coupling to the gauge fields may also be interesting. 

In the inflationary models based on the above effective action |122[ I123j 
the dilaton background is not at all constant, but undergoes an accelerated 
evolution from the string perturbative vacuum {ip = — oo) towards the strong 
coupling regime, where it is expected to remain frozen at its present value. 
The peculiar feature of this string cosmological scenario (sometimes called 
pre-big bang |123j ) is that not only the curvature evolves but also the gauge 
coupling. Suppose, for the moment that the gauge fields are set to zero. The 

^^To assume that the internal dimensions are constant during the radiation and matter 
dominated epoch is not strictly necessary. If the internal dimensions have a time variation 
during the radiation phase we must anyway impose the BBN bounds on their variation 
IllTl I118j . The tiny variation allowed by BBN implies that 6(77) must be effectively 
constant for practical purposes. 




(5.15) 
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phase of growing curvature and dilaton coupling {H > 0, ip > 0), driven 
by the kinetic energy of the dilaton field, is correctly described in terms 
of the lowest order string effective action only up to the conformal time 
7] = r]s sX which the curvature reaches the string scale Hg = = 
is the fundamental length of string theory). A first important parameter 
of this cosmological model is thus the value (fig attained by the dilaton at 
f] = rjs. Provided such value is sufficiently negative (i.e. provided the coupling 
g = e"^/^ is sufficiently small to be still in the pcrturbative region dXi] = rjs), it 
is also arbitrary, since there is no perturbative potential to break invariance 
under shifts of (p. For rj > rjs high-derivatives terms (higher orders in a') 
become important in the string effective action, and the background enters 
a genuinely "stringy" phase of unknown duration. An assumption of string 
cosmology is that the stringy phase eventually ends at some conformal time r]i 
in the strong coupling regime. At this time the dilaton, feeling a non-trivial 
potential, freezes to its present constant value (p — (pi and the standard 
radiation-dominated era starts. The total duration rji/rjs, or the total red- 
shift Zs encountered during the stringy epoch (i.e. between rj^ and rji), will 
be the second crucial parameter besides Lps entering our discussion. For the 
purpose of this paper, two parameters are enough to specify completely our 
model of background, if we accept that during the string phase the curvature 
stays controlled by the string scale, that is if ~ gMp ~ Xj^ {Mp is the 
Planck mass) for rjs < r) < rji. 

During the string era (p and H are approximately constant, while, during 
the dilaton-driven epoch 

a = (-i)°, q; = — Wl-E, (5.16) 

v3 

Here S = represents the possible effect of internal dimensions whose 

radii bi shrink like {—t)^'- for t — >■ 0_ (for the sake of definiteness we show 
in the figure the case E = 0). The shape of the coupling curve corresponds 
to the fact that the dilaton is constant during the radiation era, that (p is 
approximately constant during the string era, and that it evolves like 

g{v) = e-/' = a\ -^ = ^(3 + ^=^), (5-17) 

during the dilaton-driven era. 

An interesting possibility, in the pre-big bang context is that the exit 
to the phase of decelerated expansion and decreasing curvature takes place 
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without any string tension correction. Recently a model of this kind has been 
proposed |124j . The idea consists in adding a non-local dilaton potential 
which is invariant under scale factor duality. The evolution equations of the 
metric and of the dilaton will then become, in 4 space-time dimensions, 



-2 „ ^t2 



0, 



2(p — ip 



dV 

3H' + V- — 

cnp 



0, 



(5.18) 



where = — 3 log a is the shifted dilaton. A particular solution to these 
equations will be given by 



a{t) 



+ Vr2 + 1 



--log(l + r2) + (^o, 



where 



r 



t 

To' 



tn 



(5.19) 
(5.20) 

(5.21) 
(5.22) 



This solution interpolates smoothly between two self-dual solutions. For 
t ^ — oo the background superinflates while for t +oo the decelerated 
FRW limit is recovered. 



5.2 Gauge fields in FRW Universes 

In a FRW metric of the type (j5.H) the evolution equations of the gauge 
fields are invariant under Weyl rescaling of the four-dimensional metric. This 
property is often named conformal invariance. Consider, for simplicity, the 
evolution equation of sourceless Maxwell fields: 



0. (5.23) 



Since the metric ()5.H) is conformally fiat, the time dependence arising in 
Eq. ()5.23p from y/—G, is always compensated by the controvariant indices 
of the fields strength. Maxwell fields in conformally fiat metrics obey the 
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same equations obtained in the case of Minkowski space-time provided the 
conformal time coordinate is used and provided the fields are appropriately 
rescaled. This property holds not only for Maxwell fields but also for chiral 
fermions in FRW backgrounds. 

Using Weyl invariance, the generalization of MHD equations to the case 
of conformally fiat FRW space-times can then be easily obtained: 

dB ^ ^ 

+ Vx^ = (5.24) 

OT] 

V-^=0, V-i3 = 0, (5.25) 
- - ^ d£ 

VxB = A7rJ + —, (5.26) 
or] 

J = a{£ + vx B), (5.27) 

where 

B = a^B, £ = a^E, J = a^J, a = a^a, (5.28) 

are the curved space fields expressed as a function of the corresponding flat 
space quantities. 

In the case of a radiation dominated background, the Navier-Stokes equa- 
tion takes the form 



> + p)v\' + v-V[{p + p)v\ = -V 



p + 



\B 
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in 



+ ^^-^^ + {p + p)iyV'v, (5.29) 



where the incompressible closure, i.e. V ■ v = 0, has been already adopted 
and where the vector identity V x B x B = —^V\B\'^+ [B ■ V]B has been 
used. In Eq. (I5.29|l p = a^p and p = a^p. In terms of the rescaled pressure 
and energy densities the continuity equation becomes 

1^ + (37 - m-p + (7 + 1) V ■ [-pv] = 0, (5.30) 

where a barotropic equation of state has been assumed for the background 
fluid, i.e. p = 7p. The energy density is homogeneous and hence the last 
term in Eq. ()5.30|) can be consistently dropped within the incompressible 
closure. However, one may ought to linearize Eq. ()5.30|) and this is why the 
last term may be kept for future considerations. If 7 = 1/3, then Eq. ()5.30p 
implies p' = 0. For generic barotropic index Eq. ()5.29|) becomes 

|^ + (l-37)W=#^ + .VV, (5.31) 

dr] ^ ' 47r(7 + l)p ^ ^ 
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where the magnetic pressure term has been dropped since it is neghgible with 
respect to p. 

Recalhng that for L ^ L^j the magnetic diffusivity equation can be writ- 
ten, using vector identities, as 

^ + [v-V]B=[B-V]v + ^V^B. (5.32) 

dT] Ana 

Eqs. (I5.24j) - (j5.29|) form a closed system sharing analogies with their flat 
space counterpart. Thanks to conformal invariance, various solutions can 
be lifted from fiat to curved space like, for instance the following fully 
nonlinear solution of Eqs. ()5.32p and ()5.29p 

v = ± ,^ (5.33) 
VP+P VP 

valid if thermal and magnetic diffusivities are neglected and if p ^ 

Eqs. ()5.24j) - ()5.29j) present also some differences with respect to the flat 
space case: since the Universe is expanding, the relative balance between 
diffusion scales may change as time goes by. All the propagation of MHD ex- 
citations can be then generalized to the curved-space case. In various papers 
jl25| 1126^ I127| I127j this program has been achieved (see also |129^ I13()j for 
earlier attempts). In \^\ the evolution of magnetic fields with non- vanishing 
magnetic helicity has been studied using the resistive MHD approximation 
in curved space. In |131j the evolution of non-linear Alfven waves with the 
purpose of deriving possible observable effects |132j . 

In order to analyze the high-frequency spectrum of plasma excitations, the 
Vlasov-Landau approach can be generalized to curved spaces [133j. There are 
two options in order to achieve this goal. The first option is to assume that 
the magnetic field breaks the isotropy of space-time. This approach leads to 
the study of plasma effects in a magnetic field in full analogy with what is 
done in terrestrial tokamaks. The second choice is to discuss magnetic fields 
which are stochastic and which do not break the isotropy. Later in the present 
review the effects of magnetic fields breaking the isotropy of space-time will 
be considered but, for the moment, the attention will be concentrated on the 
second approach. 

A typical problem which should be discussed within the kinetic approach 
is the relaxation of charge and current density fluctuations in a relativistic 
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plasma |134j . The problem is then to compute what magnetic field is induced 
from charge and current density fiuctuations. 

Consider an equilibrium homogeneous and isotropic conducting plasma, 
characterized by a distribution function 

/o(p) = ^e-^ (5.34) 

common for both positively and negatively charged ultrarelativistic particles 
(for example, electrons and positrons). The normalization is chosen in such 
a way that / d^pfoip) = Ug. Within the kinetic approach, the initial charge 
and current density fiuctuations are related to the perturbations of the initial 
value of the distribution function around the equilibrium distribution ()5.34|) . 
i.e 

f+{x,P,v) = fo{p) + Sf+{x,P,v), 

f-{x,P,v) = fo{p) + Sf-{x,p,r]), (5.35) 

where + refers to positrons and — to electrons, and p is the conformal mo- 
mentum. The Vlasov equation defining the curved-space evolution of the 
perturbed distributions can be written as 

+ v ■ + e{S + v X B) ■ = [^—] , (5.36) 
or] ox op \ or] J coU 

^ + v-^~e{£ + vxB)-^=(^\ (5 37) 

dr] dx dp \ dr] /cou' 

where the two terms appearing at the right hand side of each equation are 
the collision terms. This system of equation represents the curved space 
extension of the Vlasov-Landau approach to plasma fiuctuations jH21 IBS]- 
All particle number densities here are related to the comoving volume. 



Notice that, in general v = p/ ^p^ + m\c?. In the ultra-relativistic limit 
V = p/\p\ and the Vlasov equations are conformally invariant. This implies 
that, provided we use conformal coordinates and rescaled gauge fields, the 
system of equations which we would have in fiat space [HI] looks exactly the 
same as the one we are discussing in a curved FRW (spatially fiat) background 

[ns]. 

The evolution equations of the gauge fields coupled to the two Vlasov 
equations can be written as 



V.£ = 4.e/dV|/,(f,p.„)-/_(f.p.„)l, 
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V X £ + B' = 0, 

V X B-£' = A7re J £pv[f+{x,p,ri) - f_{x,p,ri)]. (5.38) 

In order to illustrate the physical relevance of this description, Eqs. (|5.38|) 
will be used later in this Section to study the fate of charge and current 
density fluctuations. 

5.3 Physical scales in the evolution of magnetic fields 

When MHD is discussed in curved spce-times, the evolution of the geometry 
reflects in the time evolution of the physical scales of the problem. 



5.3.1 Before matter-radiation equality 

From Eq. ()5.32p . the comoving momentum corresponding to the magnetic 
diffusivity scale is given by the relation 

n = aH. (5.39) 



Ana 7] 

In the case of a relativistic plasma at typical temperatures T ^ 1 MeV, the 
fiat-space conductivity (see Eq. ()5.28|) ) is given by 

<yc - ^, (5.40) 

where a ~ al^/T"^ is the cross section Then, it can be shown |135j that 

T 

Oc = (To (5.41) 

where (Tq is a slowly increasing function of the temperature. In |135j it 
has been convincingly argued that ctq ~ 0.06 at T ~ 100 MeV while it 
becomes 0.6 at T ~ 100 GeV (for a recent analysis of conductivity in high 
temperature QED see [136J). Recalling Eq. ()5.28|1 the physical diffusivity 
scale uja- = k„/a{ri) will be, in units of the Hubble parameter. 



H V \Hj ■ " ^* VMi 



72.2x.o^^-/^f^V . (5.42) 
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The cross section is denoted with a while the conductivity is a. 
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Discounting for the effect due to the possible variations of the effective rel- 
ativistic degrees of freedom, the (physical) magnetic diffusivity length scale 
is much smaller than the Hubble radius at the corresponding epoch. For 
instance, if T ~ 10 MeV, then ~ uj~^ ~ 10~^°iJ~^. This means that 
there are ten decades length-scales where the magnetic field spectrum will 
not experience diffusion due to the finite value of the conductivity. 

The estimate of the magnetic diffusivity scale will now be repeated in a 
different range of temperatures, namely for T < 0.2me when weak interac- 
tions have fallen out of thermal equilibrium. In this case the relevant degrees 
of freedom are the three neutrino species and the photon, i.e. = 3.36. In 
this case the conductivity will be given by 

cr ~ (5.43) 



a, = —[ — ]'\, (5.44) 



which leads to 



if t;th ~ {T/m^f/^ and a ~ alJT^. Eq. (lOll leads to the known T^'^ 
dependence of the conductivity from the temperature, as expected from a 
non-relativistic plasma. The physical magnetic diffusivity momentum is, in 
this case 




11 __l/4/T^v4^Mp^l/2 



9:"'[-] [-^] ■ (5.45) 



H v a Kmc- J V H XnipJ \ T 



Eqs. fl5.42|) and ()5.45|) refer to two different ranges of temperatures. The 
magnetic diffusivity scale in units of the Hubble parameter obtained in Eq. 
fl5.45|) is smaller than the one given by in Eq. ()5.42|) . Hence, by lowering 
the temperature across the weak equilibration temperature makes the con- 
ductivity length scale larger in units of the Hubble radius : inertial range 
of length scales (i.e. the region where conductivity effects are negligible) is 
larger above the weak equilibration temperature. Summarizing 

K{T)= \J* — Lh{T), T>MeV 



72.2 VMp. 

. rp s 1/2 / 7^ X -1/4 

Mr)^^(^) (-) L«(T), T<MeV<5.46) 
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In more general terms the inertial range is the interval of scales (either in momentum 
or in real space) where the dynamics is independent on the scales of dissipation so that, 
in this range, the difffusivities can be taken to zero. 
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where Lj^(T) = H-\T). 

In the evolution of the plasma it is sometimes relevant to estimate the 
thermal diffusivity scale setting a bound on the coherence of the velocity 
field and not of the magnetic field. There are however phenomena, like the 
propagation of MHD waves, where the excitations of a background magnetic 
field are coupled to the excitations of the velocity field. Generally speaking 
the typical correlation scale of the velocity field is always smaller, in various 
physical situations, than the typical correlation scale of the magnetic fields. 
The range of physical momenta affected by thermal diffusivity is bounded 
from above by the thermal diffusivity scale which can be obtained from the 
following relations 

fcV ~ —L(t]) c^- = n = aH, (5.47) 
5 1] 

were is the photon mean free path, i.e. 

i-y ~ (5.48) 

UeCT 

The mean free path measures the efficiency of a given process to transfer 
momentum. From Eqs. ()5.29j) and ()5.47j) the typical thermal diffusivity 
momentum is: 

where Udis = ^difr/o is the physical momentum as opposed to the comov- 
ing momentum kdm- The scale defined in Eq. fl5.49p is sometimes called 
Silk scale. Notice that the mean free path given in Eq. ()5.48|) increases as 
a^{rj). Since the physical scale of an inhomogeneous region increases as a{rj), 
the photon mean free path can become larger than the physical size of an 
inhomogeneous region. This is the so-called free streaming regime. 

If T < 0.2me, the photon mean free path is given by Eq. ()5.48|) where 
ax = 1.7 X lO^GeV"^ ~ aem^sT^ is the Thompson cross section and 

ne ~ 6.6 X 10"^(fib/io) a^cT^ (5.50) 

is the electron number density depending both on the critical fractions of 
baryons, f2b and on the ionization fraction Xe- The ionization fraction is the 
ratio between the number density of protons and the total number density 
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of baryons. The value of Xg drops from 1 to 10~^ in a short time around 
recombination. Eq. (j5.5(J|) is simply derived by noticing that, after electron 
positron annihilation, Ue — rjn^ where r/ = 2.74 x 10~^f2b/z.o is the baryon- 
to-photon ratio, the free parameter of big-bang nucleosynthesis calculations 
[138, 1321 Enillll]. Thus 

^(T).8,8xl0-(^) (^) {—)t-\ (5,51) 

For temperatures larger than the MeV the photon mean free path can 
be estimated by recalling that, above the MeV, the Thompson cross section 
is replaced by the Klein-Nishina cross section (i.e. a ~ T~^GeV~^). Thus, 
from Eq. (j5.48|) . the photon mean free path is 

£,(T) ^ 0.1 GeV-^ (5.52) 

As before, it is useful to refer ujais to the Hubble rate and to compute the 
typical diffusion distances for T > 1 MeV 

and for T < 1 MeV 

^S(r) 1.03 X l0-f^V"Y^V"V^V'\/^. (5.54) 



Lh{T) ■ VO.02/ VO.5/ VlO.75/ V T 

At temperatures (9(MeV), corresponding roughly to e"*" e~ annihilation, Eqs. 
fl5.53|) and ()5.54|) differ by a factor 10^. This drop of five orders of magnitude 
of the thermal diffusivity scale reflects the drop (by ten orders of magnitude) 
of the photon mean free path as it can be appreciated by comparing Eqs. 
fICTll and at T ^ O(MeV). 

For temperatures smaller than the electroweak temperature (i.e. T ~ 100 
GeV) down to the temperature of neutrino decoupling, neutrinos are the 
species with the largest mean free path and, therefore, the most efficient 
momentum transporters. Neglecting possible variations of the number of 
relativistic species the neutrino mean free path will be simply given by ^y = 
Gp^T^^(niep + nqua)^^ where niep and Uq^n are the leptons and quark number 
densities. Using this estimate the diffusion scale will be 

- 10-(^)^'^« (5.55) 
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which should be compared with Eq. ()5.53|) . Clearly L^j'g > L^^g. 

It seems then difficult to get observable effects from inhomogeneities in- 
duced by MHD excitations. In |137j the damping rates of Alfven and magne- 
tosonic waves have been computed and it has been shown that MHD modes 
suffer significant damping (see, however, |13Hll32j ). To avoid confusion, the 
results reported in Ref. |137j apply to the damping of MHD waves and not 
to the damping of magnetic fields themselves as also correctly pointed out in 



5.3.2 After matter-radiation equality 

In standard big bang cosmology there are three rather close (but rather 
different) epochs, 

• the matter-radiation equality taking place at Teq ^ 5.5 {Qq h^) eV, i.e. 



• the decoupling of radiation from matter occurring at Tdec — 0.26 eV 



• the recombination occurring, for a typical value of the critical baryon 
fraction, at T^ec ^ 0.3 eV, i.e. tree ^ 4.3 x lO^^j^^o K)'^'^ sec. 
After plasma recombines, the baryons are no longer prevented from moving 
by viscosity and radiation pressure. Thus, the only term distinguishing the 
evolution of the baryons from the evolution of a cold dark matter fiuid is the 
presence of the Lorentz force term in the Navier-Stokes equation. In this sit- 
uation, it is interesting to investigate the interplay between the forced MHD 
regime and the generation of gravitational inhomogeneities. It was actually 
argued long ago by Wasserman |143j and recently revisited by Coles |144j 
that large-scale magnetic fields may have an important impact on structure 
formation. Assuming that, for illustrative purposes, the matter fiuid consists 
only of baryons and cold dark matter, for t > tree the evolution of the system 
is described by the following set of non relativistic equations 



ma. 



4 X 10^° (fio hl)-^ sec; 



i.e. tdec ^ 5.6 X 10^2 (fio /l2)-l/2 



sec 




-V$ + 



V xBxB 
47ra4(t)pb ' 



(5.56) 



(5.57) 



= 7 

2M^ 



(5.58) 



56 



where $ is the Newtonian potential. Eqs. (|5.56|) come from Eq. (|5.31|) 
for 7 = in the case of homogeneous pressure density. Eq. (|5.57p can be 
obtained, in the same hmit, from Eq. Eq. ()5.3H) by setting to zero the 
Lorentz force. 

From Eq. ()5.30|) . by hnearizing the continuity equations for the two 
species, the relation between the divergence of the velocity field and the 
density contrast can be obtained 

5'y, = -V- tJb, 5c = -V ■ Vc, (5.59) 

where (5b = 5pb/Pb and 6c = Spc/ Pc Thus, taking the divergence of Eqs. ()5.56p 
and (|5.57p and using, in the obtained equations, Eqs. (|5.58|) and (|5.59|) the 
following system can be obtained 

.// . V ■ [(V X i3) X i3] 

K + - ^PmS = 0, (5.61) 

where 6 = fib^b + ^c^c and Pm = Pb + Pc is the total matter density. Multi- 
plying Eq. ()5.60|) by fib and Eq. ()5.6Hl by Qc, the sum of the two resulting 
equations can be written as 

where it has been used that the background geometry is spatially flat. Going 
to cosmic time, Eq. ()5.62|) reads 

^■,l,-A,,V-|(Vxflxgl. (5.63) 

3t 3^2 47rpba6(t) ^ ^ 

The homogeneous term of Eq. ()5.63|) can be easily solved. The source term 
strongly depends upon the magnetic field configuration [HH] ■ However, as far 
as time evolution is concerned, B is constant in the inertial range of scales, 
i.e. L > Lfj. Hence, B = Bq and B{t) = Bo/a^it). Assuming a stochastic 
magnetic field the solution of equation ()5.63|) can be obtained integrating 
from tree and recalling that Phif) ~ Pb(^roc)a'~^- The result is 



J.2 I 7^ 

5{t) = fib- ' 



4vrpb(trcc)a^(^rec)A 



9ft ^2/3 3^^^^^^ 3 



lOVtrnc/ ^ 5 \ t J 2 



(5.64) 
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where Ab is the typical length scale of variation of the magnetic field. The 
total density contrast grows as t^^^ with an initial condition set by the Lorentz 
force, i. e. 

S{t.ec) ^ n^-, If'fl}; (5-65) 

To seed galaxy directly the effective 5(trec) should be of the order off 10^^ 
for a fiat Universe, implying that 

|4ec| ^ G. (5.66) 

Also the converse is true. Namely one can read Eq. ()5.66|) from left to right. 
In this case Eq. ()5.66|) will define the typical scale affected by the presence 
of the magnetic field. This scale is sometimes called magnetic Jeans scale 
[101 : it can be constructed from the Jeans length 

Aj = Cs , (5.67) 



where Cg is replaced by Ca, i.e. the Alfven velocity. 



5.4 Charge and current density fluctuations 

As a simple application of kinetic techniques, it is useful to understand which 
magnetic field can be obtained from an initial inhomogeneous distribution of 
charge and current density fluctuations. 

Consider then the system of Eqs. (j5.38p under the assumption that, 
initially, electric and magnetic fields vanish: 

Sf±{x,P,'no) 0, £{x,r]o) = B{x,r]o) = 0. (5.68) 

Using Eqs. ()5.36|) - ()5.38|) the induced magnetic fields can be obtained.. This 
problem is discussed, in greater detail, in jl34j . 
Subtracting Eqs. (|5.36p and (j5.37|) . and defining 

f{x,p,7]) =6U{x,p,ri) -6f_{x,p,7]) (5.69) 

the Vlasov-Landau system can be written as 

-9dP) + ^ik■v- uj)f^,Jp) + 2e4 ■ ^ = -z>4^, (5.70) 
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tk . £^,^ = Ane J f^,Jp)d'p, (5.71) 
2fc-i%^ = 0, (5.72) 
4 = -k X 4^., (5.73) 

- + -Mk ■ 4) = 4vre / d'p v f^,Jp}, (5.74) 

where ^^(p) is the initial profile of the distribution function and where z/ is the 
typical collision frequency. The subscripts ku denote the Fourier transform 
with respect to space and the Laplace transform with respect to time of the 
corresponding quantity. 

After enforcing the Gauss constraint at the initial time, the electric field 
can be separated in its polarizations parallel and transverse to the direction 
of propagation of the fluctuation. The transverse current provides a source 
for the evolution of transverse electric field fluctuations 

tuil - ^)£l = Ane J d'pf^^Jpjv^, (5.75) 

whereas the charge fluctuations provide a source for the evolution of longi- 
tudinal electric field fluctuations 

tk-4=A7ie I d'pfrip). (5.76) 



'few J fct^ 

In Eqs. fl5.75j) and ()5.76j) the longitudinal and transverse part of the electric 
field fluctuations have been defined: 

4 = 4-^(^-4), 4. = ^(^1-^)- (5-77) 

The solution of Eq. ()5.70p is given by 

hJp) = — -Mp) - 2^-^-4^]. (5.78) 

i[k ■ V — UJ — iv) op 

where we used that dfo/dp = vdfo/dp. The longitudinal and transverse 
components of the electric fluctuations can be obtained by inserting Eq. 
(Kim into Eqs. (IK7H|) - (fH7?)|) 

fe-' tk ey J "^iik-v-uj-iv) 
St = / d'pv-^^ , (5.80) 
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where ey and e± are, respectively, the longitudinal and transverse part of the 
polarization tensor 

e,ik,.)^l-'-^[d^^^ f , (5.81) 

e±{k,co) = 1 - - / d'p-^^^^—-^. (5.82) 
^ J {k ■ V - 00 - tu) op 

Now, the general expression for the generated magnetic field is 

Br = ^ ^ / d^plvxk]^-^ . (5.83) 

'^'^ u;h^{k,u;)-k^J \k-v-u-iv) ^ ' 

The space-time evolution of the magnetic fluctuations can be determined 
by performing the inverse Laplace and Fourier transforms: 

B{x,r^) = f ^ ""^^^ ,^ / d'ke^'-'lvxk] f d'p-^^ -. 

J uj^e±{k,uj) - k^ J J . k - CO - 

(5.84) 

In order to perform this integral, the explicit relations for the polariza- 
tion tensors should be given. They depend on the equilibrium distribution 
function fo{p), which we take to be the expression given in Eq. ()5.34j) ^^ 

Then we have for transverse polarization 

and for the longitudinal polarization: 

<„ft.) = l + ^{2 + — (5.86) 

Consider now the case of very small momenta k <^uj and ^ z/, relevant 
for long-ranged magnetic fields. Then the computation of the integral ()5.84p 
in the large time limit and with the use of explicit form of the transverse 
polarization tensor in (j5.85p gives ^^: 

^ T ^ ^ 

B{x, 7]) ~ exp (-fc^r^/fj) V x J, (5.87) 

47raem'^g 



^""^ Notice that most of our considerations can be easily extended to the case of a Bose- 
Einstein or Fermi-Dirac distribution. What is important, in our context, is the analytical 
structure of the polarization tensors and this is the same for different distributions jl25| . 

^^For small k, the equation uj'^e±{k, lj) — fc^ = defining the poles of the inverse Laplace 
transform implies uj ^ ik'^ /a. 
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where a is the plasma conductivity in the relaxation time approximation, 



o = (5.88) 
and initial electric current is given by 

J{x) = f d'pvg^ip) . (5.89) 



The obtained results assumed that the linearization of the Vlasov equa- 
tion is consistent with the physical assumptions of our problem. This is 
indeed the case. In order to safely linearize the Vlasov equation we have to 
make sure that the perturbed distribution function of the charge fluctuations 
is always smaller than the first order of the perturbative expansion (given by 
the distribution of Eq. l5.H4j) . In other words we have to make sure that 

|5/+(f,p,r/)|«/o(p), |5/_(f,p,r7)|«/o(p). (5.90) 

These conditions imply that 

e^o. . 9 Hp) 
\k ■ V — iu — iD\ dp 

If the relativistic plasma frequency is defined, 



« /o(p). (5.91) 



= (6.92) 

the condition expressed by Eq. ()5.91|) can be restated, for modes k < Up, as 
< UqT (where we essentially took the square modulus of Eq. ()5.91|l ). 
This last inequality expresses the fact that the energy density associated with 
the gauge field fluctuations should always be smaller than the critical energy 
density stored in radiation. The linear treatment of the Vlasov equation is 
certainly accurate provided the typical modes of the the field are smaller 
than the plasma frequency and provided the energy density in electric and 
magnetic fields is smaller than T^, i.e. the energy density stored in the 
radiation background. 
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5.5 Inverse cascades 

A relevant topic in the evolution of large scale magnetic fields concerns the 
concepts of direct and inverse cascades. A direct cascade is a process, oc- 
curring in a plasma, where energy is transferred from large to small length 
scales. An inverse cascade is a process where the energy transfer goes from 
small scales to large length scales. One can also generalize the the concept 
of energy cascade to the cascade of any conserved quantity in the plasma, 
like, for instance, the helicity. Thus, in general terms, the transfer process of 
a conserved quantity is a cascade. 

The concept of cascade (either direct or inverse) is related with the con- 
cept of turbulence, i.e. the class of phenomena taking place in fluids and plas- 
mas at high Reynolds numbers, as anticipated in Section 4. Non-magnetized 
fluids become turbulent at high kinetic Reynolds numbers. Similarly, magne- 
tized fluids should become turbulent for sufficiently large magnetic Reynolds 
numbers. The experimental evidence of this occurrence is still poor. One of 
the reasons is that it is very difficult to reach, with terrestrial plasmas, the 
physical situation where the magnetic and the kinetic Reynolds numbers are 
both large but, in such a way that their ratio is also large i.e. 

R 

R„,>1, R>1, Pr = -^>1. (5.93) 

R 

The physical regime expressed through Eqs. ()5.93|1 rather common in the 
early Universe. Thus, MHD turbulence is probably one of the key aspects of 
magnetized plasma dynamics at very high temperatures and densities. Con- 
sider, for instance, the plasma at the electroweak epoch when the temperature 
was of the order of 100 GeV. One can compute the Reynolds numbers and 
the Prandtl number from their definitions given in Eqs. ()4.40|1 - ()4.42|1 . In 
particular, 

Rm ~ 10^^ R = 10^\ Pr ~ 10^ (5.94) 

which can be obtained from Eqs. ()4.40|) - ()4.42j) using as fiducial parameters 
V ~ 0.1, aT/a, v ^ (aT)"^ and L ~ 0.01 E'^ ~ 0.03 cm for T ~ 100 GeV. 

At high Reynolds and Prandtl numbers the evolution of the large-scale 
magnetic fields is not obvious since the inertial regime is large. There are 
then two scales separated by a huge gap: the scale of the horizon (3 cm 
at the electroweak epoch) and the diffusivity scale (about 10"*^ cm). One 
would like to know, given an initial spectrum of the magnetic fields, what 
is the evolution of the spectrum, i.e. in what modes the energy has been 
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transferred. The problem is that it is difficult to simulate systems with huge 
hierarchies of scales. 

Suppose, for instance, that large scale magnetic fields are generated at 
some time to and with a given spectral dependence. The energy spectrum 
of magnetic fields may receive, for instance, the dominant contribution for 
large comoving momentum k (close to the cut-off provided by the magnetic 
diffusivity scale). The "initial" spectrum at the time is sometimes called 
injection spectrum. The problem will then be to know the how the spectrum 
is modified at a later time t. If inverse cascade occurs, the magnetic energy 
density present in the ultra-violet modes may be transferred to the infra-red 
modes. In spite of the fact that inverse cascades are a powerful dynamical 
principle, it is debatable under which conditions they may arise. 

The possible occurrence of inverse cascades in a magnetized Universe has 
been put forward, in a series of papers, by Brandenburg, Enqvist and Olesen 
[145, 11461 11471 1148j (see also |J;49J for an excellent review covering also in a 
concise way different mechanisms for the generation of large-scale magnetic 
fields). 

MHD simulations in 2 + 1 dimensions |1451 1146j seem to support the idea 
of an inverse cascade in a radiation dominated Universe. Two-dimensional 
MHD simulations are very interesting but also peculiar. In 2 + 1 dimen- 
sions, the magnetic helicity does not exist while the magnetic flux is well 
defined and conserved in the ideal MHD limit. It is therefore impossible to 
implement 2 + 1-dimensional MHD simulations where the magnetic flux and 
the magnetic helicity are constrained to be constant in the ideal limit. This 
technical problem is related to a deeper physical difference between magne- 
tized plasmas with and without helicity. Consider then, for the moment, the 
situation where the helicity is strictly zero. It is indeed possible to argue 
that, in this case, 2 + 1-dimensional MHD simulations can capture some of 
the relevant features of 3 + 1-dimensional MHD evolution. 

In |147t I148l| the scaling properties of MHD equations have been investi- 
gated in the force-free regime (i.e. in the approximation of negligible Lorentz 
force). In the absence of helical components at the onset of the MHD evolu- 
tion a rather simple argument has been proposed by Olesen |147j in order to 
decide under which conditions inverse cascade may arise. Consider the mag- 
netic diffusivity and Navier-Stokes equations in the limit where JT" x ;B ~ 0. 
Then, the evolution equations are invariant under the following similarity 
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transformations: 



X — > Lx, 7] ^ L ^rj, -u — > L'^v, 

B^L'B, z/^L^+V, a^L-^-'a. (5.95) 

For simplicity it is also useful to concentrate the attention on the inertial 
range of momenta where the dynamics is independent on the scale of dissi- 
pation. 

In order to understand correctly the argument under discussion the vol- 
ume average of the magnetic energy density has to be introduced (see |15Uj 
for a lucid discussion on the various possible averages of large-scale magnetic 
fields and see [151J for a complementary point of view). A generic rotation- 
ally and parity invariant two-point function for the magnetic inhomogeneities 
can be written as 

g{f,r]) = {B{x,T])B{y,r])), f = x - y. (5.96) 

Defining now 

g{r,ri) = J £kg{k,ri)e-'''-\ (5.97) 
From Eqs. (|5.96j) - (j5.97|) the following relation can be obtained 

V J dkQk{ri) = J d'xd'y{B{x,7])B{y,r])), (5.98) 

where Qk{v) is related to g{k, rj) and V is the averaging volume. 

The idea is now to study the scaling properties of Qkiv) under the sim- 
ilarity transformations given in Eq. ()5.95|) . Following |147j . it is possible to 
show that 

a\v)Qkiv) = k'^Fix) 

X = k^^+'^^/^T]. (5.99) 

with a = — 1 — 2e. In Eq. (|5.99p Qkiji) is related to the magnetic energy 
density and F{x) is a function of the single argument x = {k^^^°'^/'^7]). If, 
initially, the spectrum of the magnetic energy density is A;", then at later time 
the evolution will be dictated by F[x). Since a^{rii)Qk{rii) is approximately 
constant, because of flux conservation, then the comoving momentum will 
evolve, approximately, as 

k ~ ^-2/(3+a) ^ ^-l/(3+a)^ (5.100) 
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while the physical momentum will scale like 

UO ^ ^ ^-(5+a)/(2a+6)_ (5.101) 

The conclusion of this argument is the following: 

• if q; < — 3 the cascade will be direct (or forward); 

• if a > —3 the comoving momentum moves toward the infra-red and 
the cascade will be inverse; 

• if a = —3, then, e = 1 and, from Eqs. ()5.95|) - ()5.99|) . the momentum 
and time dependence will be uncorrelated. 

The Olesen scaling argument has interesting implications: 

• if the initial spectrum of the magnetic energy density is concentrated 
at high momenta then, an inverse cascade is likely; 

• from Eqs. ()5.100|) and ()5.101|) the typical correlation scale of the mag- 
netic field may, in principle, evolve faster than in the case where the 
inverse cascade is not present. 

In a FRW Universe the correlation scale evolves as ^{r|) = i{rio)a{'r]) from 
a given initial time rjo where magnetic energy density is injected with spectral 
dependence fc". If a > -3, Eqs. (|5.1(l(ljl - (j5.101jl imply that 

2 a + 5 

= ^iVo)a{v) (-) ^ ^iVo) (^) (5.102) 

where the second equality follows in a radiation dominated phase of expan- 
sion. In flat space, the analog of Eqs. ()5.100|) - ()5.101|) would imply that 
i{t) ~ i{to){t/to)^/^°+^\ Finally, in the parametrization of Olesen [147^ the 
case of a Gaussian random field corresponds to a = 2, i.e. e = —3/2. In this 
case the correlation scale will evolve as 

i{r])=£{r]o)airi)(^Y\ (5.103) 

In the context of hydrodynamics the scaling law expressed through Eq. ()5.103p 
is also well known [153J. It is appropriate, at this point, to recall that Hogan 
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[152j was the first one to discuss the evolution of the correlation scale pro- 
duced by a Gaussian random injection spectrum (i.e. a = 2). Hogan was 
also the first one to think of the first order phase transitions as a possible 
source of large scale magnetic fields |152j . 

An argument similar to the one discussed in the case of MHD can be 
discussed in the case of ordinary hydrodynamics. In this case, the analog of 
Qk{v) will be related to the kinetic energy density whose spectrum, in the 
case of fully developed turbulence goes as i.e. the usual Kolmogorov 

spectrum. 

The evolution of the correlation scale derived in Eq. ()5.103|) can be ob- 
tained through another class of arguments usually employed in order to dis- 
cuss inverse cascades in MHD: the mechanism of selective decay jHZl which 
can be triggered in the context of magnetic and hydrodynamical turbulence. 
The basic idea is that modes with large wave-numbers decay faster than 
modes whose wavenumber is comparatively small. Son jl42j applied these 
considerations to the evolution of large-scale magnetic fields at high temper- 
atures and his results (when the magnetic field configuration has, initially, 
vanishing helicity) are in agreement with the ones obtained in |147t I148j . In 
|154j the topics discussed in |147| I142j have been revisited and in |155j the 
possible occurrence of inverse cascades has been criticized using renormaliza- 
tion group techniques. 

The main assumption of the arguments presented so far is that the evo- 
lution of the plasma occurs in an inertial regime. As discussed above in this 
Section, at high temperatures the inertial regime is rather wide both from 
the magnetic end thermal points of view. However, in the case where the 
magnetic and thermal diffusivity scales are finite the arguments relying on 
the inertial regime cannot be, strictly speaking, applied. In this situation 
numerical simulations should be used in order to understand if the diffusive 
regime either precedes or follows the inverse cascade. 

If the magnetic and thermal diffusivity scales are finite, two conceptually 
different possibilities can be envisaged. One possibility is that the inverse 
cascade terminates with a diffusive regime. The other possibility is that 
the onset of diffusion prevents the occurrence of the inverse cascade. This 
dilemma can be solved by numerical simulations along the lines of the ones 
already mentioned |1451 114b'] . The numerical simulations presented in |145| 
I146j are in (2 + 1) dimensions. The results support the occurrence of inverse 
cascade in the presence of large kinetic viscosity. The moment at which 
the cascade stops depends on the specific value of the viscosity coefficient. 
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Figure 5: The time evolution of the correlation scale is illustrated in the 
case of a FRW Universe dominated by radiation and in units io = 1. The 
case a = 2 corresponds to the behaviour obtained in (j5.1U3|) . The case 
a = 2, corresponding to Eq. (j5.1(J4|) . refers to the situation of inverse cascade 
with an initial helical spectrum. For comparison the case a = —3/2 is 
also reported. This last evolution takes place if the injection spectrum is 
of Iroshnikov-Kraichnan type. 

In jl56| I157j (3 + l)-dimensional MHD simulations have shown evidence of 
inverse cascade. Another possible limitation is that in the early Universe 
the magnetic and kinetic Reynolds numbers are both large and also their 
ratio, the Prandtl number, is very large. This feature seems to be difficult 
to include in simulations. 

The considerations presented so far do not take into account the possi- 
bility that the initial magnetic field configuration has non vanishing helicity. 
If magnetic helicity is present the typical features of MHD turbulence can 
quahtatively change |HZI and it was argued that the occurrence of inverse cas- 
cades may be even more likely |158| I159j (see also |16Uj ) . The conservation 
of the magnetic fiux is now supplemented, in the magnetic sector of the evo- 
lution, by the conservation of the helicity obtained in Eqs. ()4.25|) and ()4.27|) . 
The conservation of the helicity is extremely important in order to derive 
the MHD analog of Kolmogorov spectrum. In fact, the heuristic derivation 
of the Kolmogorov spectrum is based on the conservation of the kinetic en- 
ergy density. If the conservation of the magnetic helicity is postulated, the 
same qualitative argument leading to the Kolmogorov spectrum will lead 
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in the case of fully developed MHD turbulence to the so-called Iroshnikov- 
Kraichnan spectrum [S7| ^1 , 162] whose specific form, is, in terms of the 
quantity introduced in Eq. ()5.99|) Qkivo) ~ k~^/'^. 

Suppose, for simplicity, that the injection spectrum has a non-vanishing 
helicity and suppose that the helicity has the same sign everywhere. In this 
case A-B ~ i\B\'^. Since the helicity is conserved the magnetic field will scale 
according to Eqs. (j5.95j) but with e = — 1/2 corresponding to a = 0. Hence, 
the correlation scale will evolve, in this 

2 5 

l{r,) = i{r^o)a{v) (-) ' ^ ^ivo) (^Y, (5.104) 

i.e. faster than in the case of Eq. (j5.103p . It is difficult to draw general 
conclusions without specifying the nature and the features of the mechanism 
producing the injection configuration of the magnetic field. The results crit- 
ically summarized in this Section will be important when discussing specific 
ideas on the origin of large-scale magnetic fields. In Fig. El the cases of dif- 
ferent injection spectra are compared. It is clear that an injection spectrum 
given by a Gaussian random field leads to a rather inefficient inverse cas- 
cade. Since the magnetic helicity vanishes at the onset of the evolution, the 
rate of energy transfer from small to large scales is rather slow. On the con- 
trary if the injection spectrum is helical, the rate increases significantly. For 
comparison the situation where the injection spectrum corresponds to the 
Iroshnikov-Kraichnan case, i.e. the case of fully developed MHD turbulence, 
is also reported. The dashed line of Fig. illustrates the evolution of the 
correlation scale dictated by the the Universe expansion. 

5.6 Hypermagnetic fields 

At small temperatures and small densities of the different fermionic charges 
the SUl{'2) ® Uy{X) is broken down to the f/em(l) and the long range fields 
which can survive in the plasma are the ordinary magnetic fields. For suf- 
ficiently high temperatures (and for sufficiently high values of the various 
fermionic charges) the SUl{2) (8> ^y(l) is restored and non-screened vector 
modes correspond to hypermagnetic fields. In fact, Abelian electric fields 
decay within a typical time scale l/cr where a is the conductivity. The long- 
ranged non- Abelian magnetic fields (corresponding, for instance, to the color 
SU{3) or to the weak SU{2)) cannot exist because at high temperatures 
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the non-Abelian interactions induce a "magnetic" mass gap ~ g'^T where g 
is the gauge couphng constant. Also the non-Abehan electric fields decay 
because of the finite value of the conductivity as it occurs for Abelian elec- 
tric fields. Therefore, the only long scale field that can exist in the plasma 
for enough time must be associated with some Abelian U(l) group. This 
statement, valid to all orders in perturbation theory, has been confirmed 
non-perturbatively for the electroweak theory by recent lattice studies in 
|165j . Under normal conditions (i.e. small temperatures and small densities 
of the different fermionic charges) the SU(2) xU(l)y symmetry is "broken" 
down to U(1)£;m, the massless field corresponding to U(1)em is the ordinary 
photon and the only long-lived field in the plasma is the ordinary magnetic 
one. At sufficiently high temperatures, T > Tc, the SU(2)xU(l)y symmetry 
is "restored", and non-screened vector modes correspond to the U(l)y 
hypercharge group. Hence, if primordial fields existed at T > Tc, they did 
correspond to hypercharge rather than to U(1)£;m- 

At the electroweak epoch the typical size of the Hubble radius is of the or- 
der of 3 cm . The typical diffusion scale is of the order of 10~^ cm. Therefore, 
over roughly eight orders of magnitude hypermagnetic fields can be present 
in the plasma without being dissipated |163j . 

5.6.1 Anomalous MHD equations 

The evolution of hypermagnetic fields can be obtained from the anomalous 
magnetohydrodynamical (AMHD) equations. The AMHD equations gener- 
alize the treatment of plasma effects involving hypermagnetic fields to the 
case of finite fermionic density j^^. 

There are essential differences between the interactions of magnetic fields 
and the ones of hypermagnetic fields with matter. The ordinary electromag- 
netic field has a vector-like coupling to the fermions, while the coupling of the 
hypercharge fields is chiral. Thus, if hyperelectric {Ey) and hypermagnetic 
{Hy) fields are present simultaneously, they cause a variation of the fermionic 

/2 — * — * 

number according to the anomaly equation, ~ ji^HyEy (here g' the hy- 
percharge gauge coupling constant). Now, the presence of non-homogeneous 
hypermagnetic fields in the EW plasma with (hyper) conductivity ac always 
implies the existence of a related electric field, Ey ~ x Hy. Since for 

a general stochastic magnetic background {{Hy ■ V x HyY) 7^ 0, the non- 
uniform hypermagnetic field may absorb or release fermions and produce, 
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ultimately, baryon and lepton density perturbations because of the anomaly 
equation. 

The behaviour of cold fermionic matter with non-zero anomalous Abelian 
charges was considered in |166| where it was pointed out that the anoma- 
lous fermionic matter is unstable against the creation of Abelian gauge field 
with non-zero Chern-Simons number, which eats up fermions because of the 
anomaly. The right electron number density may be converted to the hyper- 
charge field because of a similar effect. Also the opposite effect is possible: 
hypercharge fields may be converted into fermions in a hot environment. 

The electroweak plasma in complete thermal equilibrium at a temperature 
T can be characterized hj Uf chemical potentials /ij, i = 1, ...,nf correspond- 
ing to the exactly conserved global charges 

Ni = Li-— (5.105) 

(Lj is the lepton number of the i-th generation, B is the baryon number, 
and rif is the number of fermionic generations). One should also introduce 
a chemical potential /ly corresponding to weak hypercharge, but its value 
is fixed from the requirement of the hypercharge neutrality of the plasma, 
{Y) = 0. 

It is interesting to study this plasma slightly out of thermal equilibrium, 
for instance in the situation where a non-uniform distribution of the hyper- 
magnetic field is present. Because of the anomaly, this field is coupled to the 
fermionic number densities. In principle, different chemical potentials can 
be assigned to all the fermionic degrees of freedom of the electroweak theory 
(45 if Uf = 3) and the coupled system of Boltzmann-type equations for these 
chemical potentials and the hypercharge fields may be written. Since we are 
interested in the slow processes in the plasma, this is not necessary. If the 
coupling, corresponding to some slow process, is switched off, then the elec- 
troweak theory acquires an extra conserved charge and a further chemical 
potential should be added to the system given in Eq. ()5.105|) . 

An interesting observation (which turns out to be quite important in our 
context) has been made in |lti7| Ilti9j . where it was noticed that pertur- 
bative reactions with right-handed electron chirality flip are out of thermal 
equilibrium at temperatures higher than some temperature Tr.^^ Thus, the 

^■^This temperature depends on the particle physics model, see also the discussion re- 
ported in Section 5. In the MSM Tb. ~ 80 TeV [M| HEHl HEHl ■ 
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number of right electrons is perturbatively conserved at temperatures T > 
and the chemical potential /i/j can be introduced for it. On the other hand, 
this charge is not conserved because of the Abehan anomaly, 

d,fR = -^3^,.3>'^^ (5.106) 

and it is therefore coupled to the hypermagnetic field. Here y and y are, 
respectively, the f/y(l) hypercharge field strengths and their duals, g' is the 
associated gauge coupling and yR = —2 is the hypercharge of the right elec- 
tron. 

Now we are ready to derive the anomalous MHD equations in fiat space 
|lti3| Ilti4] . The effective Lagrangian density describing the dynamics of the 
gauge fields at finite fermionic density is |170j : 

CY,e^ = -\^Y^pY"^-^J^Y'' + lJie,,uY''Y\ = (5.107) 

{g is the determinant of the metric defined in (j5.1|) : ¥^13 = ^[0X13], Vq, is the 
covariant derivative with respect to the metric ()5.1|) [notice that in the metric 
()5.H) ViaVa] = d[aYi3]\] g' is the Abelian coupling constant). The first and the 
last terms in Eq. ()5.107j) are nothing but the curved space generalization of 
the fiat-space effective Lagrangian for the hypercharge fields at finite fermion 
density |l(iH^ Il(i4j . is the ohmic current. The equations of motion for the 
hyperelectric and hypermagnetic fields are then 

or) or] TT^ 

V ■ Tiy = 0, V ■ = 0, 

V-Jy = 0, Jy = (t{Sy + v xHy), a = aca{r]), (5.108) 

with the same notations introduced in the case of the conventional MHD 
equations. 

To Eqs. ()5.108p . the evolution equation of the right electron chemical 
potential, accounting for the anomalous and perturbative non-conservation 
of the right electron number density (n^j), must be added: 

^ = -£-2^Y ■ -Hy - nriR - ni). (5.109) 
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where T is the perturbative chirahty-changing rate, F = T^, ra^ is the equi- 
hbrium value of the right electron number density, and the term proportional 
to Ey ■ Hy is the right electron anomaly contribution. 

Finally, the relationship between the right electron number density and 
the chemical potential must be specified. This relation depends upon the 
particle content of the theory. In the case of the Minimal Standard Models 
(MSM) the evolution equation of the chemical potential becomes |l()4j 

At finite hyperconductivity (in what we would call, in a MHD context, 
resistive approximation) we have that from Eq. (|5.1U8|) the induced hyper- 
electric field is not exactly orthogonal to the hypermagnetic one and, more- 
over, an extra "fermionic" current comes in the game thanks to the fact that 
we are working at finite chemical potential. Therefore in our context the 
resistive Ohm law can be written as 

^ Jy ^ ^ 

Cy = V X riy ~ 

cr 

-{—liRaily + yx'Hy\-vx'Hy, a' = ^ . (5.111) 

a \ TT J 4:71 

In the bracket appearing in Eq. (j5.111|) we can identify two different 
contributions. One is associated with the curl of the magnetic field. We will 
call this the MHD contribution, since it appears in the same way in ordi- 
nary plasmas. The other contribution contains the chemical potential and it 
is directly proportional to the magnetic field and to the chemical potential. 
This is a typical finite density effect. In fact the extra Ohmic current simply 
describes the possibility that the energy sitting in real fermionic degrees of 
freedom can be transferred to the hypermagnetic field. It may be of some 
interest to notice the analogy between the first term of Eq. ()5.111|) and 
the typical form of the ohmic current discussed in Eq. (j4.17|) appearing in 
the context of the dynamo mechanism. In the latter case the presence of a 
current (proportional to the vorticity through the a dynamo term) was indi- 
cating that large length scales magnetic fields could grow by eating up fiuid 
vortices. By inserting Ey obtained from the generalized Ohm law ()5.111|) 
in the evolution equations (|5.1(J8|) of the hypercharge fields, we obtain the 
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generalized form of the magnetic diffusivity equation: 



Aaa' 



V X (/ifiTYy) + V X ({7 X 7?) + -V^Hy- (5.112) 



drj 



In order to be consistent with our resistive approach we neglected terms con- 
taining time derivatives of the electric field, which are sub-leading provided 
the conductivity is finite. In our considerations we will also make a fur- 
ther simplification, namely we will assume that the EW plasma is (globally) 
parity-invariant and that, therefore, no global vorticity is present. There- 
fore, since the length scale of variation of the bulk velocity field is much 
shorter than the correlation distance of the hypermagnetic field, the infrared 
modes of the hypercharge will be practically unaffected by the velocity of the 
plasma, which will be neglected when the large-scale part of the hypercharge 
is concerned. This corresponds to the usual MHD treatment of a mirror 
symmetric plasma (see, e.g. Eq. ()4.32|) - ()4.34p . when a = 0). 

Eqs. (I5.112j) and ()5.11()j) form a set of AMHD equations for the hyper- 
charge magnetic field and right electron chemical potential in the expanding 
Universe. 

5.6.2 Hypermagnetic knots 

The Abelian nature of the hypercharge field does not forbid that the hyper- 
magnetic flux lines should have a trivial topological structure. This situation 
is similar to what already encountered in the case of conventional electro- 
magnetic fields with the important difference that the evolution equations of 
hypermagnetic fields are different from the ones of ordinary magnetic fields. 
After a swift summary of the properties of hypermagnetic knots (based on 
|17H I172p . some interesting applications of these hypercharge profiles will 
be reviewed. 

In the gauge Yq = 0, V • F = 0, an example of topologically non-trivial 
configuration of the hypercharge field is the Chern-Simons wave |173t I174| 



YyiZ,t) 



Y(t) sin koz, 
Y{t) cos koz, 



0. 



(5.113) 
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This particular configuration is not homogeneous but it describes a hyper- 
magnetic knot with homogeneous hehcity and Chern-Simons number density 

HyV xHy = koH^it), 

"cs = -J^5v-.y = ^ff'W. (5,114) 

where Hy = V xY, H{t) = koY{t); g' is the f/(l)y couphng. 

It is possible to construct hypermagnetic knot configurations with finite 
energy and helicity which are localized in space and within typical distance 
scale Lg. Let us consider in fact the following configuration in spherical 
coordinates 11721 



y;(7^,( 



2 Bo cos e 
2Bo sin e 



ttL, [7^2 + 1]^ 

YJn,e) = — 5.115 

' ^ vrL, [7^2 + l]2' ^ ^ 

where TZ = r/Lg is the rescaled radius and Bq is some dimensionless am- 
plitude and n is just an integer number whose physical interpretation will 
become clear in a moment. The hypermagnetic field can be easily computed 
from the previous expression and it is 

Hq{TZ,9) = tt;—— n-nsin^, 

H- ^-Tu m 8fio 7^sing fKTTfi\ 
i7,(7e,^) = -^^^^. (5.116) 

The poloidal and toroidal components of H can be usefully expressed as 
Hp = HrCr + HeCe and Tit = 'H^e^. The Chern-Simons number is finite and 
it is given by 



327r2yy" 327r2yo [7^2 + 1]^ 327r2 
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The total helicity of the configuration can also be computed 



256 Bl 
7rL2 



n 



j 

Jo 




n 



(5.118) 



We can compute also the total energy of the field 



E = -f d^x\n 




{n" + 1). 



(5.119) 



2Jv 



2 L 



s 



and we discover that it is proportional to n^. This means that one way of 
increasing the total energy of the field is to increase the number of knots and 
twists in the fiux lines. 

This type of configurations can be also obtained by projecting a non- 
Abelian SU(2) (vacuum) gauge field on a fixed electromagnetic direction |176j 

The resulting profile of the knot depends upon an arbitrary function of 
the radial distance. 

These configurations have been also studied in jl78| 1179] . In particular, 
in |179j , the relaxation of HK has been investigated with a technique different 
from the one employed in |17H I172j but with similar results. The problem 
of scattering of fermions in the background of hypermagnetic fields has been 
also studied in [MHITKT] . 

Hypermagnetic knots with large correlation scale can be also generated 
dynamically provided an unknotted hypermagnetic background is already 
present. 

Let us assume that dynamical pseudoscalar particles are evolving in the 
background geometry given by Eq. (j5.ip . The pseudoscalars are not a source 
of the background (i.e. they do not affect the time evolution of the scale 
factor) but, nonetheless, they evolve according to their specific dynamics 
and can excite other degrees of freedom. 

The action describing the interaction of a dynamical pseudoscalar with 
hypercharge fields can be written as 



^ In order to interpret these solutions it is very interesting to make use of the Clebsh de- 
composition. The inipHcations of this decomposition (beyond the hydrodynamical context, 
where it was originally discovered) have been recently discussed (see |177j and references 
therein) . 



S = j d'xy/^^^g-^d^^dp^ - - -Y^pY^P + c^y„^y"/5] . (5.120) 
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This action is quite generic. In the case Vi^ip) = {rn'^/'2)ip^ Eq. (j5.12U|) is 
nothing but the curved space generahzation of the model usually employed 
in direct searches of axionic particles |1821 118JI I184j . The constant in front 
of the anomaly is a model-dependent factor. For example, in the case of 
axionic particles , for large temperatures T > mw, the Abelian gauge fields 
present in Eq. (j5.12(jp will be hypercharge fields and c = c^yo^' / (27r) where 
a' = g''^/4:7c and c^y is a numerical factor of order 1 which can be computed 
(in a specific axion scenario) by knowing the Peccei-Quinn charges of all the 
fermions present in the model 1186j . For small temperatures T < mw 
we have that the Abelian fields present in the action ()5.120j) will coincide 
with ordinary electromagnetic fields and c = c^^aem/27r where ttem is the 
fine structure constant and c^^ is again a numerical factor. 

The coupled system of equations describing the evolution of the pseu- 
doscalars and of the Abelian gauge fields can be easily derived by varying 
the action with respect to ijj and Y^, 



'-9 



--d„ 



+ 



dV 



di/j AM 



Y nY"'^ 
J- afi-'- 1 



VuY^" = — V.V^r 



0, 



(5.121) 



where, 



--d„ 



-9Y 



-.d„ 



(5.122) 



are the usual covariant derivatives defined from the background FRW metric 
Eqs. ()5.121|) can be written in terms of the physical gauge fields 

r + 2W - V'ij + a'^ = ■ Hy, 

dip a'' M 



V ■nY = o, V xSy + n^ 

c 



y-0. 



V-£y = j^ViJ ■ Hy, 



^ xnY = s'- 



M 



iP'Hy + X 



(5.123) 



We want now to study the amplification of gauge field fluctuations in- 
duced by the time evolution of if). Then, the evolution equation for the 
hypermagnetic fluctuations Tiy can be obtained by linearizing Eqs. ()5.123|) . 
We will assume that any background gauge field is absent. In the linearisa- 
tion procedure we will also assume that the pseudoscalar field can be treated 



76 



as completely homogeneous (i.e.|V'?/'| <^ 4'')- This seems to be natural if, 
prior to the radiation dominated epoch, an inflationary phase diluted the 
gradients of the pseudoscalar. 

In this approximation, the result of the linearization can be simply written 
in terms of the vector potentials in the gauge = and V ■ V = 0: 

y" - V^y + ^^'V X 3^ = 0, (5.124) 
dV 

ij + 3Hij + — = 0, (5.125) 
dip 

By combining the evolution equations for the gauge fields we can find a 
decoupled evolution equation for Hy, 

H'y -V^Hy + xHy = 0. (5.126) 

From this equation is already apparent that the pseudo-scalar vertex 
induces an interaction in the two physical polarizations of the hypermag- 
netic field. Giving initial conditions which are such that Tiy ^ with 
Hy ■ V X Hy = a profile with Hy ■ V x TCy can be generated provided 



6 Two approaches to the Origin 

Since the flux is conserved the ratio between the magnetic energy density, 
PBiL,rj) and the energy density sitting in radiation, p^{ri) is almost constant 
and therefore, in terms of this quantity (which is only scale dependent but not 
time dependent), the most simplistic dynamo requirement can be rephrased 
(see Section 3) as 

ML) = > 10-'\ L-lMpc, (6.1) 

to be compared with the value re ~ 10~^ which would lead to the galactic 
magnetic field only thanks to the collapse and without the need of dynamo 
action (this would be the case when the magnetic field is fully primordial). 
Notice that Eq. ()6.1|) assumes that the magnetic fiux is exactly frozen into 
the plasma element. 
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As previously pointed out during the discussion of the dynamo mecha- 
nism, the requirement expressed by Eq. is unreahstic: it would corre- 
spond to thirty e-folds of amplification during galactic rotation and perfect 
flux freezing during the collapse of the protogalaxy. In a more realistic sit- 
uation situation, taking into account the effectively achievable amplification 
of the dynamo action the requirement would rather be, prior to collapse, 

= ^^i^ > 10-23, L~lMpc. (6.2) 

In spite of the richness of the theoretical models, the mechanisms for magnetic 
field generation can be divided, broadly speaking, into two categories: as- 
trophysical and cosmological. The cosmological mechanisms can be divided, 
in turn, into causal mechanisms (where the magnetic seeds are produced at 
a given time inside the horizon) and inflationary mechanisms where correla- 
tions in the magnetic field are produced outside the horizon. 



6.1 Inhomogeneous MHD equations 

The first attempts in this direction have been made by Biermann |187j and 
Harrison |189| 1190^ . The Biermann battery mechanism is easy to un- 
derstand from the form of the generalized Ohm law discussed in Eq. ()4.17|) . 
MHD equations are linear and homogeneous in the magnetic field intensities. 
The idea is then to look for a natural source term which seemed to be pro- 
vided by the thermoelectric current already introduced in the context of the 
generalized Ohm law . Consider, indeed, Eq. (j4.17|) in the approximation 
where the Hall term is neglected (the argument can be also generalized to 
the case of non vanishing Hall term). Then, the Ohmic electric field will not 
be the one simply obtained in resistive MHD but it will get a contribution 
from the thermoelectric term: 

- V X 1 - , , 

E = -vx B + VPe- (6.3) 

Ana erie 

Using the other MHD equations in the incompressible approximation the new 
form of the magnetic diffusivity equation can be derived: 

r9R 1 1 

— - V X (tJ X 5) - -—V^B = -Vne x VPe- (6.4) 

ot Ana en| 
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The magnetic diffusivity equation has now a source terms which does not 
vanish provided the gradients of the charge density of the electrons is not 
parallel to the pressure gradients. As a consequence of the presence of the 
source term, the magnetic field will grow linearly in time until the thermo- 
electric source term is comparable with the dynamo term. The estimate of 
the magnetic field intensity prior to the onset of galactic rotation will then 
be 

|5| ~ ^|Vne X VPel, (6.5) 

leading to magnetic fields 0(10^^'^) G over typical scales of the order of 10 
kpc (see also |188j for a recent analysis in a similar framework). 

In a complementary perspective, Harrison |189| 11901 1191] discussed the 
possibility that the pressure gradients are strictly vanishing, i.e. Vpe = 
0. This is a reasonable assumption, for instance, in a radiation-dominated 
stage of expansion where pressure gradients are expected to be negligible by 
the global homogeneity of the background. In this situation, an evolution 
equation for the vorticity can be derived 

^{a'u + —B) = (6.6) 
or] nip Airamp 

where u = V x v and rup is the ion mass. If we now postulate that some 
vorticity was present prior to decoupling, then Eq. ()6.6|) can be solved and 
the magnetic field can be related to the initial vorticity as 

mp ^ /ai^"^ 



B . (6.7) 

e V a / 

It is clear that if the estimate of the vorticity is made prior to equality (as 
originally suggested by Harrison |1 89 j ) of after decoupling as also suggested, a 
bit later[192] can change even by two orders of magnitude. Prior to equality 
\Lo{t) ~ 0.1/t and, therefore, |-Beq| ~ lO^^^G. If a similar estimate is made 
after decouphng the typical value of the generated magnetic field is of the 
order of 10~^* G. So, in this context, the problem of the origin of magnetic 
fields is circumvented by postulating an appropriate form of vorticity whose 
origin must be explained. 

The idea was employed, later on, in the context of the physics of topo- 
logical defects. Vachaspati and Vilenkin |193j have suggested that cosmic 
strings with small scale structure may be a source of the wanted vorticity. 
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The argument is, in short, that since matter flow in baryonic wakes is turbu- 
lent, velocity gradients will be induced in the flow by the small-scale wiggles 
which produced the required vorticity. Furthermore, dynamical friction be- 
tween cosmic strings and matter may provide a further source of vorticity 
|194j . There have been also studies trying to generate large-scale magnetic 
flelds in the context of superconducting cosmic strings (see, for instance, |195j 
and references therein). 

Recently the possible generation of large-scale magnetic flelds prior to hy- 
drogen recombination has been discussed in |196tll97] . The vorticity required 
in order to produce the magnetic flelds is generated, according to [196 , by 
the photon diffusion at second order in the temperature fluctuations. In a 
similar perspective Hogan |198j got much less optimistic estimates which, ac- 
cording to |196tll97j . should be attributed to different approximation schemes 
employed in the analysis. 

6.2 Inflationary mechanisms 

Different astrophysical objects, of different physical sizes, have comparable 
magnetic flelds. This coincidence is hard to explain in the context of causal 
mechanisms of generation. The analogy with structure formation, already 
presented in the Introduction, is here useful. In the late seventies, prior to 
the formulation of the inflationary paradigm, the initial conditions for the 
density contrast were taken as primordial input. Later on, in the context 
of inflationary models, the primordial spectrum of curvature and density 
fluctuations could be calculated. 

During inflation, in fact, flelds of various spin are present and they can 
be excited by the dynamical evolution of the geometry. In the context of 
inflationary models of particular relevance are the scalar and tensor fluc- 
tuations of the geometry, corresponding respectively, to fluctuations of the 
scalar curvature and to gravitational waves. Gravitational waves and curva- 
ture perturbations obey evolution equations which are not invariant under 
Weyl rescaling of the four-dimensional metric. Then the quantum mechan- 
ical fluctuations in the corresponding flelds will simply be amplifled by the 
evolution of the background geometry. 

It is interesting to speculate that large scale magnetic flelds could be 
produced thanks to a similar mechanism. The major obstruction to this 
type of models is that gauge flelds are not amplifled thanks to the evolution 
of the background geometry the reason being that the evolution equations of 



80 



gauge fields are invariant under Weyl rescaling of tlie metric. 

Wlien conformal invariance is broken, by some means, one is often led 
to estimate the amplitude of gauge field fluctuations arising as a result of 
the breaking. Clearly the detailed amount of amplified gauge fields will be 
specific of the particular model. In the following various ways of breaking 
conformal invariance will be listed without getting through the details of the 
calculations. Before doing so it is anyway instructive to introduce some gen- 
eral considerations stressing the analogy between the production of magnetic 
inhomogeneities and the production of gravitational inhomogeneities whose 
late time evolution leads the anisotropies in the CMB. 

The tensor modes of the geometry are described by a rank-two (transverse 
and traceless) tensor in three spatial dimensions i.e. 



ds^ = a^{vi)[drf — {6ij + hij)dx'^dx^], h\ 
obeying the equation, in Fourier space. 



0, 



(6.8) 



0, 



(6.9) 



where Hk = ah^ and hj^ is the Fourier mode of each polarization. In this 
equation the "pump field" is simply given by the scale factor. When k"^ 3> 
\a"/a\ the mode is said to be adiabatically damped: in fact, in this regime, 
\hk\ — a~^, i.e. decreasing in an expanding Universe. In the opposite regime, 
i.e. -C \a"/a\ the mode is super-adiabatically amplified. In fact, during 
a de Sitter or quasi-de Sitter fi^ ~ '^('7) and hk is constant. Hence, for the 
whole time the given mode spends under the "potential barrier" of Eq. (j6.9p . 
it is amplified. 

A similar equation holds for the canonical normal mode for the scalar 
fluctuations of the geometry, i.e. 



Vk + 



Vk = 0. 



(6.10) 



The variable Vk is defined as 



V - ""^ 

K-k — ■ 

Z 



(6.11) 



^^The theory of cosmological perturbations is assumed in the fohowing considerations. 
Comprehensive accounts of the relevant topics can be found in [199. .200, (see also 1201) 1. 
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where 5(pk is the fluctuation of the inflaton, TZk is the curvature perturbation 
and if) is the scalar fluctuation of the geometry in the conformally Newtonian 
gauge |199t 1200] . If the inflaton has an exponential potential zij]) ~ a{rj). 

In the case of gauge fields, each of the two polarization of the appropri- 
ately rescaled vector potential obeys, the following equation 

A = 0. (6.12) 

In specific models, g{ri) may be associated with the gauge coupling. However 
g~^{ri) can also be viewed as a generic pump field arising as a result of the 
breaking of conformal invariance. The time dependence of the potential term 
is rather common to different models: it goes to zero for |?7| oo as {u"^ — 
l/4)r]~'^. The numerical coefficient appearing in the potential determines 
the strength and spectrum of the amplified gauge fields whose subsequent 
evolution has to be however computed at finite conductivity. In Fig. IHlthe 



(g")"g ' 




Figure 6: The effective potential appearing in Eq. ()6.12|) is illustrated in 
general terms. In this example the pump field leading to amplification of 
electromagnetic quantum fluctuations is assumed to be constant prior to the 
onset of BBN. On the vertical axis few relevant physical frequencies (to be 
compared with the height of the potential) have been reported. 

typical form of the potential term appearing in Eq. ()6.12|) is illustrated. 
Different comoving frequencies go under the barrier at different times and 
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the amount of amplification is roughly proportional to the time spent under 
the barrier. Clearly, given the generic form of the barrier, smaller frequencies 
are more amplified than the frequencies comparable with the height of the 
barrier at rji. In Fig. IHl the explicit numerical value of the height of the 
barrier corresponds to a (present) frequency of 10^ Hz which can be realized 
if the pump field goes to a constant right after the end of a conventional 
inflationary phase followed by a radiation-dominated stage of expansion. 

In the case of gravitational waves, and, with some caveats, in the case of 
scalar metric fluctuations, the "potential barriers" appearing, respectively, in 
Eqs. fj6.9p and ()6.10|) may be related with the inverse Hubble radius. Hence, 
in the case of metric fluctuations, a mode which is under the barrier is also, 
with a swift terminology, outside the horizon (see Fig. [7j) . This is the reason 
why, following the conventional nomenclature, in Fig. IHlthe moment when a 
given scale gets under the barrier has been denoted by r/ex (i-e. horizon exit). 
According to the same convention, the moment when a given scale gets out 
the potential barrier is labeled by 7]^^ (i-e. horizon re-entry). In Fig. IHlfew 




10 20 30 40 50 60 70 



Figure 7: The evolution of a given physical wavelength is illustrated in the 
case of an inflationary model with minimal duration. The region above the 
full curve (denoting the Hubble radius rn) = H^^) corresponds, for the tensor 
and scalar modes of the geometry, to the region where a given comoving 
wavenumber k is below the "potential barrier" appearing in Eqs. (j6.9|) and 

fiFnnii . 

relevant frequencies have been compared with the height of the barrier. Con- 
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sider, for instance, ujq i.e. the typical scale of gravitational collapse which is 
of the order of 1 Mpc, i.e. 10"^'* Hz. In Fig. IHlthe physical frequency has been 
directly reported. Another interesting frequency is Ua ~ mHz corresponding 
to the present value of the magnetic diffusivity momentum. The amplifica- 
tion caused by the parametric amplification of the vacuum fluctuations can 
be computed by solving Eq. (j6.12j) in the different regimes. 

' ^ ^ Vexy 

H g^iv')] , Vex <r] < Vre 

Ak = -^[c+(A;)e-^'='' + c_(A;)e*^^] , r, > r^^e, (6.13) 

where ( C^, D^, and c±{k) are integrations constants). 

The mixing coefficients c±{k), determining the parametric amplification 
of a mode k'^ < \V{rii)\, computed by matching these various branches of the 
solution reported in Eq. ()6.13|) . One finds: 

2,fce-^^(''-^^-)c±(A;) = ^^(-^ ^ ik) ± %(_^ + ^k) 

Qre Qre Qex Qex 

± ^_(_^ + ,A;)(-^T2fc) r'(?2^r/. (6.14) 

9ex9re 9 ex 9re J Vex 

Similar calculations can be performed in order to obtain the spectrum of 
scalar and tensor fluctuations of the geometry |204| I2U5| |206j . 

Suppose now to make a simple estimate. Assume, for instance, that g is 
evolving prior to rji according to the dynamics dictated by a given model. 
Suppose also that after 771 the Universe is suddenly dominated by radiation, 
and 9' ^ for rj > rji. In this situation all the modes reenter during radiation 
and the amplification will be roughly given, to leading order by 

|c_| ~ — . (6.15) 

9ex 

If the function g{ri) is identified with the evolving gauge coupling this result 
suggest that in order to have large amplification, g{7]) has to grow from 
smaller to larger values. This is what happens, for instance, in the case of 
pre-big bang models where g ~ e"^/^ and ip is the for-dimensional dilaton field 
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Since we ought to estimate the amphfication of an initial quantum me- 
chanical fluctuations, a fully quantum mechanical treatment is certainly ap- 
propriate also in view of the discussion of the correlation properties of the 
obtained fluctuations. This analysis has been performed in |2U7j where the 
squeezing properties of the obtained photons have also been discussed. 

The perturbed effective Lagrangian density 



l2 



L^T]) = d xC{x, 1]), 



9 Vfl' / 

(6.16) 

describes the evolution of the two ( a = (8>, ©) transverse degrees of freedom 
defined by the Coulomb gauge condition Aq = and V ■ A = (the prime 
denotes differentiation with respect to conformal time). The fields Aa = gAa 
have kinetic terms with canonical normalization and the time evolution given 
in Eq. ()6.12|) stems from the Euler-Lagrange equations derived from Eq. 
()6.16|1 . By functionally deriving the the action the canonically conjugated 



momenta can be obtained leading to the Hamiltonian density and to the 
associated Hamiltonian 



n{x,ri) = Iy1 



T^l + {diA, 



2'-A^n^ 



H{r]) = J d^xHix). 



(6.17) 



The operators corresponding to the classical polarizations appearing in the 
Hamiltonian density 



Aa{k,7])e'' 



ik-x 




Aa{x,7]) 
Aa{k,T]) 
7Taix,r]) -- 



(6.18) 



obey canonical commutation relations and the associated creation and anni- 
hilation operators satisfy [dk,a, d^j^ 

The (two-modes) Hamiltonian contains a free part and the effect of the 
variation of the coupling constant is encoded in the (Hermitian) interaction 
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term which is quadratic in the creation and annihilation operators whose 
evolution equations, read, in the Heisenberg picture 

ddk,a _ ^ g' ^] 
dv 9 ' 

d^i a g' . /« 1n^ 

—7^ = ikal^ a^k,a- (6.19) 

dr] ' g 

The general solution of the previous system of equations can be written in 
terms of a Bogoliubov-Valatin transformation, 

dkAv) = fJ'kAv)bk,a + ^kAv)b-k,a 

dk,aiv) = fJ'l,ai^)iia + K,ai^)b-k,a (6.20) 

where dk,a{0) = and d^k,a{^) ~ ^-k,a- Notice that the Bogoliubov co- 
efficients are the quantum analog of the mixing coefficients discussed in the 
semiclassical approach to the problem. 

Unitarity requires that the two complex functions ^k{v) and I'kiv) are 
subjected to the condition \fik{v)\^ ~ Wk{v)\^ = 1 which also implies that 
IJ^kiv) and I'kiv) can be parameterized in terms of one real amplitude and 
two real phases 



l^k 



e*^* cosh rfc, = e'^^'^"'^"^ sinh r^, (6.21) 



(r is sometimes called squeezing parameter and (pk is the squeezing phase; 
from now on we will drop the subscript labeling each polarization if not 
strictly necessary). The total number of produced photons 

{0-kOk\di{r])dki'n) + alfca_fe|OfcO_fc) = 2 n^. (6.22) 

is expressed in terms of = sinh^r^, i.e. the mean number of produced 
photon pairs in the mode k. Inserting Eqs. ()6.20|) - ()6.22|) into Eqs. ()6.19p we 
can derive a closed system involving only the fik and the related phases: 

-^ = -2/nfe^cos20fc, 6.23 
drj g 

— = -A; + — — — sm20fc, (6.24) 
dv 9 j[nk) 

—k-\ — — smzcpk, (6.25) 



dr] g duk 
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where /(rifc) = ^Ukiuk + 1). 

In quantum optics [208, 20^ the coherence properties of hght fields have 
been a subject of intensive investigations for nearly half a century. In the 
present context the multiparticle states described so fare are nothing but 
squeezed states of the electromagnetic field [208^ l2U9j . In fact, up to now 
the Heisenberg description has been adopted. In the Schrodinger picture the 
quantum mechanical states obtained as a result of the time evolution are 
exactly squeezed vacuum states |2U8| l2U9j . 

Magnetic fields over galactic scales have typical frequency of the order 
]^g-i4_]^g-i5 which clcarly fall well outside the optical range. Thus, the 
analogy with quantum optics is only technical. The same quantum opti- 
cal analogy has been successfully exploited in particle |21()j and heavy-ions 
physics |211j of pion correlations (in order to measure the size of the strongly 
interacting region) and in the phenomenological analysis of hadronic multi- 
plicity distributions. 

The interference between the amplitudes of the magnetic fields (Young 
interferometry |212j . in a quantum optical language) estimates the first order 
coherence of the magnetic background at different spatial locations making 
use of the two-point correlation function whose trace over the physical po- 
larizations and for coincidental spatial locations is related to the magnetic 
energy density. Eqs. ()6.23|) - ()6.25|) can be solved once the pump field is 
specified but general expressions can be also obtained |2U7j . 

6.2.1 Conventional inflationary models 

In conventional infiationary models it is very difficult to produce large scale 
magnetic fields with phenomenologically relevant strength. This potential 
difficulty has been scrutinized in various investigations [ 2U21 12131 12141 I^T3] . 

Turner and Widrow |213j listed a series of field theory models in de Sitter 
space with the purpose of finding a natural way of breaking conformal in- 
variance. The first suggestion was that conformal invariance may be broken, 
at an effective level, through the coupling of photons to the geometry |21tij . 
Typically, the breaking of conformal invariance occurs through products of 
gauge-field strengths and curvature tensors, i.e. 

^F^.F^^i?'^^"^ l-R^.F^'PF'^g^p, l-F^pF'^^R (6.26) 

rrt'^ 

where m is the appropriate mass scale; R^uap and R^^ are the Riemann and 
Ricci tensors and R is the Ricci scalar. If the evolution of gauge fields is 
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studied during a phase of de Sitter (or quasi-de Sittter) expansion, then the 
amphfication of the vacuum fluctuations induced by the couphngs hsted in 
Eq. ()6.26|) is minute. The price in order to get large amplification should be, 
according to ^213j, an explicit breaking of gauge-invariance by direct coupling 
of the vector potential to the Ricci tensor or to the Ricci scalar, i.e. 

RA^A^, R^uA^Ar (6.27) 

In |213j two other different models were proposed (but not scrutinized in 
detail) namely scalar electrodynamics and the axionic coupling to the Abelian 
field strength. 

Dolgov [215J considered the possible breaking of conformal invar iance due 
to the trace anomaly. The idea is that the conformal invariance of gauge fields 
is broken by the triangle diagram where two photons in the external lines 
couple to the graviton through a loop of fermions. The local contribution 
to the effective action leads to the vertex {y/^Y^Ta/sF'^'^ where e is a 
numerical coefficient depending upon the number of scalars and fermions 
present in the theory. The evolution equation for the gauge fields, can be 
written, in Fourier space, as 

A'l + ^HA'k + k'^Ak = 0, (6.28) 
8 

and it can be shown that only if e > the gauge fields are amplified. Fur- 
thermore, only is e ~ 8 substantial amplification of gauge fields is possible. 

In a series of papers |217| I218| I219j the possible effect of the axionic 
coupling to the amplification of gauge fields has been investigated. The idea 
is here that conformal invariance is broken through the explicit coupling of 
a pseudo-scalar field to the gauge field (see Section 5), i.e. 

V^C^7«em^i^a/3^"^ (6.29) 

where F"^ is the dual field strength and where numerical factor of 

order one. Consider now the case of a standard pseudoscalar potential, for 
instance m'^ip'^, evolving in a de Sitter (or quasi-de Sitter space-time). It 
can be shown, rather generically, that the vertex given in Eq. (j6.29|) leads 
to negligible amplification at large length-scales. The coupled system of 
evolution equations to be solved in order to get the amplified field is similar 
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to Eqs. (j5.123|) already introduced in the duscussion of hypermagnetic fields 



B" - V^B - ^^V^'V xB = 0, (6.30) 
^p" + 2n^' + m^c?i) = 0. (6.31) 

From Eq. ()6.30p . there is a maximally amplified physical frequency 

2nM 27r ^ ^ 

where the second equality follows from i/j ^ a cosmt (i.e. tpmax ~ mM). 

The amplification for uj ~ ujm&x is of the order of exp [maem/ (2vrif)] where 
H is the Hubble parameter during the de Sitter phase of expansion. From 
the above expressions one can argue that the modes which are substantially 
amplifed are the ones for which ujrna.x ^ H. The modes interesting for the 
large-scale magnetic fields are the ones which are in the opposite range, i.e. 
i^max -C H. Clearly, by lowering the curvature scale of the problem the 
produced seeds may be larger and the conclusions much less pessimistic |i219j . 

Another interesting idea pointed out by Ratra |214j is that the electro- 
magnetic field may be directly coupled to the inflaton field. In this case the 
coupling is specified through a parameter a, i.e. e°"^ Fa^F"'^ where (p is the 
inflaton field in Planck units. In order to get sizable large-scale magnetic 
fields the effective gauge coupling must be larger than one during inflation 
(recall that ip is large, in Planck units, at the onset of inflation). 

In |22Uj it has been suggested that the evolution of the Abelian gauge cou- 
pling during inflation induce the growth of the two-point function of magnetic 
inhomogeneities. This model is different from the one previously discussed 
|214j . Here the dynamics of the gauge coupling is not related to the dynamics 
of the inflaton which is not coupled to the Abelian field strength. In partic- 
ular, rs(Mpc) can be as large as 10~^^. In |22Uj the MHD equations have 
been generalized to the case of evolving gauge coupling. Recently a scenario 
similar to |22Uj has been discussed in |221j . 

In the perspective of generating large scale magnetic fields Gasperini 
|222j suggested to consider the possible mixing between the photon and the 
graviphoton field appearing in supergravity theories (see also, in a related 
context |223j ). The graviphoton is the massive vector component of the 
gravitational supermultiplet and its interaction with the photon is specified 
by an interaction term of the type XF^^G^^ where G^i, is the filed strength 



89 



of the massive vector. Large-scale magnetic fields with rs(Mpc) > 10~^^ can 
be obtained if A ~ 0{1) and for a mass of the vector m ~ lO^TeV. 

Bertolami and Mota |224j argue that if Lorentz invariance is sponta- 
neously broken, then photons acquire naturally a coupling to the geometry 
which is not gauge-invariant and which is similar to the coupling considered 
in inSl- 

Finally Davis and Dimopoulos |225j considered the possibility of phase 
transitions taking place during infiation. They found that sizable large-scale 
magnetic fields can be generated provided the phase transition occurs in the 
last 5 e-foldings of the infiationary stage of expansion. 

6.2.2 Abelian Higgs model 

While the coupling of electromagnetic field to the metric and to the charged 
fields is conformally invariant, the coupling of the charged scalar field to 
gravity is not. Thus, vacuum fiuctuations of the charged scalar field can 
be amplified during infiation at super-horizon scales, leading potentially to 
non-trivial correlations of the electric currents and charges over cosmological 
distances. The fiuctuations of electric currents, in turn, may induce magnetic 
fields through Maxwell equations at the corresponding scales. The role of the 
charged scalar field may be played by the Higgs boson which couples to the 
hypercharge field above the electroweak phase transition; the generated hy- 
percharged field is converted into ordinary magnetic field at the temperatures 
of the order of electroweak scale. 

No detailed computations were carried out in |213j in order to substanti- 
ate this idea. The suggestion of [2l3j was further developed quite recently in 
|226j for the standard electroweak theory with an optimistic conclusion that 
large scale magnetic fields can be indeed generated. In ^227] a supersymmet- 
ric model was studied. In |134j , previous treatments have been further scru- 
tinized by computing, with higher accuracy, the amplification of the charged 
scalar field and the damping induced by the conductivity. It turns out that 
the resulting magnetic fields are insufficient in order to provide reasonable 
seeds for the dynamo amplification. 

Introducing appropriately rescaled fields the action of the Abelian-Higgs 
model in a conformally flat FRW space-time can be written as 



/ 



a' 



S 



(6.33) 



a 



90 



Now, since the evolution equation of the charged scalar is not conformally 
invariant, current density and charge density fluctuations will be induced. 
Then, by employing a Vlasov-Landau description similar to the one in- 
troduced in Section 5, the resulting magnetic field will be of the order of 
Bdec ~ 10~''°T|gj, which is, by far, too small. Later it has been proposed that 
much larger magnetic fields may be obtained in the context of the Abelian- 
Higgs model |228j (see however |229j for a detailed criticism of this proposal). 

6.2.3 Internal dimensions 

If internal dimensions are dynamical, then Weyl invariance may be natu- 
rally broken |23Uj . Consider a pure electromagnetic fiuctuation decoupled 
from the sources, representing an electromagnetic wave propagating in the 
c?-dimensional external space such that = A^{x, rj), Aa = 0. In the metric 
given in Eq. (j5.1|) the evolution equation of the gauge field fiuctuations can 
be written as 



where Fajs = V[qA^] is the gauge field strength and G is the determinant of 
the D dimensional metric. Notice that if n = the space-time is isotropic 
and, therefore, the Maxwell's equations can be reduced (by trivial rescaling) 
to the fiat space equations. If n 7^ we have that the evolution equation of 
the electromagnetic fiuctuations propagating in the external (i-dimensional 
manifold will receive a contribution from the internal dimensions which can- 
not be rescaled away |233j In the radiation gauge {Aq = and VjA* = 0) 
the evolution the vector potentials can be written as 



The vector potentials Ai are already rescaled with respect to the (conformally 
fiat) d + 1 dimensional metric. In terms of the canonical normal modes of 
oscillations Ai = b^^'^Ai the previous equation can be written in a simpler 

^^Notice that the electromagnetic field couples only to the internal dimensions through 
the determinant of the Z3-dimensional metric. In string theories, quite generically, the 
one-form fields are also coupled to the dilaton field. This case has been already analyzed 
in the context of string inspired cosmological scenarios and will be discussed later. 

^^For a discussion of gauges in curved spaces see |232| . 




(6.34) 




(6.35) 
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form, namely 

2 

— * 71 71 

A': - V{v)A - VU. = 0, V{v) = -J"' + -T'. (6.36) 

In order to estimate the amplification of the gauge fields induced by the 
evolution of the internal geometry we shall consider the background metric of 
Eq. ()5.1Hj) in the case of maximally symmetric subspaces = ^ah = Sab- 

Suppose now that the background geometry evolves along three different 
epochs. During the first phase (taking place for ] — oo, — t^i]) the evolution is 
truly multidimensional. At r] = —7]i the multidimensional dynamics is con- 
tinuously matched to a radiation dominated phase turning, after decoupling, 
into a matter dominated regime of expansion. During the radiation and 
matter dominated stages the internal dimensions are fixed to their (present) 
constant size in order not to conflict with possible bounds arising both at 
the BBN time and during the matter-dominated epoch. The evolution of the 
external dimensions does not affect the amplification of the gauge fields as it 
can be argued from Eq. ()6.36|) :in the limit n — (i.e. conformally invariant 
background) Eq. (|6.36|) reduces to the fiat space equation. A background 
with the features described above has been introduced in Eq. (|5.14j) . 

Defining, respectively, Bbbn and 6o as the size of the internal dimensions 
at the BBN time and at the present epoch, the maximal variation allowed to 
the internal scale factor from the BBN time can be expressed as bsBN/bo ~ 
1 + e where |e| < 10^^ jll61 1117| I118j . The bounds on the variation of the 
internal dimensions during the matter dominated epoch are even stronger. 
Denoting with an over-dot the derivation with respect to the cosmic time 
coordinate, we have that \b/b\ < 10~^Hq where Hq is the present value of 
the Hubble parameter |116j . The fact that the time evolution of internal 
dimensions is so tightly constrained for temperatures lower of 1 MeV does 
not forbid that they could have been dynamical prior to that epoch j231j . 

In the parameterization of Eq. ()5.14|) the internal dimensions grow (in 
conformal time) for A < and they shrink for A > 0. 

^®To assume that the internal dimensions are constant during the radiation and matter 
dominated epoch is not strictly necessary. If the internal dimensions have a time variation 
during the radiation phase we must anyway impose the BBN bounds on their variation 
[1161 11171 1118| . The tiny variation allowed by BBN implies that b{ri) must be effectively 
constant for practical purposes. 
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By inserting this background into Eq. (j6.36p we obtain that for t] < —rji 



V{r,) 



nX 
4^ 



{n\ - 2) 



(6.37) 



whereas V{ri) — > for t] > —rji. Since V{ri) goes to zero for rj ±oo we 
can define, in both hmits, a Fourier expansion of Ai in terms of two distinct 
orthonormal sets of modes. The amphfication of the quantum mechanical 
fluctuations of the gauge fields can then be computed using the standard 
techniques (see, for instance, |234t I235t I236j ). In |23Uj it has been shown that 
in simple models of dimensional decoupling magnetic fields can be generated 
with strength compatible with the bound of Eq. ()6.2|1 . In particular there 
is the interesting possibility that large-scale magnetic fields are produced in 
the case when internal dimensions also expand while the external ones also 
expand but at a different rate. 

It should be mentioned that the interplay between gauge fields and large 
(or infinite) extra-dimensions is still under investigation. There is, at the 
moment, no definite model leading to the generation of large-scale magnetic 
field in a framework of infinite internal dimensions. The main reason is 
that it is difficult to build reasonable cosmological models with large extra- 
dimensions. More specifically, one would like to have a model where gauge 
fields arise naturally in the bulk. In six- dimensions interesting models can be 
constructed |239j where localized gauge zero modes may be naturally present 
(in analogy with what happens in the case of six-dimensional warped models 



6.2.4 String cosmological models 

Large scale magnetic fields may also be produced in the context of string 
cosmological models |2U2t I289j (see also |238j for the analysis of the ampli- 
fication of magnetic fields during a phase of coherent dilaton oscillations^^). 
In many respects the case of string cosmological models is related to the one 
where the gauge coupling is dynamical. However, in the string cosmological 
context, the dynamics of the gauge coupling and of the geometry are con- 
nected by the solutions of the low-energy /5-functions. The basic evolution 
of the background and of the gauge coupling has been discussed in Section 

^^The possible production of (short scale) magnetic fields by parametric resonance ex- 
plored in ,238, , has been subsequently analyzed in the standard inflationary contex | 2U3| 
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5. In order to achieve a large amplification of the quantum fluctuations of 
the gauge fields, the gauge coupling should be sufficiently small when the 
typical scale of the gravitational collapse hits the "potential barrier" . In par- 
ticular, gexiujc) < 10~^^ is the dynamo requirement has to be satisfied. It is 
difficult to produce large scale magnetic fields with reasonable amphtudes if 
the pre-big bang phase matches immediately to the post-big bang evolution 
|237j . However, phenomenologically consistent pre-big bang models lead to 
a sufficiently long stringy phase when the dilaton and the curvature scale 
are roughly constant. If this phase is included the conditions expressed by 
Eqs. ()6.H) and ()6.2|) are easily satisfied. Furthermore, there are regions in 
the parameter space of the model where r^iuJc) ~ 10~^ which implies that 
the galactic magnetic field may be fully primordial \2[)'2\ 1289] , i.e. no dynamo 
action will be required. 

6.3 Inside the Hubble radius 

The experimental evidence concerning large-scale magnetic fields suggests 
that magnetic fields should have similar strength over different length-scales. 
In this sense inflationary mechanisms seem to provide a rather natural ex- 
planation for the largeness of the correlation scale. At the same time, the 
typical amplitude of the obtained seeds is, in various models, rather minute. 

Primordial magnetic flelds can also be generated through physical mech- 
anisms operating inside the Hubble radius at a given physical time. Partic- 
ularly interesting moments in the life of the Universe are the epoch of the 
electroweak phase transition (EWPT) or the QCD phase transition where 
magnetic flelds may be generated according to different physical ideas. In 
the following the different proposals emerged so far will be reviewed. 

6.3.1 Phase transitions 

At the time of the EWPT the typical size of the Hubble radius is of the 
order of 3 cm and the temperature is roughly 100 GeV. Before getting into 
the details of the possible electroweak origin of large-scale magnetic flelds 
it is useful to present a kinematical argument based on the evolution of the 
correlation scale of the magnetic flelds |142j . 

Suppose that, thanks to some mechanism, sufficiently large magnetic 
fields compatible with the critical density bound are generated inside the 
Hubble radius at the electroweak epoch. Assuming that the typical coherence 
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length of the generated magnetic fields is maximal, the present correlation 
scale will certainly be much larger but, unfortunately, it does not seem to 
be as large as the Mpc scale at the epoch of the gravitational collapse. As 
already mentioned, the growth of the correlation scale may be enhanced, by 
various processes such as inverse cascade and helical inverse cascade. For 
instance, if the injection spectrum generated at the electroweak epoch is 
Gaussian and random a simple estimate shows that the present correlation 
scale is of the order of 100 AU which is already larger than what one would 
get only from the trivial expansion of length-scales (i.e. 1 AU) [142j. If, in a 
complementary perspective, the injection spectrum is strongly helical, then 
the typical correlation scale can even be of the order of 100 pc but still too 
small than the typical scale of the gravitational collapse. 

Large-scale magnetic fields can be generated at the electroweak epoch in 
various ways. Consider first the case when the phase transition is strongly 
first order. 

Hogan |152j originally suggested the idea that magnetic fields can be 
generated during first-order phase transitions. Since during the phase transi- 
tion there are gradients in the radiation temperature, similar thermoelectric 
source terms of MHD equations (which were discussed in the context of the 
Biermann mechanism) may arise. The magnetic fields, initially concentrated 
on the surface of the bubbles, are expelled when bubbles collide thanks to 
the finite value of the conductivity. 

The idea that charge separation can be generated during first-order phase 
transitions has been exploited in |242j . The suggestion is again that there 
are baryon number gradients at the phase boundaries leading to thermoelec- 
tric terms. In the process of bubble nucleation and collisions turbulence is 
then produced. In spite of the fact that the produced fields are sizable, the 
correlation scale, as previously pointed out, is constrained to be smaller than 
100 pc. 

In a first-order phase transition the phases of the complex order param- 
eter of the nucleated bubbles are not correlated. When the bubbles undergo 
collisions a phase gradient arises leading to a source terms for the evolution 
equation of the gauge fields. Kibble and Vilenkin |243j proposed a gauge- 
invariant difference between the phases of the Higgs field in the two bubbles. 
This idea has been investigated in the context of the Abelian-Higgs model 
|2441 1245t I246j . The collision of two spherical bubbles in the Abelian-Higgs 
model leads to a magnetic field which is localized in the region at the inter- 
section of the two bubbles. The estimate of the strength of the field depends 
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crucially upon the velocity of the bubble wall. The extension of this idea to 
the case of the standard model SU{2) x f/(l) has been discussed in |247j . 
A relevant aspect to be mentioned is that the photon field in the broken 
phase of the electroweak theory should be properly defined. In |248j it has 
been shown that the definition employed in |247j is equivalent to the one 
previously discussed in |249j . 

It has been argued by Vachaspati |25Uj that magnetic fields can be gener- 
ated at the electroweak time even if the phase transition is of second order. 
The observation is that, provided the Higgs field fluctuates, electromagnetic 
fields may be produced since the gradients of the Higgs field appear in the 
definition of the photon field in terms of the hypercharge and SU{2) fields. 
Two of the arguments proposed in |25Uj have been scrutinized in subsequent 
discussions. The first argument is related to the averaging which should be 
performed in order to get to the magnetic field relevant for the MHD seeds. 
Enqvist and Olesen j251j noticed that if line averaging is relevant the ob- 
tained magnetic field is rather strong. However |251j (see also |15Uj ) volume 
averaging is the one relevant for MHD seeds. 

The second point is related to the fact that the discussion of |25Uj was 
performed in terms of gauge-dependent quantities. The problem is then to 
give a gauge-invariant definition of the photon field in terms of the standard 
model fields. As already mentioned this problem has been addressed in |247j 
and the proposed gauge-invariant is equivalent |248j to the one proposed in 

In |252j a mechanism for the generation of magnetic fields at the elec- 
troweak epoch has been proposed in connection with the AMHD equations. 
The idea is to study the conversion of the right-electron chemical potential 
into hypercharge fields. In this context the baryon asymmetry is produced 
at some epoch prior to the electroweak phase. The obtained magnetic fields 
are rather strong (i.e. |-B| ~ 10^^ G at the EW epoch) but over a small scale 
, i.e. 10~^H~J dangerously close to the diffusivity scale. 

The final point to be mentioned is that, probably, the electroweak phase 
transition is neither first order nor second order but it is of even higher or- 
der at least in the context of the minimal standard model |253| 1254] . This 
conclusion has been reached using non-perturbative techniques and the rele- 
vant point, in the present context, is that for Higgs boson masses larger than 
mw the phase transition seems to disappear and it is possible to pass from 
the symmetric to the broken phase without hitting any first or second order 
phase transition. 
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There have been also ideas concerniing the a possible generation of mag- 
netic fields at the time of the QCD phase transition occurring roughly at 
T ~ 140MeV, i.e. at the moment when free quarks combine to form colorless 
hadrons. The mechanism here is always related to the idea of Biermann with 
thermoelectric currents developed at the QCD time. Since the strange quark 
is heavier than the up and down quarks there may be the possibility that the 
quarks develop a net positive charge which is compensated by the electric 
charge in the leptonic sector. Again, invoking the dynamics of a first-order 
phase transition, it is argued that the shocks affect leptons and quarks in a 
different way so that electric currents are developed as the bubble wall moves 
in the quark-gluon plasma. In |255j the magnetic field has been estimated to 
be ~ G at the time of the QCD phase transition and with typical scale 
of the order of the meter at the same epoch. 

In \2bb\ 1257] it has been pointed out that, probably, the magnetic fields 
generated at the time of QCD phase transition may be much stronger than 
the ones estimated in |255j . The authors of |256| 1257] argue that strong 
magnetic fields may be generated when the broken and symmetric phase of 
the theory coexist. The magnetic fields generated at the boundaries between 
quark and hadron phases can be, according to the authors, as large as 10^ G 
over scales of the order of the meter at the time of the QCD phase transition. 

Recently, in a series of papers, Boyanovsky, de Vega and Simionato |258| 
12591 l26Uj studied the generation of large scale magnetic fields during a phase 
transition taking place in the radiation dominated epoch. The setting is 
a theory of charged scalar fields coupled to an Abelian gauge field that 
undergoes a phase transition at a critical temperature much larger than the 
electroweak scale. Using non-equilibrium field theory techniques the authors 
argue that during the scaling regime (when the back-reaction effects are 
dominant) large scale magnetogenesis is possible. The claim is that the 
minimal dynamo requirement of Eq. ()6.1|) is achievable at the electroweak 
scale. Furthermore, much larger magnetic fields can be expected if the scaling 
regime can be extended below the QCD scale. 

6.4 Mixed mechanisms 

From the above discussion it is apparent that the mechanisms producing 
magnetic fields at some moment in the life of the Universe (and inside the 
Hubble radius) have to justify the large correlation scale of the magnetic seeds 
required by subsequent MHD considerations. The possibility of turbulent 
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behaviour in the early Universe, even though physically justified, seems to be 
difficult to quantify even if rather reasonable quantitative estimates appeared 
so far. On the other hand, inffationary scenarios can lead to magnetic seeds 
with large correlation scale. 

A useful point to be stressed is that there may exist "mixed" mecha- 
nisms. Consider the following example which has been proposed in |171j and 
|172j . Suppose that during the inffationary phase a stochastic background of 
gauge ffelds is produced thanks to the breaking of conformal invariance. The 
spectral amphtude of the generated ffeld will be constrained by the critical 
density bound and by all the other dynamical considerations related to the 
finite values of the diffusion scales. 

The generated magnetic ffeld spectrum will be distributed over a large 
interval of frequencies. In particular, the energy spectrum of the magnetic 
ffeld will be non negligible at some intermediate scale, like, for instance, the 
electroweak scale, i.e. the physical frequency Wew ~ -f^ew In other words, 
if magnetic ffelds are generated thanks to some inffationary mechanism, the 
ampliffed magnetic inhomogeneities will be reentering at different epochs 
during the radiation dominated phase and the scales which left the horizon 
later during inffation will reenter earlier. 

Given the background of abehan gauge ffelds, the dynamics of the elec- 
troweak epoch may add interesting physical features. For instance the mag- 
netic ffeld of inffationary origin will be, as discussed, topologicaly trivial: its 
magnetic gyrotropy will vanish. However, during the EW epoch, thanks to a 
pseudoscalar vertex due either to a pseudo goldstone boson or to the chem- 
ical pootential, the topologically trivial hypermagnetic background may be 
turned into a topologically non-trivial background and large magnetic helic- 
ity may be generated. 

7 Effects of primordial magnetic fields 

In the context of various mechanisms, the magnetic ffelds generated at large 
at the scale of the protogalactic collapse may be rather small. In this Section 
the possible effects of magnetic present over much smaller length-scales will 
be analyzed. Two important scales are the Hubble radius at the electroweak 
epoch and at the big-bang nucleosynthesis epoch. 



98 



7.1 Electroweak epoch 



Hypermagnetic fields present for temperatures T > lOOGeV have a twofold 
effect: 

• they may affect the phase diagram of the electroweak theory; 

• they may offer a mechanism to seed the baryon asymmetry of the Uni- 
verse (BAU). 

7.1.1 EW phase diagram 

The physical picture of the possible effects of magnetic fields on the elec- 
troweak phase diagram is exactly the same as the macroscopic description of 
superconductors in the presence of an external magnetic field. The normal- 
superconducting phase transition, being of second order in the absence of 
magnetic fields, becomes of first order if a magnetic field is present. The rea- 
son for this is the Meissner effect, i.e. the fact that the magnetic field cannot 
propagate inside a superconducting cavity, and, therefore, creates an extra 
pressure acting on the normal-superconducting boundary. Our consideration 
below explores this simple picture. 

Consider the plain domain wall that separates broken and symmetric 
phase at some temperature T, in the presence of a uniform hypercharge 
magnetic field Yj. Far from the domain wall, in the symmetric phase, the 
non-Abehan SU(2) field strength {Wf) is equal to zero, because of a non- 
perturbativc mass gap generation. Inside the broken phase, the massive Zj 
combination of Yj and H^?, 

Zj = cos dwW^ - sin OwYj (7.1) 

must be equal to zero, while the massless combination, corresponding to 
photon survives. The matching of the fields on the boundary gives 

AJ^ = y j cos Ow Thus, an extra pressure ^\'Hy\'^ siv? 0^ acts on the do- 
main wall from the symmetric side. At the critical temperature it must be 
compensated by the vacuum pressure of the scalar field. If we neglect loop 
corrections associated with the presence of magnetic fields, then the condition 
that determines the critical temperature is: 

\\HY\^sm^9w = V{Q,T,)-V{^^^,T,) , (7.2) 
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where V{ip, T) is the effective potential in the absence of magnetic field, ipmin 
is the location of the minimum of the potential at temperature T. 

The above consideration was dealing with the uniform magnetic fields. 
Clearly, it remains valid when the typical distance scale of inhomogeneities 
of the magnetic field are larger than the typical bubble size. This is the 
case for bubbles smaller than the magnetic diffusion scale, and, in particular, 
at the onset of the bubble nucleation. Thus, the estimate of the critical 
temperature coming from ()7.2|) is applicable. For bubbles larger than the 
diffusivity scale, the presence of a stochastic magnetic field will considerably 
modify their evolution. In particular, the spherical form of the bubbles is 
very likely to be spoiled. 

These considerations were presented in |lb3j . Later |2(jlj perturbative es- 
timates were performed in order to corroborate the proposed picture. In |262j 
a full non-perturbative analysis of the phhase diagram of the electroweak the- 
ory in the presence of an hyermagnetic background has been performed. As 
previously discussed |253l I254j for values of the Higggs boson mass larger 
than the W boson mass the electroweak phase diagram seems to exhibit a 
cross-over regime. The inclusion of a constant hypermagnetic background 
with typical strength \Hy\/T^ < 0.3 does modify the electroweak phase di- 
agram but does not seem to make the phase transition strongly first order 
for rriH > my/ as expected from perturbative considerations. Furthermore 
for |ffy|/T^ > 0.3 (but still compatible with the critical density bound) a 
new (inhomogeneous) phase has been observed. The analysis performed in 
|.262] included a net hypermagnetic flux. It would be interesting to repeat the 
same calculation in the presence of a non-vanishing hypermagnetic helicity 
(or gyrotropy). 

7.1.2 Baryon asymmetry 

Depending upon the topology of the flux lines, hypermagnetic fields can have 
two distinct effects: 

• topologically non-trivial configurations of the hypermagnetic flux lines 
lead to the formation of hypermagnetic knots whose decay might seed 
the Baryon Asymmetry of the Universe (BAU); 

• even if the electrowak plasma has vanishing hypermagnetic gyrotropy, 
i.e. {Hy ■ V X Hy) = 0, matter-antimatter fluctuations may be gener- 
ated. 
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Consider first the situation where the electroweak plasma contains, for T > 
Tc a network of hypermagnetic knots of the type of the ones described in 
Eqs. fl5.113|) - ()5.116|) . In Section 5 these configurations have been named 
hypermagnetic knots (HK) and they are Chern-Simons condensates carrying 
a non-vanishing (averaged) hypermagnetic helicity. As discussed, HK can 
be dynamically generated from a background of hypermagnetic fields with 
trivial topology provided a (time-dependent) pseudo-scalar is present in the 
plasma (see, for instance, Eq. ()5.120p ). 

In order to seed the BAU a network of HK should be present at high 
temperatures. In fact for temperatures larger than Tc the fermionic number 
is stored both in HK and in real fermions. For T < T^, the HK should release 
real fermions since the ordinary magnetic fields (present after EW symmetry 
breaking) do not carry fermionic number. If the EWPT is strongly first order 
the decay of the HK can offer some seeds for the BAU generation. This 
last condition can be met in the minimal supersymmetric standard model 
(MSSM). 

The integration of the U{1)y anomaly equation gives the CS number 
density carried by the HK which is in turn related to the density of baryonic 
number hb for the case of nj fermionic generations. In fact, using Eqs. 
f)5.109|) - ()5.111|) and after some algebra |163l I172j it can be shown that 

s ^ 27iac s T + Th ry 47r ^ ^ 

{g' is the U{1)y coupling and s = (2/45)7r^5f*T^ is the entropy density; g^,, 
at Tc, is 106.75 in the MSM; Mq = Mp/1.66^). In Eq. T is the 

perturbative rate of the right electron chirality fiip processes (i.e. scattering 
of right electrons with the Higgs and gauge bosons and with the top quarks 
because of their large Yukawa coupling) which are the slowest reactions in 
the plasma and 

783 a'2 \Hy\^ 



H 



22 (T,7r2 T2 



(7.4) 



is the rate of right electron dilution induced by the presence of a hypermag- 
netic field. In the MSM we have that F < Th whereas in the MSSM F can 
naturally be larger than Th- Unfortunately, in the MSM a hypermagnetic 
field can modify the phase diagram of the phase transition but cannot make 
the phase transition strongly first order for large masses of the Higgs boson. 
Therefore, we will concentrate on the case F > F// and we will show that in 
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the opposite limit the BAU will be anyway small even if some (presently un- 
known) mechanism would make the EWPT strongly first order in the MSM. 

It is interesting to notice that, in this scenario, the value of the BAU is 
determined by various particle physics parameters but also by the ratio of 
the hypermagnetic energy density over the energy density sitting in radiation 
during the electroweak epoch. 

Consider now the complementary situation where the electroweak plasma, 
for T > Tc, is filled with topologically trivial hypermagnetic fields. In this 
case, fluctuations in the baryon to entropy ratio will be induced since, in 
any case {{Hy ■ V x HyY) 7^ 0. These fluctuations are of isocurvature type 
and can be related to the spectrum of hypermagnetic fields at the EWPT. 
Defining as 

A(r, Q = ^(5^^) (^ + r,te)), (7.5) 

the fluctuations in the ratio of baryon number density ub to the entropy den- 
sity s att = tc |164j . the value of A(r, tc) can be related to the hypermagnetic 
spectrum which is determined in terms of its amplitude ^ and its slope e. A 
physically realistic situation corresponds to the case in which the Green's 
functions of the magnetic hypercharge fields decay at large distance (i. e. 
e > 0) and this would imply either "blue" ( e > ) or "violet" (e 3> 1) energy 
spectra. The fluctuations of the baryon to entropy ratio generated at the 
electroweak epoch may survive until the BBN epoch |164j . The possibility of 
survival of these fluctuations is related to their typical scale which must ex- 
ceed the neutron diffusion scale appropriately blue-shifted at the electroweak 
epoch. The implications of these fluctuations will be discussed in a moment. 

7.2 Big-bang nucleosynthesis epoch 

Large scale magnetic fields possibly present at the BBN epoch can have an 
impact on the light nuclei formation. By reversing the argument, the success 
of BBN can be used in order to bound the magnetic energy density possibly 
present at the time of formation of light nuclei. Magnetic fields, at the BBN 
epoch may have four distinct effects: 

• they can enhance the rate of the Universe expansion at the correspond- 
ing epoch in a way proportional to pb', 
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• they can affect the reaction rate in a way proportional to aemPe; 

• they can affect the phase space of the electrons; 

• they can induce, prior to the formation of the light nuclei, matter- 
antimatter fluctuations. 

While the first three effects are direct, the fourth effect is mainly caused by 
hypermagnetic fields present even before BBN. 

The idea that magnetic fields may increase the Universe expansion and, 
consequently, affect directly the abundance of ^He (which is directly sensitive 
to the expansion rate) has been pointed out long ago by Greenstein |2ti4j and 
by Matese and O'Connel |255j . From a qualitative point of view the effects 
over the expansion should be leading in comparison with the effects over the 
rate of interactions (which are suppressed by ttem)- 

Thanks to detailed numerical simulations performed indipendently by dif- 
ferent groups (i.e. Cheng et al. |266tl267j : Kernan, Starkman and Vachaspati 
[268"^, I269j : Grasso and Rubinstein \270\ I271j ). the common opinion is that 
the Universe expansion is the leading effect even if this conclusion has been 
reached only recently and not without some disagreements on the details of 
the calculations (see, for instance, j269j ). 

In order to prevent the Universe from expanding too fast at the BBN 
epoch pb < 0.27pi. where p^, is the energy density contributed by the stan- 
dard three light neutrinos with masses much smaller than the MeV. The 
BBN bound is physically relevant in many situations since it is a bit more 
constraining than the critical density bound. 

The BBN bounds discussed so far are derived assuming a stochastic mag- 
netic field at the BBN epoch. This means that the isotropy of the geometry 
is not affected. However, it is not unreasonable to think of the possibility 
that a magnetic field with a preferred direction may break the isotropy of 
the background. Thus bounds on the isotropy of the expansion at the BBN 
epoch may be turned into bounds on the shear induced by the presence of a 
magnetic field (see, for instance, jZ721l2ZS] and P7I1I77H] 1 

Let us now come to the last point. The matter-antimatter fluctuations 
created at the electroweak epoch thanks to the presence of hypermagnetic 
fields may survive down to the epoch of BBN. Since the fluctuations in the 
baryon to entropy ratio are, in general, not positive definite, they will induce 
fluctuations in the baryon to photon ratio, at the BBN epoch. The possible 
effect of matter-antimatter fluctuations on BBN depends on the typical scale 
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(i/r„) 



Figure 8: The value of the baryon number fluctuations for different parame- 
ters of the hj^ermagnetic background ^ and e is reported. 

of of the baryon to entropy ratio at the electroweak epoch. Recalling that for 
T ~ Tc ~ 100 GeV the size of the electrowek horizon is approximately 3 cm, 
fluctuations whose scale is well inside the EW horizon at have dissipated 
by the BBN time through (anti) neutron diffusion. The neutron diffusion 
scale at Tc is 

r,(Tj= 0.3 cm. (7.6) 

The neutron diffusion scale at T = 1 keV is 10^ m, while, today it is 10~^ 
pc, i.e. of the order of the astronomical unit. Matter-antimatter fluctuations 
smaller than 10^ m annihilate before neutrino decoupling and have no effect 
on BBN. Two possibilities can then be envisaged. We could require that 
the matter-antimatter fluctuations (for scales r > rn) are small. This will 
then imply a bound, in the (.^,e) plane on the strength of the hypermagnetic 
background. In Fig. |H1 some typical baryon number fluctuations have been 
reported for different choices of the parameters ^ and e. In Fig. El such an 
exclusion plot is reported with the full line. With the dashed line the bound 
imphed by the increase in the expansion rate (i.e. pn < 0.2p^) is reported. 
Finally with the dot-dashed line the critical density bound is illustrated for 
the same hypermagnetic background. The second possibility is to study 
the effects of large matter-matter domains. These studies led to a slightly 
different scenario of BBN |276j , namely BBN with matter-antimatter regions. 
This analysis has been performed in a series of papers by Rehm and Jedamzik 
I2771I2H21 and by Kurki-Suonio and Sihvola |2IH1 1211 (see also [2101 IMl)- 
The idea is to discuss BBN in the presence of spherically symmetric regions 
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Figure 9: The parameter space of the hypermagnetic background in the case 
A(r, tc) <nB/s for r > rn (full line). 

of anti-matter characterized by their radius ta and by the parameter R, i.e. 
the matter/antimatter ratio. Furthermore, in this scenario the net baryon- 
to-photon ratio, r], is positive definite and non zero. Antimatter domains 
larger than 10^ m at 1 keV may survive until BBN and their dissipation has 
been analyzed in detail |279j . Antimatter domains in the range 

10^ m <rA < 10^ m (7.7) 

at 1 keV annihilate before BBN for temperatures between 70 keV and 1 MeV. 
Since the antineutrons annihilate on neutrons, the neutron to proton ratio 
gets smaller. As a consequence, the ^He abundance gets reduced if compared 
to the standard BBN scenario. The maximal scale of matter-antimatter 
fluctuations is determined by the constraints following from possible distor- 
tions of the CMB spectrum. The largest scale is of the order of 100 pc 
(today), corresponding to 10^^ m at 1 keV. Suppose that matter-antimatter 
regions are present in the range of Eq. (|7.7p . Then the abundance of ^He 
get reduced. The yield of ^He are reported as a function of R, the matter- 
antimatter ratio and r^. Now, we do know that by adding extra- relativistic 
species the '^He can be increased since the Universe expansion gets larger. 
Then the conclusion is that BBN with matter-antimatter domains allows for 
a larger number of extra-relativistic species if compared to the standard BBN 
scenario. This observation may have implications for the upper bounds on 
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Figure 10: From \2Ki\ the ^He yield is illustrated in the {R, ta) plane for 
r] = 6 X 10~^°. As the matter/antimatter ratio decreases, we recover the 
standard ^He yield. 



the stochastic GW backgrounds of cosmological origin |283j since the extra- 
relativistic species present at the BBN epoch can indeed be interpreted as 
relic gravitons. 



8 CMB physics an magnetic fields 

The vacuum fluctuations of the gauge fields present during the inflationary 
stage of expansion may be amplified if conformal invariance is broken. It 
is then plausible that magnetic inhomogeneities are amplified not only at 
the scale of protogalactic collapse but also at larger length-scales, agnetic 
fields can be generated over all physical scales compatible Large-scale mag- 
netic fields may then have various interesting implications on the physics of 
the CMB and of its anisotropies. The following possible effects have been 
discussed through the years: 

• distortion of the Planckian spectrum; 

• shift of the polarization plane (provided the CMB is linearly polarized); 

• shift in the position of the first Doppler peak; 
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• generic increase of the amount of the (primary) anisotropy. 

On top of these effects, magnetic fields can also modify the evolution of the 
tensor fluctuations of the geometry for typical length scales much smaller 
than the ones probed by CMB anisotropy experiments. 

The possible distortions of the Planckian spectrum of CMB are usually 
discussed in terms of a chemical potential which is bounded, by experimental 
data, to me < 9 x 10~^. Magnetic field dissipation at high red-shift implies 
the presence of a chemical potential. Hence bounds on the distorsion of the 
Planckian spectrum of the CMB can be turned into bounds on the magnetic 
field strength at various scales. In particular j287j the obtained bound are 
such that -B < 3 X 10^^ G for comoving coherence lengths between 0.4 kpc 
and 500 kpc. 

Large scale magnetic fields can also afffect the poosition of the Doppler 
peak. In |288j this analysis has been performed in a non-relativistic approxi- 
mation where the scalar perturbations of the geometry obey linearized New- 
tonian equations of motion. It has been found that, in this approximation, 
the effect of the presence of the magnetic fields is an effective renormalization 
of the speed of sound of the baryons. 

8.1 Large-scale magnetic fields and CMB anisotropies 

Large scale magnetic fields present prior to the recombination epoch may 
act as a source term in the evolution equation of the cosmological perturba- 
tions. If magnetic fields are created over large length scales (possibly even 
larger than the Mpc) it is rather plausible that they may be present after 
matter-radiation equality (but before decoupling) when primordial fluctua- 
tions are imprinted on the CMB anisotropies. In fact, the relevant modes 
determining the large scale temperature anisotropies are the ones that are 
outside the horizon for rj^q < r] < r/dcc where rj^q and r^dec are, respectively, 
the equality and the decoupling times. After equality there are two comple- 
mentary approaches which can be used in order to address phenomenological 
implications of large scale magnetic fields: 

• large scale magnetic fields are completely homogeneous outside the hori- 
zon after equality; 

• large scale magnetic fields are inhomogeneous after equality. 
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The idea that homogeneous magnetic fields may affect the CMB anisotropies 
was originally pointed out by Zeldovich 2841 1285j and further scrutinized 
by Grishchuk, Doroshkevich and Novikov [286]. Since this idea implies that 
large-scale magnetic fields are weakly gravitating the details of the discussion 
(and of the possible generalization of this idea) will be postponed to the 
following Section. 

The suggestion that inhomogeneous magnetic fields created during infla- 
tion may affect the CMB anisotropies was originally proposed in |289j . In 
|289j it was noticed that large scale magnetic fields produced during infiation 
may also lead to large fiuctuations for modes which are outside the horizon 
after equality. In |289j it was noticed that either the produced magnetic fields 
may be used to seed directly the CMB anisotropies, or, in a complementary 
persepective, the CMB data can be used to put constraints on the production 
mechanism. 

If large-scale magnetic fields are inhomogeneous, their energy momentum 
tensor will have, in general terms, scalar, vector and tensor modes 

ST,. = 5T(^^ + 5T^r + 5TiI\ (8.1) 

which are decoupled and which act as source terms for the scalar, vector and 
tensor modes of the geometry 

S9,u = SglfJ + 6gZ^ + 6g'^^J, (8.2) 

whose gauge- invariant description is well explored |199| l20U] (see also |201j ). 
The effect of the scalar component of large scale magnetic fields should be 
responsible, according to the suggestion of |289j . of the large scale tempera- 
ture anisotropies. The direct seeding of large-scale temperature anisotropies 
seems unlikely. In fact, the spectrum of large scale gauge fields leading to 
plausible values of the magnetic energy density over the scale of protogalactic 
collapse is typically smaller than the value required by experimental data. 
However, even if the primordial spectrum of magnetic fields would be tai- 
lored in an appropriate way, the initial conditions for the plasma evolution 
compatible with the presence of a sizable (but undercritical) magnetic field 
are of isocurvature type |290j. 

In order to illustrate this point consider the evolution equations for the 
gauge-invariant system of scalar perturbations of the geometry expressed in 
terms of the two Bardeen potential $ and ^ corresponding, in the Newtonian 
gauge, to the fiuctuations of the temporal and spatial component of the 
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metric |2(JU|, I2(J1] . The linearization of the (00), (Oi) and components 
of the Einstein equations leads, respectively, to: 

V^*-3H(H* + *') = ^^ + ^^T„». (8.3) 

1 2 2 2 

— didH^ -m = — '^5%^ + pS, (8.5) 

2 ^ ' 2Ml ' 2Mp ' ^ ^ 

where 5 = 5p/p and w = p/p is the usual barotropic index, and 



5V 



- (8-6) 

are the fluctuations of the energy-momentum tensor of the magnetic field. 
Eqs. ()9.10p - ()8.5p have to be supplemented by the perturbation of the covari- 
ant conservation equations for the fluid sources: 

6' - {l + w)V\-3^'{l + w) = 0, (8.7) 

'ill 

u' + n(l~3w)u 5-$ = 0. (8.8) 

w + I 

If the magnetic field is force-free, i.e. i3 x V x i3 the above system simplifies. 
Notice that Eq. ()8.8p has been already written in the case where the Lorentz 
force is absent. However, the forcing term appears, in the resistive MHD 
approximation, in Eq. (|8.6|) . The term i3 x V x i3 is suppressed by the 
conductivity, however, it can be also set exactly to zero. In this case, vector 
identities imply 

iBxVxB), = B^diB^ - Bmd'^B, = 0, (8.9) 

i.e. BmdiB"^ = Bmd'^Bi. With this identity in mind the system of Eqs. 
fl8.3p - ()8.5p and ()8.7|) - ()8.8p can be written, in Fourier space, as 
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wp6k 



Bl{k) 



Bl{k), 



6'k + A;2(l + w)uk - ?>%{w + 1) = 0, 

w 

u'^ + H{1 - 3w)Uk r5fc - $i 



w + 1 



where 



Blik) 



£pB{p) -Bik-p). 



(8.11) 

(8.12) 

(8.13) 
(8.14) 

(8.15) 



Combining Eqs. (jSTHIl - fjOTjl with Eq. (j8T2|l the following decoupled 
equation can be obtained 



^'fc' + 37^(1 + w)% + k^w^^ 



Bl{k) 



n\i-3w) -2n' 

2fc2 



— w 



(8.16) 



Consider, for simplicity, the dark-matter radiation fluid after equality and 
in the presence of a force-free magnetic field. In this case Eq. ()8.16p can be 
immediately integrated: 



Mk)-^M'--^(^-^ 



5 \ 7] J Snk^M^ V V 



^o{k) 



5 \ 1] J IGnk'^Mp \ r] J 



The solutions given in Eq. ()8.17p determine the source of the density and 
velocity fluctuations in the plasma, i.e. 



5; - -V'u, - W 



0, 



u' 



$ = 0, 

m - 3^' = 0, 
HUr^ - $ = 0, 



(8.18) 

(8.19) 

(8.20) 
(8.21) 



The solution for the velocity fields and density contrasts can be easily ob- 
tained by integrating Eqs. (j8.18|) - (j8.2ip with the source terms determined 
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by Eq. (|8.17|) . The purpose is not to integrate here this system (see for 
instance j291j for the standard case without magnetic fields). Rather it is 
important to notice that there are two quahtatively different situation. The 
first situation is the one where, 

i*»wi«8Sii(T)'- 

namely the case where the constant mode of the Bardeen potential is negli- 
gible if compared to the contribution of the magnetic field. The derivative 
of the Bardeen potential will be, however, non vanishing. In this case the 
CMB anisotropies are said to be seeded by isocurvature initial conditions for 
the fluid of radiation and dark-matter present after equality. In the opposite 
case, nemely 

^o(^) fVeq 



the constant mode of the Bardeen potential is provided, for instance, by 
inflation and the magnetic field represent a further parameter to be taken 
into account in the analysis of CMB anisotropies. Assuming, roughly, that 
|^o(fc)| ~ 10"^, Eq. ()8.2Hj) implies that for the typical scale of the horizon at 
decoupling the critical fraction of magnetic energy density should be smaller 
than 10~^. In the case of Eq. ()8.23|) the leading order relations among the 
different hydrodynamical fluctuations are well known and can be obtained 
from Eqs. ()8.18|l - ()8.21|1 if the magnetic field is approximately negligible, i.e. 

ku, ~ ku^, S, ~ (4/3)5^. (8.24) 

This result has a simple physical interpretation, and implies the adiabaticity 
of the fluid perturbations. The entropy per matter particle is indeed pro- 
portional to S" = T^/rim, where is the number density of matter particles 
and T is the radiation temperature. The associated entropy fluctuation, 6S, 
satisfies 

where we used the fact that pr ~ T'^ and that pm = mrim, where m is the 
typical mass of the particles in the matter fluid. Equation ()8.24|) thus implies 
6S/S = 0, in agreement with the adiabaticity of the fluctuations. 

After having computed the corrections induced by the magnetic field on 
the (leading) adiabatic initial condition the temperature anisotropies can be 
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computed using the Sachs- Wolfe effect |292j . In terms of the gauge- invariant 
variables introduced in the present analysis, the various contributions to the 
Sachs- Wolfe effect, along the n direction, can be written as 



5T 

— {n,r]o,Xo) = 



5.^ 
- + n 



('7dec,^('7dec)) 



'VO 



(8.26) 

where tiq is the present time, and x{ri) = Xq — n{rj — rjo) is the unperturbed 
photon position at the time t] for an observer in xq. The term {Tb is the pecu- 
liar velocity of the baryonic matter component. In Fig. ^2 the Ct are plotted 




Figure 11: The spectrum of is illustrated for a fiducial set of parameters 
(/lo = 0.65, fib = 0.04733, = 0.7, fi^ = 0.25267) and for fiat (full line, 
n = 1), slightly red (dashed hne, n = 0.9) and slightly blue (dotted line, 
n = 1.02) spectral indices for the constant (adiabatic) mode. 

in the case of models with adiabatic initial conditions. The data points re- 
ported in Figs. HU and 1121 are those from COBE ^,^7^, BOOMERANG 
[M], DASI MAXIMA and ARCHEOPS [298J. Notice that the 

data reported in |298j fill the "gap" between the last COBE points and the 
points of 1223 IM IM- In Fig. [HI the recent WMAP data are reported 
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Figure 12: The WMAP data are reported. 



Recall, that the temperature fluctuations are customarily discussed in 
terms of their Legendre transform 

ST 

— (fo, n, r]o) = aemixo)Yemin), (8.27) 

£,Tn 

where the coefficients agm define the angular power spectrum Ci by 

O^im ■ O^l'm' I — ^U'^mm'Ce, (8.28) 



and determine the two-point (scalar) correlation function of the temperature 
fluctuations, namely 

(ST ST \ / \ 

— {n) — {n')\ = J2 («^m4m')^^m(^)^^'m'(^^') 

= -^J2i'^i+l)CePeicos^). (8.29) 
47r ^ 

In Fig. ^2 and ^1 the experimental points are presented in terms of the 
spectra. 

In Fig. ^2 and El the curves fitting the data are obtained by imposing 
adiabatic initial conditions. If magnetic fields would seed directly the CMB 
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anisotropies a characteristic "hump" would appear for £ < 100 which is not 
compatible with the experimental data. If, on the contrary, the condition 
fl8.23|) holds, then the modifications due to a magnetic field can be numeri- 
cally computed for a given magnetic spectral index and for a given amplitude. 
This analysis has been performed by Koh and Lee |302j . The claim is that 
by modifying the amplitude of the magnetic field in a way compatible with 
the cosmological constraint the effect on the scalar can be as large as the 
10% for a given magnetic spectral index. 

The main problem is order to detect large-scale magnetic fields from the 
spectrum of temperature fluctuations is the parameter space. On top of the 
usual parameters common in the CMB analysis at least two new parameters 
should be introduced, namely the magnetic spectral index and the amphtude 
of the magnetic field. The proof of this statement is that there is no avail- 
able analysis of the WMAP data including also the presence of large scale 
magnetic fields in the initial conditions according to the lines presented in 
this Section. 

To set bounds on primordial magnetic fields from CMB anisotropies it 
is often assumed that the magnetic field is fully homogeneous. In this case 
the magnetic field amplitude has been bounded to be 5 < 10~^ G at the 
decoupling time |303j (see also j305j for a complementary analysis always in 
the case of a homogeneous magnetic field). 

There is the hope that since magnetic fields contribute not only to the 
scalar fluctuations, but also to the vector and tensor modes of the geometry 
useful informations cann be also obtained from the analysis of vector and 
tensor power spectra. In |304j temperature and polarization power spectra 
induced by vector and tensor perturbations have been computed by assuning 
a power-law spectrum of magnetic inhomogeneities. In [30B[ I307| I308j the 
tensor modes of the geometry have been specifically investigated. In j309| 
13101 1311] it is argued that magnetic fields can induce anisotropies in the 
polarization for rather small length scales, i.e. C. > 1000. 

As already pointed out in the present Section, magnetic fields may be 
assumed to be force free in various cases. This is an approximation which 
may also be relaxed. However, even if magnetic fields are assumed to be 
force-free there is no reason to assume that their associated magnetic helicity 
(or gyrotropy) vanishes. In |312j the posssible effects of helical magnetic 
fields on CMB physics have been investigated. If helical flows are present at 
recombination, they would produce parity violating temperature-polarization 
correlations. The magnitude of helical flows induced by helical magnetic 
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fields turns out to be unobservably small, but better prospects of constraining 
helical magnetic fields come from maps of Faraday rotation measurements of 
the CMB ^312] . 

8.2 Faraday rotation of CMB 

Large scale magnetic fields present at the decoupling epoch can also depo- 
larize CMB ^313j. The polarization of the CMB represents a very interesting 
observable which has been extensively investigated in the past both from the 
theoretical |314j and experimental points of view |315t I316j . Forthcoming 
satellite missions like PLANCK |317j seem to be able to achieve a level of 
sensitivity which will enrich decisively our experimental knowledge of the 
CMB polarization with new direct measurements. 

If the background geometry of the universe is homogeneous but not 
isotropic the CMB is naturally polarized |314j . This phenomenon occurs, for 
example, in Bianchi-type I models |318j . On the other hand if the background 
geometry is homogeneous and isotropic (like in the Friedmann-Robertson- 
Walker case) it seems very reasonable that the CMB acquires a small de- 
gree of linear polarization provided the radiation field has a non-vanishing 
quadrupole component at the moment of last scattering |319| l32Uj . 

Before decoupling photons, baryons and electrons form a unique fluid 
which possesses only monopole and dipole moments, but not quadrupole. 
Needless to say, in a homogeneous and isotropic model of FRW type a pos- 
sible source of linear polarization for the CMB becomes efficient only at the 
decoupling and therefore a small degree of linear polarization seems a firmly 
established theoretical option which will be (hopefully) subjected to direct 
tests in the near future. The linear polarization of the CMB is a very promis- 
ing laboratory in order to directly probe the speculated existence of a large 
scale magnetic field (coherent over the horizon size at the decoupling) which 
might actually rotate (through the Faraday effect) the polarization plane of 
the CMB. 

Consider, for instance, a linearly polarized electromagnetic wave of phys- 
ical frequency uj travelling along the x direction in a plasma of ions and 
electrons together with a magnetic field {B) oriented along an arbitrary di- 
rection ( which might coincide with x in the simplest case). If we let the 
polarization vector at the origin (a; = ?/ = 2; = 0,t = 0)be directed along the 
y axis, after the wave has traveled a length Ax, the corresponding angular 
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shift (Aa) in the polarization plane will be 



2m \ LU J 



(8.30) 



From Eq. ()8.30|) by stochastically averaging over all the possible orientations 
of B and by assuming that the last scattering surface is infinitely thin (i.e. 
that AxfeUe — (7t^ where a-r is the Thompson cross section) we get an 
expression connecting the RMS of the rotation angle to the magnitude of B 

at t — trfec 



(in the previous equation we imphcitly assumed that the frequency of the 
incident electro-magnetic radiation is centered around the maximum of the 
CMB). We can easily argue from Eq. ()8.3H) that if B(tdec) > Be the expected 
rotation in the polarization plane of the CMB is non negligible. Even if we 
are not interested, at this level, in a precise estimate of Aa, we point out that 
more refined determinations of the expected Faraday rotation signal (for an 
incident frequency um ~ 30 GHz) were recently carried out |321|,l322j leading 
to a result fairly consistent with ()8.30|) . 

Then, the statement is the following. // the CMB is linearly polarized 
and if a large scale magnetic field is present at the decoupling epoch, then 
the polarization plane of the CMB can be rotated |313j . The predictions of 
different models can then be confronted with the requirements coming from 
a possible detection of depolarization of the CMB |313j . 

8.3 Relic gravitational waves 

Magnetic fields can source the evolution equations of the fluctuations of the 
geometry also over length-scales much smaller than the ones where CMB 
anisotropy experiments are conduced. This suggests that magnetic inhomo- 
geneities may leave an imprint on the relc background of gravitational waves. 

If a hypermagnetic background is present for T > Tc, then, as discussed in 
Section 5 and 6, the energy momentum tensor will acquire a small anisotropic 
component which will source the evolution equation of the tensor fluctuations 
of the metric. Suppose now, that |7iy| has constant amplitude and that it is 




Be = 10-^ G, u;a/ ^ 3 X 10^° Hz 



(8.31) 
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also homogeneous. Then as argued in |323j we can easily deduce the critical 
fraction of energy density present today in relic gravitons of EW origin 

n^{to) = ^c^z;y, p,{T,)^K^T'^ (8.32) 

Pc 

{ze(^ is the redshift from the time of matter-radiation, equality). Because of 
the structure of the AMHD equations, stable hypermagnetic fields will be 
present not only for uj^^ ~ k^^/a but for all the range cUew < oj < where 
oj^ is the diffusivity frequency. Let us assume, for instance, that ~ 100 
GeV and g^, = 106.75. Then, the (present) values of cUew is 

.™(W-2.01xlO-'(^)(^)'V (8.33) 

Thus, co'o-(to) ~ lO^t^ew Suppose now that Tc ~ 100 GeV; than we will have 
that Uewito) ~ 10"^ Hz. Suppose now, that 

IHyI/T^ > 0.3, (8.34) 

as, for instance, implied by the analysis of the electroweak phase diagram in 
the presence of a magnetized background. This requirement imposes r ~ 0.1- 
0.001 and, consequently, 

^o^GW - 10"^ - 10"^ (8.35) 

Notice that this signal would occurr in a (present) frequency range between 
10^^ and 10'^ Hz. This signal satisfies the presently available phenomenolog- 
ical bounds on the graviton backgrounds of primordial origin. The pulsar 
timing bound ( which applies for present frequencies ujp ~ 10~^ Hz and 
implies /iq^gw ^ 10~^) is automatically satisfied since our hypermagnetic 
background is defined for 10~^Hz < uj < lO^Hz. The large scale bounds 
would imply /igf^cw < 7 x 10^^^ but a at much lower frequency (i.e. 10~^^ 
Hz). The signal discussed here is completely absent for frequencies u < cUew 
Notice that this signal is clearly distinguishable from other stochastic back- 
grounds occurring at much higher frequencies (GHz region) like the ones 
predicted by quintessential inflation |324l l525. 326J. It is equally distinguish- 
able from signals due to pre-big-bang cosmology (mainly in the window of 
ground based interferometers 327J). The frequency of operation of the in- 
terferometric devices (VIRGO/LIGO) is located between few Hz and 10 kHz 
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\6'27\ . The frequency of operation of LISA is well below the Hz (i.e. 10~^Hz, 
approximately). In this model the signal can be located both in the LISA 
window and in the VIRGO/LIGO window due to the hierarchy between the 
hypermagnetic diffusivity scale and the horizon scale at the phase transition. 
In Fig. ^1 the full thick line illustrates the spectrum of relic gravitational 




Figure 13: The stochastic background of GW produced by inflationary mod- 
els with fiat logarithmic energy spectrum, illustrated together with the GW 
background of hypermagnetic origin. The frequencies marked with dashed 
lines correspond to the electroweak frequency and to the hypermagnetic dif- 
fusivity frequency. 

waves produced in a conventional model for the evolution of the universe. 
The flat plateau corresponds to modes which left the horizon during the in- 
flationary stage of expansion and re-entered duriing the radiation dominated 
phase. The decreasing slope between 10~^^ and 10~^* Hz is due to modes 
leaving the horizon during inflation and re-entering during the maatter dom- 
inated stage of expansion. Clearly, the signal provided by a background of 
hypermagnetic fields can be even 7 order of magnitude larger than the infla- 
tionary prediction. The interplay between gravitational waves and large-scale 
magnetic fields has been also the subject of recent interesting investigations 
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9 Gravitating magnetic fields 



Up to now large scale magnetic fields have been treated in different frame- 
works but always within a perturbative approach. However, magnetic fields 
can also break explicitely the isotropy (but not the homogeneity) of the back- 
ground geometry 

Suppose that at some time ti the Universe becomes transparent to radi- 
ation and suppose that, at the same time, the four-dimensional background 
geometry (which we assume, for simplicity, spatially flat) has a very tiny 
amount of anisotropy in the scale factors associated with different spatial 
directions, namely 

ds^ = df - a\t)dx^ - h\t)[dy^ + dz% (9.1) 

where 5(t), as it will be clear in a moment, has to be only slightly different 
from a{t). The electromagnetic radiation propagating along the x and y axes 
will have a different temperature, namely 

T^it) = Ti^ = Tie-Z^W'^*, Ty{t) = Ti^ = Tie'/^W'^*, (9.2) 
a u 

where H{t) and F{t) are the two Hubble factors associated, respectively 
with a{t) and b{t). The temperature anisotropy associated with the electro- 
magnetic radiation propagating along two directions with different expansion 
rates can be roughly estimated, in the limit H{t) — F{t) <C 1, as 

^ ~ I [H{t) - F{t)]dt = \j s{t)d\ogt (9.3) 

where, in the second equality, we assumed that the deviations from the ra- 
diation dominated expansion can be written as F{t) ~ (1 — e{t))/2t with 
\s{t)\ -C 1. The function e can be connected with the shear parameter, i.e. 
[H(t) — F(t)]/[H(t) + 2F{t)] which measures the anisotropy in the expansion. 
In the standard context |284| 1285] Einstein equations are solved in the metric 
1)9.11) with fluid sources and in the presence of a magnetic fleld. The shear 
parameter is then connected with the magnetic fleld intensity. 

The dynamical origin of the primordial anisotropy in the expansion can 
be connected with the existence of a primordial magnetic fleld (not dynam- 
ically generated but postulated as an initial condition) or with some other 
sources of anisotropic pressure and, therefore, the possible bounds on the 
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temperature anisotropy can be translated into bounds on the time evolution 
of the anisotropic scale factors |3U3j . 

Today the amount of anisotropy in the expansion must be very small 
because of the effect we just described. In the present Section some highly 
speculative considerations on the possible role of magnetic fields in the early 
Universe will be intrtoduced. In spite of the fact that magnetic fields must 
be sufficiently small today, in the past history of the Universe they might 
have been very large even modifying the dynamics of the geometry. During 
the magnetic phase the anisotropy in the expansion will be constant (or 
even grow). However, after the magnetic phase a sufficienly long radiation 
dominated phase may isotropize the background leading to a tiny amount of 
anisotropy in the expansion. This is the basic scenario invoked long ago by 
Zeldovich [2H1I2H3- 

In modern approaches to cosmology, the dynamics of the Universe at very 
high densities is often discussed in terms of the low-energy string effective 
action. It is then useful to analyze the Zeldovich proposal in terms of the 
modern perspective. 

In the low energy limit, the dilaton field is directly coupled to the kinetic 
term of the Maxwell field |332|| 

1 „ 1 



(9.4) 

where F^p = V[q,A^] is the Maxwell field strength and is the covariant 
derivative with respect to the String frame metric g^j^^. Notice that H^^a is the 
antisymmetric field strength. In Eq. (j9.4p the ellipses stand for a possible 
(non-perturbative) dilaton potential and for the string tension corrections 
parameterized by a' = A^. In Eq. ()9.4|) Ffj_^ can be thought as the Maxwell 
field associated with a U{1) subgroup of Eg x Eg. 

Consider a spatially flat background configuration with vanishing an- 
tysimmetric field strangth {H^ya = 0) and vanishing dilaton potential. The 
dilaton depends only on time and the metric will be taken fully anisotropic 
since we want to study possible solutions with a homogeneous magnetic field 
which is expected to break the isotropy of the background: 

ds^ = g^ydx^dx" = dt^ - a^{t)dx^ - h'^{t)dy'^ - c^{t)dz^. (9.5) 

By varying the effective action with respect to y?, g^y and the vector 
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potential we get, respectively 



R'f, S'f.R — — 

/I 2 2 



1 



1 ... „,^„ „ ... 1 



^ y g 



0. 



(9.6) 

(9.7) 
(9.8) 



where denotes covariant differentiation with respect to the metric of Eq. 
fFOI) . Inserting Eq. (|K3H|l into Eq. dHUI) we get that Eq. dHZI) can be 
expressed as 

0. (9.9) 



R; + V^W + ^F,^F 



Consider now a homogeneous magnetic field directed along the x axis. The 
generalized Maxwell equations ()9.8|) and the associated Bianchi identities 
can then be trivially solved by the field strength Fy^ = —F^y. The resulting 
system of equations ()5.58|) . ()9.9|) can than be written, in the metric of Eq. 
dnUI), as 



= + F^ + + H + F + G, 
Hif = HF + HG + H^ + H, 

F^ = HF + FG + F^ + F 



Gip = GH + FG + G^ + G 



262 c2 ' 



(9.10) 
(9.11) 

(9.12) 
(9.13) 



where Eqs. (I9.1(jj) - (j9.13|) correspond, respectively to the (00) and (ii) com- 
ponents of Eq. ()9.9|) . 

The solution of Eqs. (Pnn|l - (|m|l 



^{v) = ^0 + {T-Co + Ao)r] + log 



+ e 



2Aov 



(9.14) 
(9.15) 
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where 



— 



2^2 



(9.16) 



The solution given in Eq. ()9.16|) s given in terms of a new variable (a gener- 
alized conformal time) defined by the following differential relation 



These solutions represent the generalization of pre-big bang solutions 
|123j to the case where a constant magnetic field is included. Generaliza- 
tions of the solution ()9.16p can be obtained in various models anisotropic 
(but homogeneous) models of the Bianchi class |887j like Bianchi 1, 11, 111, 
Vl_i and Vllo according to the usual classification which can be found, for 
instance, in [333 . Further generalizations of the solution ()9.1(jj) are possible 
in the Kant owski- Sachs metric with a magnetic field aligned alond the radial 
coordinate [335. 1336] . 

Magnetized cosmological solutions can be used in order to describe some 
early phase in the life of the universe. This highly speculative idea may 
be implemented in semi-realistic models provided the large anisotropy in 
the expansion decays at later epochs. In the standard context |284| 1285] 
the only source of isotropization is the evolution in a radiation dominated 
background since the shear parameter decays if the background is dominated 
by radiation. If quadratic corrections are present, an initially anisotropic 
solution becomes isotropic as a consequence of the dynamical properties of 
the evolution equations as discussed in detail in |332j . 

Gauge fields have also been studied as a source of infiationary evolution 
by Ford |338j and later by Lidsey |339t I340j . The gauge kinetic term is 
supplemented by a potential which is a function of A^A^. Both exponential 
infiation |338j and power-law infiation |339j can be realized in these models. 



The topics reported in the present review suggest various open problems 
on the nature of large-scale magnetic fields in the present and in the early 
Universe: 

• in principle the morphology of magnetic fields in spiral galaxies could 
be used in order to understand the origin of galactic magnetism but in 
practice the observations offer answers which are often contradictory; 



e "^drj = dt, ip = (p — log abc. 



(9.17) 



10 Concluding Remarks 
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• recent observations of magnetic fields inside Abell clusters confirm the 
presence of strong magnetic fields in the /i G range; 

• there is evidence that also superclusters are magnetized; 

• in spite of the fact that magnetic fields may modify the patterns of 
the CMB anisotropics and, eventually, induce anisotropics in the po- 
larization, the theoretical analysis has been rarely corroborated by the 
comparison with the available experimental data; 

• Faraday rotation of the CMB polarization would be a powerful tool 
to study the possible exsitence of large scale magnetic fields prior to 
recombination; 

• the imprints of magnetic fields on the relic background of gravitational 
waves (in the window of wide band interferometers) allows to test the 
existence of a background of Abehan gauge fields at the electroweak 
epoch. 

It would be important, for the theorist, a reasonable accuracy (within one or- 
der of magnitude) in the experimental determination of large-scale magnetic 
fields. For instance even one order of magnitude in the intensity of the intra- 
cluster magnetic field (i.e. 10~^ G rather than 10~^G) makes a difference as 
far as the problem of the origin is concerned. In the same perspective, more 
accurate determinations of the typical correlation scales of the intra-cluster 
fields would be desirable. 

There is, at the moment, no compelling reason why large-scale magnetic 
fields should not be primordial, at least to some degree: 

• in the early Universe magnetic fields are easily produced during phase 
transitions but their typical correlation scales are small; 

• inflation greatly helps in producing correlations over very large scales 
and the (early) variation of gauge couplings during the inflationary 
phase allows the generation of intense large-scale fields; 

• primordial magnetic fields may have an impact on the thermodynam- 
ical history of the Universe, but, in practice, the obtained constraints 
improve only marginally, in various interesting cases, the critical den- 
sity bound imposed on the magnetic energy density; 
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• if magnetic fields are generated during phase transitions, a growth in 
their correlation scale cannot be excluded because of turbulent phe- 
nomena whose existence can be reasonably expected (but not firmly 
predicted) on the basis of our terrestrial knowledge of magnetized plas- 
mas. 

In the next decade a progress of the empirical knowledge is expected in ap- 
parently unrelated areas like x and 7-ray astronomy, radio-astronomy, CMB 
physics, detection of relic gravitational waves. All these areas are connected 
with the existence and with the properties of large-scale magnetic fields in the 
early and in the present Universe. This connection is sometimes rather frag- 
ile since it is mediated by various theoretical assumptions. However, there is 
the hope that, in a not too distant future, some of the puzzles related to the 
origin and existence of large-scale magnetic fields may be resolved and some 
of the current theoretical guessworks firmly ruled out. 
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A Complements on MHD description 



Two fluids and one-fluid MHD equations 

In a two-fluid plasma description the charge carriers are the ions (for sim- 
plicity we can think of them as protons) and the electrons. The two fluid 
equations treat the ions and the electrons as two conducting fluids which are 
coupled as in Eqs. (j4.8|) - (j4.12|) . Given the two- fluid description, one- fluid 
variables can be defined directly in terms of the two-fluid variables 



Pm(f , t) = rUerieix, t) + TJlpHp^X, t) 

Pq(x, t) = e[np{x, t) - ne{x, t)], 

J 



e[npVp - HeVe] 



In the case of a globally neutral plasma Ug 
and electron equations then become 

d 



~ rir 



(A.l) 

Ug and pg = 0. The ion 



rripUq 



d_ 

dt 



eng{E + Vp X B) - Vpp + Cpe, (A.2) 
-eng{E + VeX B) -Vpe + Cep, (A.3) 



where Cpe and Cpe are the collision terms. In the globally neutral case the 
center of mass velocity becomes 



m„Vri + nieVe 

V = ^ , 

rup + rUe 

and the one-fluid mass and charge density conservations become 

^ + V . (p.,?) = 0, 



(A.4) 



(A.5) 
(A.6) 



Summing up Eqs. ()A.2|) and ()A.3|) leads, with some algebra involving the 
continuity equation, to the momentum transport equation in the one-fluid 
theory: 

'Qy -,1 -, -, 

v = J X B -VP, (A.7) 



Pn 
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where P = Pg + Pi. In Eq. flA.7|) the colhsion term vanishes if there are 
no neutral particles, i.e. if the plasma is fully ionized. The final equation of 
the one-fluid description is obtained by taking the difference of Eqs. ()A.2|1 
and ()A.3|) after having multiplied Eq. ()A.2|) by rrie and Eq. ()A.3|) by nip. 
This procedure is more tricky and it is discussed in standard textbooks of 
plasma physicsfSHl EE]- The key points in the derivation are that the limit 
for rrLe/mp —>■ must be taken. The problem with this procedure is that the 
subtraction of the two mentioned equations does not guarantee that viscous 
and collisional effects are negligible. The result of this procedure is the so- 
called one-fluid generalized Ohm law: 

E + vxB = -J + —{JxB-VPe). (A.8) 
a eUq 

The term J x B is nothing but the Hall current and VPe is often called 
thermoelectric term. Finally the term J /a is the resistivity term and a is 
the conductivity of the one-fluid description. In Eq. ()A.8|) the pressure has 
been taken to be isotropic. This is, however, not a direct consequence of the 
calculation presented in this Appendix but it is an assumption which may 
(and should) be relaxed in some cases. In the plasma physics literature ^1 
15^ the anisotropic pressure contribution is neglected for the simple reason 
that experiments terrestrial plasmas show that this terms is often negligible. 



Conservation laws in resistive MHD 

Consider and arbitrary closed surface S which moves with the plasma. Then, 
by definition of the bulk velocity of the plasma {v we can also write = 
V X dl drj. The (total) time derivative of the flux can therefore be expressed 
as 

B-dt= —■dt+ Bxv-dt, (A.9) 



dt Jy. Jt. dt JdT, 

where 9S is the boundary of S. Using now the properties of the vector 
products (i.e. B x v ■ dl = —v x B ■ dl) we can express v x B though the 
Ohm law given in Eq. ()4.23j) and we obtain that 

vx B ^ -E + -^V y< B (A.IO) 
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Using now Eq. ()A.9|) together with the Stokes theorem, the following expres- 
sion can be obtained 



dt Jt. 



5-dS=/ — + Vx^-rfS /VxVx5-dS. (A.ll) 

Jy,V ot \ Atxg Je 



From the Maxwell's equations the first part at the right hand side of Eq. 
(jA.llll is zero and Eq. ()4.24|1 . expressing the Alfven theorem, is recovered. 

With similar algebraic manipulations (involving the use of various vector 
identities), the conservation of the magnetic helicity can be displayed. Con- 
sider a closed volume in the plasma, then we can write that dV = d^x = 
v± ■ dT, di] = n ■ v± rfS drj where n is the unit vector normal to S (the 
boundary of V, i.e. S = dV) and v± is the component of the bulk velocity 
orthogonal to dV. The (total) time derivative of the magnetic helicity can 
now be written as 

—nM= I d^x—[A-B] + [ A-Bv^-ndT.. (A.12) 
drj Jv orj Jdv=T, 

The partial derivative at the right hand side of Eq. ()A.12j) can be made 
explicit. Then the one-fiuid MHD equations should be used recalling that 
the relation between the electromagnetic fields and the vector potential, for 
instance in the Coulomb gauge. Finally using again the Ohm law and trans- 
forming the obtained surface integrals into volume integrals (through the 
divergence theorem) Eq. ()4.27|1 . expressing the conservation of the helicity, 
can be obtained. 

In spite of the fact that the conservation of the magnetic helicity can be 
derived in a specific gauge, the magnetic helicity is indeed a gauge invariant 
quantity. Consider a gauge transformation 

A^A + V^, (A.13) 

then the magnetic helicity changes as 

d^xA -B ^ d^xA -3 + d^xV ■ [^B] (A. 14) 

(in the second term at the right hand side we used the fact that the magnetic 
field is divergence free) . By now using the divergence theorem we can express 
the volume integral as 

/ d^xV-[^B]= [ iB-nd'L. (A.15) 

Jv JdV=T. 
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— * 

Now if, as we required, B ■ ft — in dV, the integral is exactly zero and Hm 
is gauge invariant. The condition B ■ n = is not specific of a particular 
profile of the magnetic field. It can be always achieved by slicing the volume 
of integration in small flux tubes where, by deflnition, the magnetic field is 
orthogonal to the walls of the flux tube. 
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